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ABSTRACT 


Llne-integrated  refractlvity  ia  isolateO  as  the  basic  atmospheric  parameter 
to  be  determined  ^rrectlona  to  radar  ra;iging  ilata.  A  basis  for  the  (ieter- 
mination  of  Integrated  refractivity  from  measurements  of  the  absorption  of  infra¬ 
red  and  of  micro’A’ave  radiation  ty  atmos|Aeric  oxygen  and  water  vapor  is  then 
established.  Three  absorption-measuring  systems  '»re  next  described,  and  engi¬ 
neering-model  design  sjKx  iflcations  are  presented.  Each  of  the  specifications 
offers  particular  advantages  to  particular  ranging  operatioris.  It  is  estimated 
*^hat  incorporation  of  such  devices  into  ranging  systems  will  permit  real-time 
measurements  of  slant  range  and  elev^ation  angle  to  accuracies  five  to  ten  times 
tl.o8e  currently  achieved. 
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EVALUATION  -  CONTRACT  AF30(602)-3715 


1.  To  apply  corrections  for  refractive  errors  encountered  by 
radar  systems,  an  accurate  account  of  the  atmosphere  is 
essential.  Present  day  measuring  techniques,  such  as  radio¬ 
sondes  and  refractometers,  provide  only  the  loi'g  term  vertical 
refractive  structure  of  the  atmosphere",  thus  a  need  exists  for 
a  real  time  sampling  technique  capable  of  measuring  the  8ho.’'t 
term  fluctuations  of  the  refractive  index  along  a  ray  path. 
These  fluctuations  can  produce  large  refraction  errors,  espe¬ 
cially  at  low  elevation  angles. 

2.  This  contract  was  awarded  to  the  Boeing  Company  to  study 
techniques  that  will  provide  a  real  time  measut 'ment  of  the 
integrated  refractivity  along  a  ray  path.  The  conclusions 
of  this  study  showed  that  the  most  accurate  technique  for 
measuring  the  refractivity  along  a  radar  path  would  be  a 
radiation  transfer  measurement  technique  (measurement  of 
absorbed  radiation  intensity)  performed  in  the  oxygen  and 
water  vapor  absorption  bands  at  the  near  infrared  and  milli¬ 
meter  frequencies.  Also  from  this  study,  two  basic  measure¬ 
ment  systems  have  been  envisioned  to  take  into  account 
cooperative  and  uncooperative  targets. 

3.  Future  work  will  be  to  develop  and  implement  an  instrument 
package  as  dictated  by  the  results  of  this  study.  The  instru¬ 
mentation  will  then  be  used  to  collect  data  under  different 
climatic  conditions  and  correlated  with  the  meteorological 
conditions  of  the  environment  for  evaluation  as  a  real  time 
refraction  correction  system. 

DONALD  S.  LUCZAie^ 

Contract  Engineer 
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INTRODUaiON 


When  usln^  radio  waves  In  long-range  target  detection  and  tracking,  the 
limiting  factor  of  system  performance  and  capability  is  largely  dependent  on  the 
medium  (or  atmosphere)  throng  adiioh  the  waves  are  propagated.  In  the  atmos¬ 
phere,  the  homoepbere  and  the  ionosphere  are  the  regions  that  greatly  influence 
and  affect  radio-wave  propagation.  The  bomosphere  is  defined  as  that  region  of 
the  atnosphere  that  exhibits  siijstantially  uniform  composition  (i.e.,  constant 
mean  molecular  weight)  from  the  surface  of  the  Earth  upwards.  It  includes  both 
the  tioposi^eric  and  stratospheric  regions  of  the  atmosphere.  The  ionosphere 
is  defined  ao  the  region  of  the  atmosphere  characterized  by  elect r  n  densities 
sufficiently  lai  je  to  affect  radio  commumcation.  This  report  is  concerned  only 
with  effects  on  radar  system  performance  ascribable  to  the  molecular  composition 
of  the  atmosphere,  and  thus  relates  only  to  the  homospheric  effects.  Ionospheric 
effects  are  left  to  future  study. 
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Section  I 

MEASUREMENT  PRINCIPLES 


1.  BACK^.ROUND 

In  a  nonrefractlve  environment,  single- station  radar  ranging  and  tracking 
data  admit  to  very  simple  interpretations.  Target  slant  range  is  determined 
simply  as  the  product  of  the  speed  of  light,  in  vacuo,  by  the  target- to- receiver 
ray  transit  or  delay  time.  Target  bearing  is  determined  by  the  angle  of  arrival 
of  the  incoming  ranging  rays. 

Even  in  a  refractive  medium,  such  as  the  atmosphere,  the  above  concepts 
are  used  to  provide  first-order  ai^roximations  to  target  range  and  bearing. 
However,  two  effects  cause  errors  in  both  the  above  mentioned  range  param¬ 
eters.  First,  the  presence  of  the  material  medium  results  in  a  slowing  down 
or  retardation  of  the  ranging  rays;  i.e. ,  the  ranging  signals  travel  at  speeds 
somewhat  less  than  the  speed  of  light  in  vacuo.  The  product  c  t  (where  c  ^ 
is  the  vacuum  velocity  of  light,  and  At  is  the  measured  delay  time)  results  in 
an  overestimate  of  slant  range.  Second,  because  refractive  index  decreases 
with  altitude  in  the  atmosphere,  a  ray  bending  associated  with  the  refract <ve 
index  gradient  occurs,  with  the  result  that  the  ray  angles  of  arrival  are  not 
those  of  the  straight- ray  path.  The  result  is  that  the  target  is  "seen"  in  the 
wrong  direction;  vis. ,  in  the  dlroctioo  of  the  tangent  (at  the  radar  receiver)  to 
the  refracted  ray.  The  effects  of  both  ray  retardation  and  ray  bending  or 
refraction  on  the  radar  measurement  of  range  parameters  are  illustrated  in 
Figure  1.  In  the  figure,  S  — the  electrical  range  —  is  the  first-order  slant- 
range  estimate  given  by  the  previously  defined  product  c^'At.  The  true  slant 
range,  S,  is  closely  approximated*  by  the  product  c  -At,  where  c  denotes  the 
rctual  raj’  propagation  speed  in  the  medium.  The  slant- range  error,  AS,  is 

*  The  approximation  neglects  the  purely  geom«*tric  effect  on  At  attributable  to 
the  finite  differences  in  curved  versus  straight- ray  path  lengths.  For  eleva¬ 
tion  angles  s>3*,  such  effects  are  negligible. 
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Apparent  Target  P..-«ition 


D  RETARDATION  EFFECTS  ON  RANGE  PARAMETERS 


thus  definable  as  uie  difference  between  the  electrics)  and  true  slant  ranges;  the 
elevation  angle  error,  is  deQnable  as  the  difference  between  the  straight- 
rky  direction,  €  ,  and  that  of  the  appc-ent  ray,  .  That  is. 


L3  -3^  -  S, 

(1) 

^  6^  -e. 

(2) 

Methods  of  estimating  the  magnitudes  of  the  error  increments  defined  by 
Equations  1  and  2  are  needed.  Historically,  the  refractive  bending  problem  baa 
received  the  g^atest  amount  of  attention  for  the  longest  period  of  time.  Past 
work  by  various  researchers  (for  example,  those  cited  in  References  1,  2,  3, 
and  4)  has  shown  that  it  is  possible  to  estimate  the  atmospheric  refraction  of 
radio  waves  to  a  respectable  precision.  The  technique  is  to  apply  Snell's  law 
of  refraction  to  effect  a  ray  tracing  through  an  atmosphere  assumed  homogeneous 
in  spherical  shells  surrounding  the  Earth,  and  in  which  refractive  index  decreases 
more  or  less  raonotonically  (and  very  nearly  exponentially)  with  bei^t.  Only 
since  the  advent  of  increasingly  sophisticated  techniques  for  measuring  delay 
times  to  higher  degrees  of  precision  have  ray- retardation  efiec's  assumed  pro¬ 
portions  significant  en^xigh  to  warrant  a  similrr  amount  o'*  attention.  Concur¬ 
rently  arith  the  increased  precision  of  time-delay  data,  however,  has  come  a 
need  for  matching  increases  la  specifying  bearing- angle  data  to  higher  degrees 
of  accura  V.  The  point  has  been  resebod  in  vhin  latter  demand  where,  even  dis¬ 
regarding  tor  the  moment  refractive  bending  'elevation  errors,  troublesome 
errors  in  the  measure  of  elevation  angles  arise  simply  because  of  feed-hom 
droop  or  dish  distortion  in  radars  such  as  the  AN/FPQ-S  pulse  systems  (Refer¬ 
ence  5).  These  purely  mechanical  pixiblems  are  compounded  in  the  tracking 
operation  when  target  velocity  determioations  must  depend  on  mechanical  changes 
in  antenna  azimuth  and  elevation.  The  i-oault  of  this  t  nigma  has  been  the  evolu¬ 
tion  of  ranging  techniques  wherein  the  direct  measure  of  target  bearing  angle  is 
eliminated,  which  is  accompUabed  by  trilateration  aol  ions  that  permit  the 
expression  of  bearing  angle  in  terms  of  range  and  range  difference  data  obtained 
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from  two  or  more  stations.  The  trtlatcration  mode  of  operation  with  pulse 
radar  systems  is  epitomized  at  the  Atlantic  Missile  Range  by  the  GLOTRAC 
pulse-system  tracking  network  (Reference  5).  Equivalent  results  are  obtained 
more  directly  and  more  simply  .dth  c-w  interferometer  systems  such  as  AZUSA 
and  MISTRAM  (Reference  5),  where  target  range  is  measured  at  a  "mother" 
station  by  comparing  the  phase  of  the  signal  received  from  a  beacon  transponder 
with  an  internal  phase  reference.  Range  differences  are  obtained  directly  by 
con^>arlng  idiases  of  received  signals  at  the  mother  and  at  "slave"  stations. 
These  I'ange  differences  define  receiver- to-Urget  direction  cosines.  These 
data,  in  conjunction  with  range  data,  aie  sufficient  to  determine  target  position. 

The  point  to  be  emphasized  is  the  elimination  of  the  direct  measure  of  target 
bearing  angles  in  precision  ranging  operations,  and  the  consequent  emergence  of 
the  ray- retardation  efifect  as  the  prime  atmospheric  source  of  error  in  range 
data.  That  is,  improved  techniques  for  estimating  the  retardation  effect  will 
serve  not  only  to  increase  the  accuracy  of  slant-range  measurement  but  will 
also  result  in  the  acquisition  of  improved  bearing-angle  data  as  well. 

2.  ATMOSPHERIC  CONTRIBUTIONS  TO  RANGE  PARAMETERS 


From  the  mathematical  statement  of  Fermat's  principle  of  least 
time,  the  time  taken  by  a  radio  ray  in  traveling  from  a  target  at  a  point,  P, 
of  its  trajectory,  T,  along  the  refracted  path,  R,  to  an  obser/atlon  point,  O 
^Figure  2),  is  given  oy  (Reference  6): 


where  dcr  is  an  element  of  the  curved  path,  and  c  is  the  speed  of  prepagation 
in  the  medium.  Introducing  the  refract;. e  Index, 

n  *  -c - V 

Equation  3  is  wri 


7, 
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Re'.rac’UvUy.  tfi  tiiv  excess  of  refractive  index  over  unity  m  parts  pt'r 

million,  is  by 

isi  -  (f,  -1/10'“.  (6) 

TbsTv^-fcrs. 

£.^At  “  t  lO'-*  \  '^'1 

"c  ‘'o 


Equ3ticf  7  espirs^a  the  electrical  range  (c^At  as  previously 
defiled)  as  sum  ot  tvo  ternia.  The  first  term  on  the  riglit-hand  side  of 
Equation  i  is  simply  the  length  of  the  refracted  ray  path,  or  geometric  range, 
Rj  the  second  is  the  contributioo  to  ascribable  to  an  atmospheric 
propagation  speed  not  equal  to  that  in  vacuo.  The  last  term  describes  the  re¬ 
tardation  of  the  ranging  signal  in  its  passage  through  the  material  medium.  By 

the  definiUon  contained  in  Equation  1,  the  slant-range  error  is  determined  by 

P 

A5  ~  lO 


Bean  and  Thayer  (Reference  3)  define  the  difference  between  the 
Imt  two  terms  on  the  right-harxl  side  of  Equation  8  as  the  "geometric  range 
erfor  "  A  ,'ig,  ana  show  that  if  does  not  represent  a  significant  portion  of  the 
total  range  error  except  at  very  small  target  elevation  angles  (between  zero  and 
iibout  3  degrees).  Evom  then  the  etiect  is  so  small  tARg  is  about  6  percent  of 
AS  when  C  ~  1  degree)  that  the  simplest  of  ray  curvature  correction  tech¬ 
niques—  for  example,  curvature  estimates  based  simply  on  surface  refractivity 


values  {Rafsrence  3)  suffice  to  reduce  it  to  negligible  proportions.  Bv  virtue  of 
these  two  facts,  the  difference  (R-S)  may  be  reasonably  considered  to  bo  zero, 
ami  the  integration  along  the  curved  path  replaced  by  integration  over  the  direct 
pati..  Under  theoe  conditions,  Equation  8  reduces  to  the  following  expression 
for  th<;  atmospheric  c(.ntributloc  of  slant-range  erior; 


Ab  10 


(9) 
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The  refora,  ron\  Equation  1, 


5  =  5^  “  lO  \NCs)<=^.'b.  (10) 

The  above  expression  indicates  that  to  determine  the  slant  range  between  target 
and  observation  point,  it  is  necessary  to  measure  the  ray  transit  time,  then 
multiply  this  by  the  vacuum  speed  of  iight,  and  diminish  the  result  by  a  quantitj’ 
that  represents  the  integrated  effect  of  the  slowing  of  the  ray  by  the  intervening, 
material  atmosphere. 

(2)  Elevation  angle  —  interferometric  measurement 

Target  elevation  in  radar  interferometry  is  determined  as  showm 
in  Figui-e  3.  From  the  law  of  cosines, 


cos 


s.. 


2  1 
i 


(11) 


where  Sf 


ft  0^ 


2 


(12) 


where  »  5^  -  $2  • 


In  the  above  expressions,  is  the  measurement  performed  by  the  mother 
station  at  Point  0  through  internal  phase  reference,  end  the  range  differences 
Sj)^  =  (S^  -  Sj)  and  Sjj^  =  (S^-  S„)  are  measured  by  the  comparison  of  the 
phases  of  the  signals  received  at  the  mother  and  at  the  slave  stations  located 
at  Points  1  and  2. 

A  further  relationship  to  be  noted  from  the  geometry  of  Figure  3 
is  that  the  projections  of  S  onto  the  baselines  B  and  B  are 

U  X  ^ 

<f)^  j  (13) 
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I 


and 


>6^  "  ■’^2  ) 


(14) 


respectively.  U,  without  loss  of  generality.  and  B2  are  considered  orthogonal 
and  as  defining  the  x-y  axes,  respectively,  of  tne  Cartesian  coordinate  system, 
then  the  cosine  functions  define  simply  the  x  and  y  direction  cosines  of  the  slant 


range,  S^; that  is. 


d  - =  i 
<*1  So 

-4 


aj  (15) 

h) 


Recognizing  that  the  sum  of  the  squares  of  the  three-dimensional  direction 
cosines  is  equal  to  unity: 


Y  ~  I  y 


then,  from  Equation  15, 


-  C  cci,^  (p  ^  ^  ^ 


tl6) 


where  n  now  defines  the  direction  cosine  to  the  z-axis  (the  target  zenith  ingle), 


i  ^  *  5m  ^  V-  'CPi 


(17) 


Therefore,  /i 

Equation  18  thus  defines  Urget  elevation  angle  in  terms  of  the  measured  inter- 
urometer  angles  /  ^  and  ^^2*  determine  and  /  ^ 

(Equations  11  atxi  12)  are  as  given  by  Equation  10.  Une- integrated  refractivity 
emerges  again  as  the  basic  atmospheric  parameter  to  be  acquired  to  effect  cor¬ 
rections  to  elevation-angle  measurements.  Also,  consider  that  the  range  differ¬ 
ences  (Sq  -  S^)  and  (S^  -  S2)  may  be  written,  in  line  with  Equations  9  and  10.  as 


9 


or 


(19) 


(Se  “AV  -AS,) 

Sp,  -  S  -  (AS-  -AS.  ) 

0  ^  e  0^  0  I 

and  \  -  AS^)  -  (  "  AS^) 

~  <rAS^  -  A  $2^  )  (20) 

That  is,  the  desired  range  differences  are  determined  by  the  differences  in  elec¬ 
trical  ranges,  diminished  b;'  the  differences  in  the  atraospheric  slant  range  con¬ 
tributions.  The  need  for  line- integrated  refractivity  measurement  at  each  of  the 
interferometer  stations  is  immediately  apparent. 

(3)  Elevation  angle  —  trilate ration 

A  three-station  trilateration  mode  of  operation  to  determine  target 
elevation  is  simply  a  special  case  of  the  interferometric  operational  mode.  That 
is,  referring  again  to  Figure  3,  the  parameter  requirements  —  the  slant  range 
S  and  range  differences  (S  -  S  )  and  (S  -  S  )  —  are  the  same,  and  the  data 

U  V  X  U 

is  interpreted  according  to  the  dictates  of  Equation  18.  This  will  vary  only  in 
that  range  difference  is  no  longer  measured  directly  through  mother- to- slave 
,  lase  comparisons,  but  is  obtained  by  separately  measuring  slant  range  from 
each  of  the  three  stations;  again,  therefore,  there  is  he  requirement  for  line- 
integrated  refractivity  to  be  measured  at  all  three  ranging  stations. 

Similarly,  tarv^t  elevation  can  be  precisely  measured  from  just 
t'TO'  ,anging  stations  (Figure  4).  Such  an  approach  requires  measuring  the  target 
azimuth  angle,  and  from  each  of  the  two  stations.  But  because  the 

atmosphere  exhibits  more  homogeneity  in  the  horizontal  than  in  the  vertical,  the 
precision  of  measurement  of  azimuth  as  compared  to  elevation  angles  is  corre¬ 
spondingly  greater.  The  work  of  Martin  and  Wright  (Reference  7)  shows  that  (1) 

-4 

the  ratio  of  horizontal  to  vertical  ray  bending  'S  on  the  order  of  10  ,  and  (2)  even 

for  very  low  elevation  angles  (  c  <  1)  a  ray  cannot  be  bent  significantly  out  of  its 
initial  vertical  plane  even  though  the  ray  azimuth  dlffe  a  from  an>'  horizontal 
refractivity  gradient.  Accordingly,  the  law  of  cosines  applied  to  the  geometry 
of  Figure  4  permits  the  statement  that 
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C£>S  € 


(21) 


With  atmospheric  contributions  to  errors  in  measurement  of  and  yi  ^ 
considered  negligible,  the  precision  to  which  the  slant  range,  S,  is  specified 
determines  again  the  precision  of  €.  . 

B.  Target  velocity 

Target  location  and  velocity  are  the  primary  concerns  in  tracking 
operations.  As  shown  in  Figure  5,  target  velocity  is  obtained  by  resolving  the 
target  velocity  vector,  V/  ,  into  its  radial,  and  angular,  ,  components. 
The  problem,  then,  is  to  obtain  the  measures  of  and  . 

(1)  Radial  velocity 

The  radial  component  of  target  velocity  is  given  simply  by  the  time 
rate-of-change  of  slant  range;  i.e. , 


dS 
d  t 


(22) 


Acting  on  Equation  10  in  accordance  with  Equation  22, 


Or,  changing  symbolism. 


(23) 


(24) 


where  the  dot  represents  the  time  derivative,  and  the  overbar  the  integrated 
refractivity  value.  Equation  24  states  that  radial  velocity  is  given  by  the  "elec¬ 
trical  velocity"  diminished  by  the  time  rate-of-change  of  line- Integrated 
retractivity. 

(2)  Angular  velocity 

The  angular  component  of  the  target  velocity  vector  is  simply  the 
product  of  slant  range  by  angular  velocity;  that  is. 


12 


(25) 


0 


gur«  5:  DEFINITION  OF  TARGET  VELOCITY  COMPONENTS 
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Differentiatlog  Equation  18  by  time,  Equation  25  (for  the  interferometric  or 
thre'i- station  trilateration  measurement  of  €  )  can  be  expressed  as 

S 

where,  from  Equation  9  and  Equations  11  and  19. 

1  t  '  ^]  ^  (27, 


Similarly,  from  FqL.atioos  12  and  20 


-  lo~^N  ll 


Equations  27  and  28  thus  call  for  the  measurements  of  the  time  rates  of  change 
of  the  electrical  range  measured  from  the  mother  stations,  of  the  slave-to- mother 
station  electrical- range  differences,  and  of  the  line- integrated  refractivity  as 
measured  at  each  of  the  three  sUit* ons,  in  order  to  evaluate  C  . 

For  the  two-atatxcn  mode  of  trilateration,  differentiation  of 


Equation  21  yields 


(z.  CSC  fc 


Here,  the  evaluation  of 


U--*’  sj  V  Sq  at-  f  , 

-  -t,  (/.„  ay.^)] 

)f  calls  for  measurements  of  t£e  Urn 


time  rates  of  change 


of  electrical  range,  of  the  azimuth  angles,  and  again  of  the  line- Integrated 
refractivity. 

(3)  Doppler  corrections 

The  alternate  approach  to  the  measurement  of  is  based  on  the 
Doppler  effect,  the  general  concept  of  which  is  that  the  reflection  of  a  radio  wave 
from  a  moving  object  results  in  a  sljnal  being  returned  to  a  fixed  tra'nsmittlng- 
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receiving  syntem  whose  frequency  differs  from  the  transmitted  frequency. 
Analytically,  th*8  phenomenon  is  written  as 


(30) 


where  f .  is  the  difference  between  the  received  and  the  emitted  frequency,  f, 
d 

c^  the  spt^ed  of  light  at  the  target  position,  P,  and  is  the  component  uf 
taigei  velocity  in  the  direction  of  the  refracted  ray.  Figure  6  illustrates  the 


geometry  involved.  According  to  Figure  6,  the  following  velocity  components 
may  be  ue fined; 


Vg  =  1^1  Q 
V  ™  ^0  r  ^  ) 


d) 

b;  (31) 

<-) 


where  is  as  defined  previously,  is  the  component  of  V  n  the  ray 
direction,  and  the  component  in  the  direction  of  the  apparent  path  (P*  is 
the  apparent  target  position). 

Doppler  frequency- shift  mea  mrements  provide  measures  of  V^. 

What  is  desired  ib  V  .  Therefore,  the  Doppler  measurement  error  may  be 
s 

defined  by 


AV^  ^V^-V  (32; 

J  r 

Substituting  appropriately  from  Equation  30 


Av^  ^  Ivi  c^^e-wi  O.S  ce  ^  %)  (33) 

Expanding  the  second  cosine  term, 

AV^  (VI  Q  —  \'7l  (c»j  0  (34a) 

or 

‘-o30M~tD5^)  lvla<o0i'<'^  (34b) 
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Figure  6; 


DOPPLER  VELOCITY  COMPONENT  IN  RELATION  TO  THE  RADIAL  AND 
ANGULAR  COMPONENTS 
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(35) 


^  is  a  very  small  angle.  Therefore,  to  good  ai^roxJmations, 

and  Equation  34b  transforms  to 

A\/  ^^!7la:.'ve 


From  Figure  6,  it  is  noted  that 

vy 

sua  *  cr>a  ri- 

Ivi 

Therefore,  Equation  36  may  be  written 

AV  -  P  V 

r  ^  € 


(37) 


(38) 


Equation  38  illustiates  that  the  Dof^ler  velocity  error  in  the  radial  direction 
attributable  to  refraction  is  composed  only  of  the  angular  component  ol  target 
.^ity,  and  is  thareiore  a  maximum  when  the  velocity  vector  is  perpendicular 
to  the  direction  of  the  direct  path,  and  a  minimum  when  tne  target  is  moving 
along  the  direct  llne-of-siglit  path.  Further,  notice  that  (1)  unlike  the  range 
measurement  error,  the  Doji^ler  error  is  not  cumulative  '(it  is  only  tliC  AX'alized 


error  at  a  specific  target  position  in  the  medium),  and  (^/  the  magnitude  ot  th(' 
pointing  error  angle.  <£  ,  la  not  signifinat. 

Substituting  Irom  Equations  24,  30,  and  38  appropriately  into  32 


^5  -  10 

^  €  e 


_L  lA. 

Z  f 


(30) 


Introducing,  now,  Equation  4 


gv  *s 

t  e  ^  n  f 


(40) 
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and  the  above,  tiirough  application  of  Equation  S  and  to  first  order  in  (n-1) 
is  written 


~S  -  10  M 

e 


fo 

»o 


or 


2 


-.c, 

10  N 


Z 


(41) 


(42) 


The  term  equation  is  simply  the  target  velocity  that  would 

be  measured  by  a  Doppler  system  in  vacuo.  This  quantity  is  thus  equivalent  to 
the  previously  defined.  Therefore,  Equation  42  simplifies  to 

AV^  -  -  10  M  pSg  ~  >0~^  N  (43) 


Hence,  in  terms  ('f  the  Donpier  velocity,  the  radial  "e'ocity  is 

Vg  "  ^  P  ~  ^ 

Tne  Doppler  measurement  cor.-ection  ic  that  required  for  the  raoial  velocity 
measurements  plus  the  term  that  accounts  for  refractivity  conditions  at 
the  target  position.  For  ranges  long  enough  and  elevation  angles  fc,raat  enough, 
the  latter  term  rap  uy  ’uses  its  significance  booaube  o*  the  rapid  decrease  of  N 
with  height.  In  t  ict,  after  the  target  leaves  ti.e  atmosphere,  N  -  0,  and  the 
Doppler  coriection  is  exactly  that  required  tor  radial  correction. 

C ,  summary' 

Tw^  ..tmospbeidc  effects  aci  lo  cause  errors  in  radar  ranging  and  track¬ 
ing  date:  H)  the  material  atmosphere  retards  ranging  signals  so  that  ray  propa¬ 
gation  apeeciri  are  less  in  Use  atmosphere  than  in  vacuo,  and  ^2)  gradients  (espe¬ 
cially  in  the  vertical)  in  the  refractivity  st.ructure  of  the  atmosphere  cause  a  ray- 
bending  or  refraction.  The  first  effect  results  in  errors  in  slant- range  measure- 
meiua,  the  secorJ  effect  results  in  errors  in  elevation-angle  deternunatiuns. 

In  the  inter'-rometric  or  triiate ration  mode  of  ranging  oper.ations,  the 
dliect  measure  of  target  elevation  is  elimiaated;  i.  e.  ,  target  position  is  com- 
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pletei/  specified  from  knowledge  of  receiver-to-target  slant  ranges  relative  to 
mi  Jtiple. ringing  stations  in  a  prescribed  geometric  configuration.  Moreover, 
because  of  the  purely  mechanicsil  orohlems  of  maintaining  dish  orientation  in 
radars  that  measure  elevation  angle  directly,  it  is  likely  that  accuracies  in 
elevation- angle  specifications  commensurate  with  the  potential  accuracies  of 
slant- range  determinations  are  attainable  only  with  the  multiple-station  approach. 

Improvement  in  the  accuracy  of  ranging  data  for  both  position  and 
Vslocity  demands  improvement  in  the  measurement  of  slant  range.  The  domi¬ 
nant  need  is  for  techniques  only  for  estimating  the  ray- retardation  effect,  which 
is  a  function  simp'v  of  the  integral  of  refractivity  over  the  ranging  path.  The 
difficult  problem  of  performing  the  measare  in  real  time  of  the  refractive  index 
gradients  required  for  the  corrections  of  elevation  angles  directly  measured  is 
thus  eliminated. 

3.  N-MEASUREMENT  TECHNIQUE  SURVEY 
A.  Background 

The  description  of  n  given  by  Equation  4  is  purely  phenomonological. 

It  is  the  consideration  of  the  actual  effects  of  an  impressed  electric  field  upon 
the  dielectric  constants  of  both  polar  and  nonpolar  molecules  that  establishes 
the  physical  basis  of  a  nonzero  refractivity,  and  that  relates  refractivity  to 
atmospheric  composition.  By  such  consideration  has  come  a  refractivity  equa¬ 
tion  of  the  form  (Reference  8) 

iS  -  A  (Jy  ^  4 

where 

A,  B  ,  and  B  in  the  above  expressnon  semiemnirically  determined  con 
stants;  T  is  the  temperature,  p^  is  the  partial  pressure  of  dry  air,  and  that  of 
water  vapor.  The  total  atmoepheric  pressure  is  thus 

d 
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The  pecuiaritiea  of  Equation  45  are  that  at  optical  wavelengths,  and 
B  are  effectively  zero,  and  N  is  therefore  largely  independent  of  the  atmos- 
pherlc  water-vapor  content;  the  constant  A,  however,  is  a  function  of  wavelength. 
Conversely,  at  radio  wavelengths,  refractivity  is  not  wavelength-dependent  (  A 
independent  of  wavelength)  but  is  strongly  affected  by  water  vapor  (Reference  9). 
Therefore,  it  is  convenient  and  customary  to  present  two  equations  lor  retrac- 
tivitj':  one  for  optical  wavelengths  and  the  other  foi  radio  wavelengths;  that  is. 


N 

ou 


=  (k,  - 


Id 

T 


(47) 


where  denotes  the  refractivity  at  the  optical  wavelengths  («20  microns), 
and 


hi 


K 


(48) 


for  radio  wavelengtiis.  Regarding  the  radio  wavelength  refractivity  equation: 
of  special  significance  (to  be  emphasized  later)  is  the  difference  in  form  of  the 
first  two  terms  on  the  right-hand  side  as  compared  to  the  third  term.  This 
difference  arises  because  the  first  two  terms  express  the  effects  on  refractivity 
of  the  distortions  of  electronic  charges  of  dry-gas  and  water-vapor  molecules 
under  the  influence  of  an  applied  electromagnetic  field,  whereas  the  third  term 
describes  the  effect  of  molecular  dipcle  orientation.  Of  all  the  constituents  of 
the  atmosphere,  only  water  vapor  exhibits  a  strong  dipole  moment;  therefore 
only  water  vapor  contributes  a  dipole-orientation  term  to  the  refractivity  equation. 
Also  note  thnt  the  first  terir.  in  tne  optical  and  radio  refractivity  equations 
are  idei*,.ical. 

B.  Point- measurement  technique 

Using  Equations  4  and  48,  radio  refractivity  may  be  measured  in  two 
ways.  That  is,  N  may  be  obia.ned  directly  through  a  measure  ■'f  the  ratio  c^/c, 
or  else  by  calculation  (using  Equation  48)  from  measures  of  the  atmospheric 
parameters  p^,  T,  and  e. 
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(^)  Dtroct  meaaurement  of  N 

A  number  of  deyices  are  capable  of  directly  measuring  radio  re- 
fractivity.  Although  nuiny  modifications  have  been  attempted,  moat  such  retrac- 
tometers  fall  into  three  basic  categories  that  reflect  the  essential  principles  of 
the  original  units.  The  three  t^’pe^  —  resonant  cavities,  coaxial  cavities,  and 
condensers  —  ail  involve  comparisons  between  a  sample  and  a  reference  quantity 
and  constitute  a  point  measurement.  The  most  noted  of  the  resonant  frequency- 
type  refractometers  are  those  of  Crain  (Reference  10)  and  Bimbaum  (Reference  11). 

The  former  measures  the  variation  in  frequency  difference  between  two  stabilized 
oscillators;  the  difference  is  proportional  to  the  change  in  refractivity.  The 
Birnbaum  instrument  uses  passive  cavities  instead  of  the  active  ones  employed 
by  Crain,  and  it  measures  the  time  delay  in  response  to  a  peak  signal  passing 
through  the  reference  and  sensing  cavities.  This  time  lag  represents  the  fre¬ 
quency  difference  and  is  a  measure  of  the  refractivity.  Cavity-resonating  metru 
ments  typically  operate  near  10,000  MHz,  and  both  the  Birnbaum  and  Craia 
instruments  have  a  number  of  selective  scales  tnat  pemit  nie*s8urementB  over 
a  total  raiige  of  some  400  N  units,  For  the  most  part,  the  many  mtxlificationa 
of  the  Birnbaum  unit  result  from  various  techniques  used  vo  determine  the  time 
delay  of  maximum  response.  System  noise  is  approximately  0, 1  N  unit. 

The  coaxial,  cavity  refractometer,  fxjst  represented  by  the  instru¬ 
ment  developed  by  Deani  (hoierence  Ili),u8e8  a  half-weave  resonant  aampluig 
cavii;,  cpt'i  .iting  •’t  a  frequenev  of  400  MHz.  The  basic  instixi-vu'iit  iu  ..  m'Xiified 
Crain  unit  that  has  been  designt-d  especially  for  lightweight  balloon  or  drop- sonde 
applicatu)n.s .  since  the  sampling  cavity  might  be  Icw-^ted  as  tar  as  20  miles  from 
the  reference  receiver, 

A  somewhat  different  approach  was  ir.hen  by  Hay  (Reference  13)  in 
applying  condensers  m  a  resonant  circuit,  wherein  the  capacitance  varies  accord¬ 
ing  to  the  changes  In  dielectric  constant  of  iht  gas  sanipied  between  the  plates. 

The  original  instrument,  operating  ai  10  MHz,  ised  six  plales  spaced  airout  0.25 
i.-.cries  'nart.  on  each  aide  oi  -...e  i r.struir.ert ,  The  primary  advantage  of  this  unit 
was  that  it  could  be  operated  remotely  from  its  reference  unit,  and  because  it 


was  lighter  than  cavity  resonators  it  was  ideally  suited  to  balloon  operations. 
Handicaps  occurred,  however,  from  problems  of  temperature  calibration  and 
frequency  response. 

As  early  as  1951,  Tolbert  and  Straiton  (Reference  14)  introduced 
a  design  significantly  different  from  existing  instruments;  it  was  an  early  effort 
toward  developing  integrated  refractivity  devices.  Microwave  transmissions 
were  sent  over  a  1- meter-long  measuring  and  reference  path,  and  a  comparison 
of  the  phase  shift  of  the  received  signal  was  interpreted  in  terras  proportional  to 
the  mean  refractive  index  over  the  path.  It  was  determined  that  results  obtained 
with  this  system  could  be  made  extremely  sensitive  to  small  changes  in  refrac¬ 
tivity,  although  disadvantages  appeared  in  attempts  to  thermally  balance  the 
system  and  in  instability  of  the  amplifier  components 

(2)  Calculation  of  N  from  meteorological  data 

The  constants  K  ,  K  ,  and  K  in  Equation  48  are  of  the  values 
1  *5 

(Reference  15); 
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The  K  values  contained  in  Equation  49  were  evaluated  by  Bean  from  a  survey 
of  determinations  of  tlie  dielectric  constants  of  both  di7  air  a,id  water  vapor 
conducted  by  various  researchers.  The  final  results  reriesent  the  weighted 
means  of  the  various  determinations,  the  weights  being  t.aken  inversely  propor¬ 
tional  to  the  square  of  the  probable  errors.  In  rounding  the  K  values  to  three 
figures  where  significant,  and  assuming  the  total  air  pressure  to  equal  the  sum 
of  the  dry-air  and  water-vapor  pressures,  the  foliowing  tamr  refractivity 
equation  obtains: 

N  -  7  7  4,  ^  -tr  7Z  —  ,'yJ.75A  lO^)  •  (5h) 
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The  K  values  specified  io  Equation  50  are  thoje  suggested  by  Smith  and  Weintraub 
(Reference  16),  and  are  stated  as  being  good  to  0.  5  percent  in  N  for  the  normal 
ranges  of  temperatures,  pressures,  and  humidities  encountered  in  the  atmosphere. 
In  recognition  of  their  work,  Equaticn  50  is  usually  called  the  Smith- Weintraub 
equation.  The  basic  uncertainty  of  the  refractivity  equation  is  compounded  by 
errors  of  measurement  of  p^  T,  and  e.  With  errors  common  in  radiosonde 
measurements,  the  residual  error  in  N  would,  in  fact,  be  closer  to  something 
like  1.  5  percent  ^uL■  on  the  order  of  a  few  N  units).  The  advantage  of  the  direct 
over  the  indirect  measure  of  N  is  obvious. 

(3)  Integrated  refractivity  value 

By  far  the  greatest  amount  of  effort  toward  supplying  the  integrated 
refractivity  demanded  in  ranging  operations  has  been  expended  in  searching  for 
techniques  to  interpret  point- measured  refractivity  in  terms  of  the  integrated 
value.  In  a  perfectly  homogeneous,  quiescent  atmosphere,  any  of  the  point- 
measurement  techniques  listed  above  could  be  used  to  determine  such  integrated 
data,  the  value  simply  being  the  measurement  obtained  at  one  point  c  f  a  given 
propagation  path  multiplied  by  path  length.  Tacit  here,  is  the  assumption  that 
the  single-point  measurement  provides  n.n  accurate  estimate  of  the  average 
refractivity  value  all  along  the  path.  But  usually  the  atmosphere  is  not  quiescent, 
and  refractivity  is  a  function  of  position.  The  requirement  for  the  true  integral 
value  of  N  permits  the  interpretation  that,  even  though  errorless  determination 
of  N  at  each  of,  say,  k  sampling  points  was  made,  the  average  N  computed 
from  such  data  would  still  be  subject  to  a  finite  sampling  error  to  tiie  extent  that 

5 

—  ^  ^  N(s)d^  (51) 

The  above  InequiJUy  defines  the  samp,  .g  error,  and  is  the  basic  reason  that 
even  though  statistical  accunicies  of  from  3  to  3  percent  in  the  estimation  of 
atmospheric  eflects  on  range  data  are  attainable  from  radiosonde  or  refractom- 
eter  profiles  (Reference  17).  in  real  time  a  residual  error  more  lik'"  three  times 
that  amount  Is  the  more  realistic  figure  (Reference  18). 
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C.  Normal  dlaperaion  technique 

Takiog  advantage  of  the  normal  and  completely  predictable  monotonlc 
decrease  of  the  in  lex  of  refraction  with  increasing  wavelength  (1.  e. ,  dispersion) 
in  the  optical  portion  of  the  spectrum  (Equation  47)),  Owens  and  Bender  (Refer¬ 
ence  19)  have  developed  a  technique  for  the  direct  measure  of  the  optical  line- 
integrated  refractivity.  Assume  that  time-delay  (interferometric)  determinations 
are  made  at  two  wavelengths, \  ^  and  X2.  In  tbe  optical  portion  of  the  spectrum. 
TaJce  Xj  ns  snorter  wavelength.  Then  the  difference  between  tho  range 
contributions  at  the  two  wavelengths  wU!  be,  according  to  Equation  8, 


iO~^  A  ^  ^  (i»)  ds 


where 


(52) 


(53) 


The  integrated  value  of  refractivity  at  wavelength  A„  is  thus  directly  obtained 

from  range  measurements  at  two  wavelen^ltha  in  the  portion  of  the  electromagnetic 

spectrum  exhibiting  normal  dispersion. 

The  instrumentation  consists  of  sources  of  modulated  (in  the  UHF  or 

microwave  range)  red  and  blue  light  (the  6J28-ang8trom  helium-neon  laser  line 

and  the  3660-angstrom  mercury  line,  respectively).  Both  colors  are  propagated 

over  the  ranging  path  and  th?n  reflected  hack  to  a  receiver,  where  amplitude  or 

polarization  modulation  techniques  are  used  to  pro/lde  the  measure  of  SX.  ”  SA  . 

_  1  2 

With  such  a  system,  accuracies  in  the  direct  measurement  of  N  of  better  than 
1  percent  are  deemed  feasible.  Such  accuracy  in  N  measrrement  would  result 
in  about  an  order- of- magnitude  improvement  in  the  accuracy  of  slant-range 
measurement  over  that  now  achievable  through  electrical -range  corrections 
based  on  N  estimates  from  point-measurement  data.  However,  the  technique 
cannot  be  applied  directly  to  the  case  of  radar  ranging,  because,  in  effect,  only 
the  first  term  of  the  integrated  radio  refractivity  equation  is  sc/  determined. 

The  technique  provides  no  means  whereby  the  second  and  thiid  terms  nsay  be 
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measured.  Strictly  speaking,  the  technique  is  applicable  only  to  optical  ranging 
devices  such  as  the  geodimeter  (Reference  20).  Ano  J'ur  disadvantage  is  that 
cooperative  targets  on  viiich  the  ipequired  reflector  could  be  mounted  are  manda¬ 
tory  for  operation  of  the  system  in  conjunction  with  ranging  operations. 

D.  Anomalous  dispersion  technique 

In  the  radio  region  of  the  spectrum,  refractive  index  is  largely  indepen¬ 
dent  of  wavelength.  Therefore,  the  normal  dispersion  technique  of  measuring 
line- integrated  refractivity  is  not  applicable.  In  the  neighborhood  of  atmospheric 
constituent  absorption  bands,  however,  an  anomalous  dispersion  is  associated 
with  the  moU  'vlar  resonances  of  the  particular  absorbing  constituents.  Specifi¬ 
cally,  on  the  long  wavelength  sides  of  such  bands  the  refractivities  fall  to  abnor¬ 
mally  low  values,  whereas  on  the  short  wavelength  sides  such  values  arc  abnor¬ 
mally  high.  Moreover,  the  magnitude  of  the  dispersion  is  proportional  to  the 
"eneral  level  of  refractivity  about  the  band.  This  magnitude,  in  turn,  can  be 
determined  from  the  difference  in  ray-propagation  transit  times  associated  with 
the  dispersion  across  the  band.  Specifically,  ifAt^  and  At 2  denote  ray 
transit  times  —  one  above  and  the  other  below  a  given  absorption  band  —  then 

At  =■  N(s)as  (54) 

where,  as  an  example,  if  the  transit  time  difference  is  determined  across  a  water- 
vapor  absorption  band,  the  integral  represents  the  integrated  water-vapor  refrac¬ 
tivity.  Because  total  refractivity  is  the  sura  of  a  dry-air  terra  and  a  water-vapor 
term;  i,  e. , 
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(55) 


then  a  similar  raeasureraert  across  an  absorption  band  of  one  of  the  dry-air  con¬ 
stituents  will  suffice  to  determine  tl  line- Integrated  refractivity.  This  is  the 
techaic[ue  of  N  measurement  currently  being  Investigated  by  Sullivan  and  Richard- 
eon  of  the  Mitre  Corpoiation  (Reference  21). 
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Aq  instrument  capable  of  measuring  the  differential  transit  time  due  to 
water  vapor  has  been  developed  and  field  tested.  It  consists  of  a  sou;  ;e  com¬ 
posed  of  two  phase-coherent  w-aves  at  15. 6  and  32. 2  GHz  transmitted  to  tree 
space  through  a  Cassegrainian  parabolic  antenna.  The  transmitted  energy  is 
detected  by  a  system  using  a  similar  antenna  in  the  receiver.  A  local  oscillator 
Is  incorporated  to  provide  two  IF  carrier  frequencies;  one  at  250  kHz  to  carry 
the  phase  information  of  the  lower  transmitted  frequency,  and  one  at  500  kHz  to 
carry  the  information  on  the  hi^er  frequency.  The  250-kHz  signal  is  doubled, 
and  the  phases  of  the  two  are  then  compared  by  a  null-detecting  digital  phase 
shifter.  The  amount  of  phase  shift  to  achieve  the  nuii  is  the  measure  of  the 
differential  transit  time. 

Theoretically,  accuracies  in  N  measurement  comparable  to  those 
achievable  with  the  normal  dispersion  8\  stem  can  be  anticipated  with  the  Mitre 
system.  However,  field  test  with  the  water-vapor  system  have  shown  that  the 
data  obtained  by  the  system  consist  of  an  a'^erage  level  and  a  widely  variable 
fluctuating  component  about  the  average  level  (Reference  22).  The  average 
level  is  determined  by  the  integrated  humidity  in  the  path,  whereas  the  majority 
of  the  noise  component  is  due  to  random  phase  scintillation  effects  associated 
with  the  propagation  of  electromagnetic  energy  through  a  turbulent  atmosphere. 
Consequently,  prohibitively  long  averaging  times  (with  the  attendant  loss  of 
real-time  capabilities)  are  required  to  recover  mc.mingful  data. 

Further,  the  transit  time  differentials  that  are  ascribable  to  the  anoma¬ 
lous  disjjersion,  and  which  must  be  measured,  are  on  the  order  of  picoseconds 
-12 

(10  seconds).  A  capability  for  such  measurement  severely  taxes  the  state  of 
the  art  of  microwave  electronics.  Although  th«  Mitre  system  represents  a  sig¬ 
nificant  achievement,  and  although  the  expenmenval  system  has  proven  Itself, 
there  might  be  problems  In  maintaining  the  stability  and  reliability  of  such  a 
system  in  routine  ranging  operations. 

Thi  linal  drawback  is  that,  again,  cooperative  targets  would  be  required 
when  applying  the  system  to  ranging  operations. 
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E.  Summary 

The  need  to  measure  line-integrateu  refractivitj  in  ranging  operations 
has  resulted  in  extensive  research  into  techniques  for  estimating  N  from  point- 
measurement  refractivity  data.  The  limit  of  such  techniques  has  just  about  been 
reached;  still,  uncertainties  of  from  5  to  10  times  the  magnitude  desired  remain 
in  the  real-time  specification  of  N.  The  problem  is  the  inadequacy  of  the  repre¬ 
sentation  of  an  integrated  function  as  a  finite  sum.  Moreover,  the  sheer  problem 
of  logistics  in  providing  an  adequate  number  of  point  measurements  in  real  time 
to  attain  the  desired  degree  of  precision  in  the  estimated  integrated  value  prac¬ 
tically  precludes  the  use  of  point- measure  mem  techniques  in  ranging  operations. 

Two  a(^roaches  to  the  direct  measure  of  the  true  integrated  function 
have  been  suggested.  Both  are  based  essentially  on  the  concept  that  dispersion 
(wavelength-dependent  refractivity)  results  in  finite  ray  transit  time  differences 
at  different  radiation  wavelengths.  Measures  of  the  magnitudes  of  such  diffeionces 
therefore  provide  the  measure  of  integrated  refractivity.  The  normal  dispersion 
at  optical  wavelengths  io  uaed  to  provide  information  regarding  the  optical  refrac¬ 
tivity;  the  anomalous  dispersiou  around  atmospheric  constituent  absorption  HinHs 
can  be  used  similarly  to  determine  the  radio  refractiviu .  Accuracies  in  N 
measurement  theoretically  attainable  with  each  system  would  more  than  satisfy 
the  requirements  of  ranging  operations.  However,  the  normal  dispersion  tech¬ 
nique  suffers  because  (1)  it  provides  N  data  applicable  only  to  optical  wavelengths, 
and  (2)  cooperative  targets  would  be  needed  if  the  approach  wore  applied  to  ranging 
operations.  The  anomalous  dispersion  technique  suffers  because  (1)  a  capability 
for  measuring  short-term  N  fluctuations  is  lacking,  (2)  the  system  is  borderline 
technically,  and  (3)  cooperative  targets  again  would  be  neeifed  for  its  application 
to  ranging  operations. 

4.  U.U)1ATI0N  ABSUH1>TK)N  AND  THE  I^yiAMETEKS  OF  t  INK- INTEGHATEU 
HEFIIACTIVITY 

.As  a  basis  tor  an  a.ternate  approach  to  tiie  measurement  of  line-integrated 
refractu  ity ,  consider  f’le  folloM-ing.  Dasically,  refraction  effects  constitute  a 
phenomonological  record  of  the  influence  of  the  atmosphere  on  the  phase  of 


propagated  electromagnetic  signals.  But  the  atmosphere  acts  not  only  on  wave 
phase,  but  also,  at  certain  precisely  defined  frequencies,  on  wave  ampJ-Mde. 

This  latter  action  gives  rise  to  the  familiar  phenomenon  of  energy  absorption. 

The  net  result  is  that  not  only  is  the  energy  propagation  speed  altered,  but  so 
is  the  actual  spectral  content  of  that  energy  (within  specified  frequency  intervals). 
Although  these  two  phenomena  may  have  widely  different  physical  manifestations, 
it  is  important  to  recognize  the  common  origin  of  both  —  the  interaction,  at  the 
molecular  level,  of  the  atmosphere  with  an  applied  electromagnetic  field.  It  is 
not  too  surprising  that,  in  spite  of  the  differences  mentioned,  a  phenomenological 
descript. on  of  one  effect  can  be  interpreted  in  terms  of  a  description  of  the  other. 

A.  Equation  of  radiative  transfer 

Regarding  the  absorption  process  mentioned  above:  The  totality  of 
molecular  interactions  between  radiation  and  matter  can  be  classed  as  either 
absorption  or  emission.  The  two  processes  are  diatinguished  by  the  sign  of  the 
cha.nge  of  radiant  intensity  as  a  result  of  the  interaction.  If  the  intensity  tiecreases, 
then  absorptica  occurs;  intensity  increases  are  ascribed  to  emission. 

The  fundamental  law  of  absorption  is  that  of  Bou^^-uer.  It  states  that  for 
monochromatic  radiation  the  absorption  process  is  liuear  independently  of  the 
intensity  of  the  radiation  and  in  the  amount  of  matter,  provided  that  the  phy  sical 
state  is  held  constant.  The  mathematical  statement  of  Bouguer's  law  is  as 
follows:  il  the  space  rate  of  diminution  of  radiant  intensit}-  along  the  direction 
of  propagation  is  defined,  then  the  change  m  intensity  in  traversing  an  inlinitesi- 
mal  element  of  this  path  is 

(•Wjor^tiCn  )  —  (X  ^  ots  (5b) 

^  c 

where  I  .  is  an  intrinsic  source  Intensity,  ds  is  the  infinitesimal  path  ele- 
meat,  and  is  the  previousl}  defined  space  rate  of  intensity  diminution. 

This  latter  quantity  is  usually'  called  the  "volume  abso.ption  coefficient.  " 

The  argument  that  the  abso  ptlon  process  Is  lineai  in  the  amount  ol 
matter  also  applies  with  equal  force  to  the  emission  pr.  ''ss,  .As  a  formal 
statement, 
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(57) 


where  is  defined  as  the  source  function  of  the  emitting  media. 

The  earlier  statement  that  all  Interactions  can  be  classed  as  absorption 
or  emission  can  now  be  summed  up  In  the  statement  that  an}'  change  in  intensity 
resulting  from  the  interaction  of  matter  and  raalatlon  must  be  the  cum  of  Eqi'ations 
56  and  57;  1.  e. , 

V)Sc>rpt»o>v)  T"  (58) 

Using  the  definitions  provided  by  Equations  56  and  57,  the  solution  to  Equation 
58  is  (Reference  23) 

=  l^e  ^  4  (59) 

o 

where  again  the  first  term  on  the  right-hand  side  of  Equation  59  describes  *he 
absorption  (as  governed  by  the  term's  expijnential  attenuation  factor)  of  the 
intrinsic  source  energ}-,  the  second  term  describes  the  self-emission  process. 

The  exponential  tac*^or  also  appears  in  ♦''e  second  term,  expressip^  the  fact  that 
an}’  energ}'  emitted  at  any  point  along  the  propagation  p>ath  suffers  the  same  ab- 
soiption  in  its  subsequent  travel  as  does  the  intrinsic  source  er^ergy.  While 
Equation  58  and  its  solution  (Equation  59)  sets  the  pattern  of  the  formalism  used 
in  radiation  transfer  problems,  its  physical  content  is  very  slight.  The  ph}'8ic8 
is  mainly  contained  m  the  defioitioo  of  at^  . 

B.  Specification  of  the  ab»o;:ptlon  coefficient 

So  far,  the  absorption  coefficient  has  been  given  only  a  phenomonological 
description.  The  actual  ph}’8’f3:S  of  tii«  absorptioa  process  lies  in  the  interactions 
between  the  incident  radiation  aud  the  energy  modes  of  the  individual  molecules. 
Specifically,  rediation  is  absorbed  when  it  is  of  the  proper  liequency  (as  deter¬ 
mined  by  the  Bohr  frequency  relation)  to  cause  a  molecule  to  be  raised  from  one 
molecular  state  to  .iiother  of  higher  energy. 
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In  a  general  way,  the  energy  of  a  molecular  state  may  be  writtei.  as  the 
sum  of  four  terms; 


E  -Kt  4-E 

e  V  r  t 


(60) 


where  E  ,  E  ,  E  ,  and  E,  are  the  electronic,  vibrational,  rotational,  and 
e  V  r  t 

consequently,  the  frequencies  of  the  radiation  required  to  effect  transitions  in 
the  energy  levels  associated  with  each  of  the  four  terms  vary  in  the  same  manner. 
That  is,  E^  transitions,  caused  by  the  highest  frequency  radiation,  produce  band 
spectra  that  lie  in  the  ultraviolet  and  visible  regions  of  the  electromagnetic  spec¬ 
trum.  Pure  vibration,  vibration- rotational,  and  pure- rotation  spectra  occur, 
progressively,  in  the  near  through  far  infrared.  Translational  spectra  are 
restricted  to  the  radio  regions. 

(1)  G  ral  quantum-mechanical,  molecular  absorption  coefficient 
Electric  or  magnetic  dipole  or  quadrupole  moments  are  involved  in 
the  interactions  between  matter  and  incident  electromagnetic  fields  that  result  in 
energy  absorption.  Electric  dipole  interactions  are  stronger  by  a  factor  on  the 
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ordei  of  iO  thau  magr  Jtic  dicole  interactions;  electric  dipole  interactions  are 
on  the  order  of  10  times  stronger  than  electric  quadrupole  interactions.  Elec¬ 
tric  dipole  transitions  are,  therefore,  responsible  for  the  strongest  spectral 
lines,  and  are  called  "permitted"  transitions.  Other  transitions  are  loosely 
named  "forbidden.  "  For  a  single  spectral  line,  the  absorption  coefficient  is 
related  to  the  intensity  of  the  permitted  dipole  transition,  according  to  (Peference 


translational  energies,  respectively.  It  is  further  found  that  E^ 
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where  Rjj  is  the  matrix  element  of  the  dipole  moment  connecting  two  stationary 
states,  i,  j,  of  energy  E.,  E^,  respectively;  z/ .j  is  the  frcciuency  of  the  corre¬ 
sponding  spectral  line,  "ivenb\  the  Bohr  frequency  condition  h2/>  -  E.  -  E-  ; 
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n  is  the  number  density  of  the  absorbing  molecular  specie,  T  is  the  tempera- 

fl 

ture,  h  is  the  Planck  constant,  k  is  the  Boltzmann  constant,  and  c  is  the  speed 
of  light.  The  frequency  of  the  incident  radiation  is  denoted  by  p  .  The  exponen¬ 
tial  is  termed  the  Boltzmann  factor,  and  the  factor  f(>/_,Z/  )  is  the  so-called 
"structure  function"  that  deternained  the  shape  of  the  absorption  line.  For  the 
Lorentz  pressure-broadened  line  shape,  *  the  expression  for  f([/_  ,  v  )  is 
(Reference  24); 
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where  Av  is  the  line  width;  i.e. ,  the  half- width  of  the  line  at  half  intensity. 
Near  resonance,  Equation  62,  can  be  written  as 


if  -j 

tr  !_ 


(63) 


Equation  63  is  a  valid  representation  of  the  structure  function  at  optical  and  near- 
through  middle- infrared  wavelengths.  Microwave  absorption  line  structures, 
however,  must  be  described  in  terms  of  Equation  62. 

(2)  as  a  function  of  frequent y 

Equation  61  is  the  general  (juantum- mechanical  expression  for  the 
molecular  absorption  coefficient.  As  hx',  ,  =  E,  -  E  =  -hp.,  it  follows  that 

ij  1  j 

f(l/j P  )  -  ,P  )■  To  any  given  term  i,  j  in  the  double  sum  called  tor  in 


^  Atmospheric  prcssuies  in  the  nortion  of  the  atmosphere  where  there  is  any 
significant  absorption  allow  the  Lorentzian  line  shape  to  be  assumed.  Ooppler 
line-broadening  elfects  are  significant  only  in  the  high-altitude  rarefipd  atmos¬ 
phere  tiiat  contributes  little  to  absorption,  and,  hence,  refractivity. 
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Equation  61,  there  also  a  corresponding  term  — J,  1  — whose  contribution  is 
of  the  opposite  sign,  and  which  differs  from  the  1,  j  term  only  to  the  extent  of 
the  difference  in  the  Boltzmann  factor  for  the  J-  as  opposed  to  the  i-state.  If, 
in  fact,  hy^^«kT,  a  condition  characteristic  of  microwave  wave  frequencies, 
the  effects  of  the  i,  j  and  ],  i  terms  very  nearly  cancel  each  other.  This 
difficulty  of  neaily  compensating  terms  in  the  microwave  region  of  the  absorption 
spectrum  is  avoided  with  the  following  approximation: 
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(64) 


Use  of  the  above  simplification  in  Equation  61  deflncr  now,  the  microwave 
absorption  coefficient  as 
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3y  performing  the  Indlcatcc.  cummations  over  the  contributing 
energy  states  exhibiting  appreciable  population  at  the  temperatures  encountered 
in  the  atmosphere,  the  molecular  absorption  coefficients  for  water  vapor  and  for 
oxygen  may  be  written  (Reference  24) 


7  =  t 


1 


•i-G 


7.- 


I 


at  microwave  regions,  and  (Reference  25) 


(66) 


(b) 
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(67) 


at  iofrared  regions.  In  the  above  equations,  denotes  the  resonant  frequency 

of  the  appropriate  absorption  lines.  The  units  of  the  C's  and  B's  that  appear  in 
the  equations  are  determined  by  the  units  of  radiation  transfer  measurements. 

(3)  Volume  absorption  coefficients 

Recognizing  now,  the  relationship  that  exists  between  number  den¬ 
sity  and  mass  density;  i.e. , 


=  52-. 
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where  is  the  molecular  weight  of  constituent  a,  and  N  is  Avogadros 
number,  a  s  .gnlficant  fact  emerges  from  comparison  of  Equations  66  and  67;  viz. , 
in  terms  of  the  volume  absorption  coefficient,  absorption  at  microwave  frequencies 
proceeds  in  proportion  to  the  ratio  of  absorbing  constituent  density  to  temperature, 
whereas  at  infrared  frequencies  the  proportion  is  according  to  density  alone. 

That  IS, 
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and 
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where  the  constants  are  adjusted  according  to  Equation  68. 

C .  Specification  of  line- integrated  refractivity 

Consider  now,  the  modification  of  Equation  48  according  to  the  equation 
of  state  for  dry  air. 
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and  for  water  vapor 
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where  R  and  R'  are  the  gas  constants  for  dry  air  and  for  water  vapor,  respec¬ 
tively;  denotes  the  partial  density  of  dry  air, that  of  water  vapor,  or 
absolute  humidity.  Af^lying  Equations  71  and  72  to  Equation  48,  the  expression 
for  N  becomes 

M  ^  p  t-  K;,  R 

'  ^ol  i-  3  T  (73) 

2  -2  -1  -3 

With  the  It  and  p  values  expressed  in  units  of  in  sec  "K  and  gm  m  , 
respectively,  the  K  values  of  Equation  47  are  redefined  to  a  new  set  according  to 


K,  =  KR  =  C2.2£.7£3  ^0  03  7)xd''^  y.,'' 

K^.  =  =^330.-4  ic  vv,"'  (74) 


Equation  48  can  thus  be  written  as 
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Equation  75  c^n  be  further  simplified  when  it  is  recognized  that  the  diy-air  con¬ 
stituents  exist  very  nearly  in  constant  proportion  to  one  another.  Therefore,  the 
partial  dry-air  density  can  be  expressed  in  terms  of  the  partial  density  o*  oxygen 
according  to 


(76) 


where  denotes  now  the  partial  density  of  oxygen,  and  denotes  the  oxygen 
mixing  ratio  (i,e. ,  percent  by  weight  of  oxygen  to  the  total  diy-air  mass  per  unit 
volume) .  The  numerical  value  of  is 
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The  uncertainty  indicated  in  the  above  equation  arises  Irom  the  variations  in 
atmospheric  carbon  dioxide  content  established  by  the  release  of  CO  to  the 
atmosphere  by  combustion  processes  and  by  biological  exchanges  of  CO  with 
soil  and  vegetation.  Equation(76)  at  any  rate,  permits  Equation  75  to  be  written 
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The  integrated  refractivity  is  thus 


(78) 


(79) 


The  full  impact  of  Equations  69  siaI  70  can  be  appreciated  now,  when  a 
comparison  is  made  be,, ween  these  equations  the  radiative-transfer  equation 
(59)  and  the  N  equation  (79).  Because  of  the  appearance  of  the  integrals  of  the 
ansorption  coefficient  in  thu  exponential  radiation  attenuation  factors  of  the  radi¬ 
ation  transfer  equation,  the  physics  of  radiative  transfer  defines,  in  effect, 
measurable  quantities,  ly's,  which, depending  on  the  frequency  of  rat’ietion 
considered,  are  functions  of  precisely  the  quantities  'required  for  an  evalua¬ 
tion  of  line- integrated  refractivity  as  specified  by  Equation  79.  That  is, 
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Inversions  of  Fauations  80  therefore  yield 
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which  are  the  atmospheric  input  parameters  demanded  by  Equation  79.  There¬ 
fore, ‘Equation  79  may  be  written 

where  k^'  is  defined  as  Transmission  measurements  performed  at 

infrared  frequencies  in  an  oxygen  and  in  a  water- vupor  absorption  band,  respec¬ 
tively,  therefore  determine  the  first  two  terms  of  the  line- integrated  refractivity 
equation,  and  similar  measurements  performed  in  a  water-vapor  band  at  micro- 
wave  frequencies  provides  the  measure  of  the  last  term.  This  fact  provides  the 
basis  for  the  development  of  a  radiometric  line-integral  refractometer . 

D.  Summary 

There  is  a  close  relationship  between  absorption  and  refraction  of  electro¬ 
magnetic  energy  because  both  phenomena  originate  in  the  molecular  interaction  of 
the  atmospheric  gases  with  an  impressed  electromagnetic  field.  That  is,  both 
absorption  and  refraction  proceed  in  proportion  to  atmospheric  composition.  Also, 
the  relationships  are  such  that  at  optical  frequencies  (wavelengths  less  than  20 
microns),  the  form  of  the  atmospheric  input  parameter  to  absorption  calculat'''ns 
is  precisely  that  of  the  electronic-charge  distortion  terms  of  line- integrated 
refractivity;  at  microwave  frequencies  the  absorption  parameter  form  matciies 
that  of  the  N  dipole  orientation  term.  This  fact,  plus  the  fact  that  the  dry-air 
constituents  exist  in  remarkably  crmstant  proportion  (so  that  knowledge 'of  the 
density  of  one  such  constituent  constitutes  such  knowledge  of  al'  the  rest)  provides 
the  basis  for  the  development  of  a  radiometric  line-integral  r  jfractometer.  For 
if  radiation  transfer  measurements  (measures  of  absorbed  radia  jon  intensity) 
performed  in  oxygen  and  in  water- vapor  absorption  bands  at  near- infrared  fre- 
(luencies,  and  in  a  water-vapor  band  at  microwave  frequencies,  the  data  so  ac¬ 
quired  may  be  interpreted  to  determine  precisely  all  tliree  terms  of  the  line- 
integrated  refractivity  equation. 
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5.  ACTIVE  AND  PASSIVE  RADIOMETRIC  SYSTEMS 


In  the  approach  to  the  measurement  of  line -Integrated  refractivity  just  oui- 
lined,  two  operating  modes  r>ui  be  considered.  First,  assume  that  an  Intrinsic 
source  of  monochromatic  energy  was  available  and  that  its  apparent  intensity 
could  be  continuously  monitored.  Then,  considering  for  this  case  that  Jy  =  0**, 
Equation  59  reduces  to 


(83) 


Equation  83  is  readily  inverted  to  yield 

I 

~  Jtn  (84) 
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Thus,  the  integral  of  the  absorption  coeificient  is  simply  given  as  the  negative, 
natural  logarithm  of  the  ratio  of  the  received-to-emitted-sourco  energy  intensity. 
This  is  termed  the  active  system  mode  of  operation.  The  instrument  to  determine 

n  n 

the  desired  ratio  would  be  termea  a  transmissometer. 

For  the  second  case,  suppose  an  intrinsic  source  was  not  available  for  moni¬ 
toring.  Then  measurements  of  would  depc>nd  solely  on  measurements  of  atmos¬ 
pheric  thermal  radiation.  In  this  case  Equation  59  reduces  to 


"  \  ®  (85) 

An  explicit  inversion  of  Equation  85  is  not  possible,  so  a  solution  for  the  integral 
of  the  absorption  coefficient  must  be  obtained  by  numerical  techniques.  This  is 
the  passive -system  mode  of  operation,  with  the  Instrument  termed  a  radiometer. 

*  Tills  condition  can  be  realized  in  one  of  two  ways:  (1)  specify  a  frequency  at 
which  the  thermal  emission  of  the  atmosphere  is  negligible,  or  (2)  modulate 
the  Eource  aiul  employ  a  background  rejection  mode  of  operation. 
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In  tracking  operations,  what  is  desired  is  the  measure  of  refractivity  inte¬ 
grated  precisely  over  the  path  from  target  to  receiver.  With  a  cooperative  target, 
such  a  measure  may  be  obtained  by  moonting  a  radiation  beacon  or  beacon  trans¬ 
ponder  on  the  target  itself  and  monitoring  continuously  the  apparent  beacon  inten¬ 
sity  as  the  target  is  tracked.  This  mode  of  operation  will  provide  the  greatest 
precision  in  N-measurement  and,  hence,  in  i  ^nge  parameter  measurement.  The 
active -system  approach  is  the  desired  approach,  and  it  is  applicable  to  operation 
at  missile  test  ranges. 

Suppose,  however,  noncooperative  targets*  are  to  be  Ira'^ked.  Here,  the 
passive-system  mode  of  operation  is  mandatory.  While  such  targets  are  beyond 
the  atmosphere,  such  a  system  provides  an  accurate  N-mcast  'eme  t  because  the 
passive  system  will  sense  the  emitted  radiation  contributed  by  ti.e  enure  deptn  of 
the  atmosphere.  Howe\  r,  for  targets  imbedded  within  the  atmosphere,  the  'iccu- 
racy  is  degraded  because  only  the  contribution  between  target  and  -ec''ivL  r  is  ihe 
governing  parameter  for  the  range  corrections.  A  variable  ^requ  y  radiometer 
could,  a '‘ter  the  method  ci  Westwater  (Reference  26),  provide  a  capa,oility  <'or  a 
degree  of  proflLag.  hov  jver,  so  this  error-source  could  son  ewha  min  mize. 
Residual  N  errors  (he  .ce,  residual  range -errors)  would  still  be  g.  eater  th  those 
characteristic  of  ue  ctive  system. 

So  the  set  of  circum  itances  is  as  follows.  A  desire  for  ht  greatest  ac>-  uracy 
demands  an  active  system,  but  expediency  may  require  the  passive  system.  Thero- 
fore,  the  ideal  radiometric  line-integral  refractometer  ought  to  bt  one  tha'  is 
operable  in  either  the  active  or  passive  mode  — the  active  mode  for  the  most  pre¬ 
cise  measurements  with  cooperative  targets,  the  passive  mode  such  Uat  the  meas¬ 
urement  capability  may  be  extended  to  the  case  of  non.. . (Operative  targeis. 

A.  Active  N-meaauremant  system 

Regarding  the  active-system  mode  ot  operation,  ui  >  measure  oi 
required  for  the  determination  of  line-integrated  refractlvit',  is  usually  expressed 
in  terms  of  the  transmittance,  which  is  defined  siinply  us  the  i  iilo  of  the  received- 
to-emltted  energy: 
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The  complement  of  transmittance  is  termed  absorptance: 

(87) 

Both  transmittance  and  absorptance  are  dimensionless  quantities  and  are  of 
magnitudes  that  can  vary  betw'een  zero  and  unity. 

The  absorption  process  is  not  the  only  process  that  acts  to  attenuate  a 
beam  of  radiant  energy  that  traverses  the  atmospheric  medium ,  ^Hirther  atten¬ 
uation  can  come  about  as  a  result  of  a  redirection  of  a  portion  of  that  energy, 
which  is  due  to  deformation  of  energy  wave  fronts  throuj^  encounters  with  atmos¬ 
pheric  suspensolds.  This  mechanism  of  radiation  attenuation  is  termed  scatter¬ 
ing  attentuation.  Generally,  it  also  is  a  frequency -dependent  process,  but  in 
contrast  to  the  absorption  process,  the  dependence  is  much  weaker  and  is 
monotonic  in  frequency.  The  picture  that  emerges  then  is  that  the  band  spectrum 
attributable  to  the  absorption  process,  and  which  provides  the  basis  for  the  radio- 
metric  measure  of  line -integrated  refractivity,  is  superimposed  upon  a  continuum 
spectrum  that  is  largely  attributable  to  the  scattering  process  but  which  also  con¬ 
tains  a  contribution  representing  the  effect  of  summation  of  the  wings  of  thousands 
of  absorption  lines  occurring  in  the  "selective"  absorption.  Obviously,  some 
account  must  be  taken  ol  tnis  continuum  attenuation  if  the  measure  of  transmittance 
is  to  be  interpreted  in  terms  of  the  integrals  of  the  atmospheric -constituent  partial 
'  tensities  required  to  define  line -integrated  refractivity. 

(1)  Need  for  a  reference  bund 

The  basis  for  an  accounting  for  the  continuum  attenuation  is  provided 
by  the  fact  that  the  processes  of  absorption  and  scattering  are  linear.  At  any  ra¬ 
diation  fretiuency  then,  the  total  attenuation,  or  extinction,  can  be  written  as  the 
sum  of  the  separate  absorption  and  scattering  attenuations  (Reference  23); 


Extinction  ^  Absorption  ♦  Scattering 


l«S) 


As  a  consequence  of  the  above,  transmittance  demonstrates  a  multiplication  prop¬ 
erty  such  that  total  transmittance  repret  ents  the  product  of  the  component  trans- 
mittances.  That  is. 
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where  the  subscripts  T,  a,  and  s  demark  the  total,  absorption,  and  scattering 
transmittances ,  respectively. 

Assume  that  radiatiem  intensity  measurements  were  perfornicd  not 
only  in  the  required  absorption  band  but  also  in  an  adjacent  band  lying  in  an  atmos¬ 
pheric  "window. "  From  Equations  86  and  8S  it  can  be  seen  that  the  absorption 
band  intensity  would  be  related  to  the  intrinsic  source  energy  according  to 
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where  the  subscript  "sens"  denotes  the  absorption -band,  or  sensing-band,  energy. 
The  window,  or  reference -band,  energy  intensity  would  suffer  rttenuation  due  only 
to  the  scattering  process.  Therefore, 


(91) 


The  ratio  of  the  sensing-to-reference-band  energy  intensities  would  then  be 


(92) 


If  the  scattering  transmittances  at  the  sensing  and  reference  bands  can  be  assunied 
(.Hjual,  then  Equation  92  reduces  to 
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Therefore, 
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(94) 


where  K,  the  ratio  ui  the  intrinsic  reference -fjand  energ>'  to  the  intrinsic  sensing- 
band  energy  can  be  treated  as  a  constant  parameter  of  the  measurement  system. 
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The  residual  of  the  received  sensing -to-reference-band  energy  ratio  is  thus  a 
function  basically  of  the  desired  absorption  transmittance.  The  technique  of  per¬ 
forming  the  above -defined  band  energy-ratio  measurement  not  jnly  results  in  the 
acquisition  of  the  direct  measure  of  the  absorption  transmittanc .  but  also  elimi¬ 
nates  the  need  for  a  hi^ly  calibrated  intrinsic  source.  The  only  source  require¬ 
ment  is  stability. 

(2)  Selection  of  sensing-  and  reference -band  frequencies  — clear- 
weather  system 

The  entire  electromagnetic  spectrum  —  from  the  ultraviolet  through 
tlis  radio  regions  —  is  rich  in  atmospheric  absorption  bands  that  are  potentially 
useful  for  the  development  of  a  spoctrographic ,  differential -absorption,  line- 
integral  refractometer.  The  problem  of  which  particular  bands  to  use  as  the  basis 
for  instrument  development  immediately  presents  itself.  And  arbitrary  as  the 
selection  mi^t  at  first  seem,  several  considerations  will  narrow  the  choice. 

First,  there  is  the  criterion  for  selection  provided  by  Equation  81.  That  is,  for 
tht  most  precise  measurement  of  line -integrated  refractivity ,  the  first  two  terms 
in  the  N  equation  should  be  determined  from  spectrographic  measurem'^nts  in  the 
"optical"  portion  of  the  electromagnetic  spectrum,  the  third  term  should  be  ueter- 
mined  from  similar  me-asuremenfs  performed  in  the  millimeter  region. 

(A)  "Optical"  absorption  bunds 

The  frequency  dtuondenceof  scattering  attenuation  is  extremely 
complex  and  wiU  be  thoroughly  discussed  in  later  sections.  present,  a  general 
laile  will  suffice;  For  given  particle  size  di.stributions ,  Lae  higher  the  frequency 
the  more  efficient  the  scatter.  H.e  lower-frequency  (longer  wavelength)  "optical" 

ridiatlon  thus  has  the  greater  anility  to  penetrate  che  «;tnor.pher Ic  arronol , 
thereby  increasing  the  effective  range  of  any  measurement  system.  Accordingly,  the 
chosen  absorption  bands  should  be  on  the  low-fre<iuency  end  of  the  optical  spectrum. 
Consider,  therefore,  the  solar  ab8orptlt*n  spectrum  normalized  to  (Mjak  transmis¬ 
sion  lha*  has  been  presented  by  Gales  and  is  shown  in  Figure  7.  nie  spt'clrum  ui 
Figure  7  is  characterized  by  four  strong -overtone  liaads  of  the  vibrabon-roLation 
water-vapor  spectrum,  two  weak-overtone  bands  of  the  carbon  lUv^xide  spectrum, 
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ami  one  weak  band  of  the  0-0  rotational  oxygen  spectrum.  Also  indicated  are  four 
window  regions.  Radiation  in  any  of  the  windows  eould  be  monitored  to  supply  the 
reference  signal  that  in  any  of  the  water-vapor  bands  is  the  sensor  water-vapor 
signal. 

Aceordingly,  th.ee  narrow -wavelength  bands  satisfying  the  re¬ 
quirements  for  recovery  of  the  first  two  terms  of  the  N  equation  could  be  centered 
at  wavelengths  of  1.22y,  l,32y,  and  1.  264y  for  the  reference,  water-vapor,  and 
o.xygen  or  dry-air  density  sampling  bands,  respectively.  The  water-vapor  sampl¬ 
ing  band  is  displaced  from  the  center  to  the  wing  of  the  strong  1.  37  absorption 
band  because,  over  any  sort  of  range  at  all,  virtually  all  the  1.37  radiation  would 
be  aJjsorbed  from  an  infrared  beam  and  would  leave  no  signal  to  effect  the  desired 
water -vapor  measurement. 

(B)  Microwave  absorption  bands 

In  the  microwave  absorption  spectrum  (Figure  8)  both  oxygen 
and  water  vapor  are  evident  absorbers,  .'ater-vapor  absorption  is  due  to  a  pure 
rotational  transition  at  22. 2  GHz  (1.35  cm)  and  the  nonresonant  contribution  of  other 
rotational  lines  starting  at  around  180  GHz  (0. 166  cm)  and  extending  w'ell  into  the 
infrared  (Reference  28).  Oxygen  absorption  is  due  to  a  band  of  rotational  transitions 
that  start  at  about  .53  GHz  (0.  57  cm)  and  extend  to  about  66  GHz  (0.  45  cm)  (Refer¬ 
ence  29).  The  dominance  of  the  oxygen  over  the  water -vapor  absorption  occurs 
because  of  the  large  density  ratio  of  ojgrgen  to  water  vapor,  since  the  mass  absorption 
coefficients  of  the  two  gases  are  of  the  same  order  of  magnitude.  If,  however,  the 
20-GHz  water-vapor  band  is  selected  as  the  sensing  band  for  the  third  term  in  the 
N- equation,  the  oxygen -absorption  contribution  to  the  received  signal  may  be  largely 
negated  by  referencing  the  20GHz  signal  to  a  reference  signal  in  the  water -vapor 
window  at  10  GHz,  Note  that  oxygen  absorption  at  these  two  frequencies  differs  in 
magnitude  only  slightly. 

(.3)  Selection  of  sen.sing-  and  reference -band  frequencies  — cloudy - 
weather  system 

(A)  Degradation  of  IR  measurements  in  clouds 

The  hjlirid  system  just  described  — the  combination  infrared 
and  microwave  active  system  — will  provide  the  greatest  accuracy  in  N 
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measurenienl.  For  missile  test  range  operations  (r-’''’racterized  by  (1)  the  avail¬ 
ability  of  cooperative  targets ,  and  (2)  conformance  of  test  schedules  to  ideal 
weather  conditions),  this  is  the  best  approach  to  real-time  measurement  of  N. 

The  greatest  drawback  to  the  system  is  the  inability  of  the  specified  infrared  wave¬ 
lengths  to  penetrate  appreciable  cloud  depths.  The  scattering  process  is  the  domi¬ 
nant  mechanism  of  extinction  of  infrared  energy  in  clouds.  Moreover,  tlie  radii 
of  the  drops  that  nominate  the  scattering  and  extinction  characteristics  of  clouds 
are  in  the  range  of  5  to  20^  (Reference  30).  This  means  that  the  Mie  scattering- 
area  ratio  appropriate  to  cloud  drops  can  be  approximated  by  its  asymptotic  value 
K  (a/x)  =  constant  =  2  (Reference  301,  The  scattering  coefficient  for  clouds  can, 
therefore,  be  expressed  as 

~  2.  'H  (X.  r\  (95) 

where  a  is  the  median  droplet  radius,  n  the  droplet-number  concentration.  The 
nonselective  scattering  coefficient  implied  by  Equation  95,  when  applied  to  the 
scattering  process  in  clouds,  is  confirmed  by  the  fact  that  light  extinction  in  clouds 
is  shown  experimentally  to  be  virtually  constant  throughout  the  ultraviolet,  visual, 
and  near-infrared  regions  of  the  electromagnetic  spectrum  (Re.  .rence  30). 

The  two  basic  parameters  required  for  estimating  the  extinction 
coefficients  of  clouds,  cloud-droplet  concentra'ion  and  median  droplet  size,  are 
functions  both  of  cloud  type  and  of  cloud  origin  (Reference  31).  In  general,  layer 
clouds  (clouds  of  tlie  genus  "stratus")  tend  to  have  lower  droplet  concentrations 
than  do  the  fair-weather  convective -type  (cumulus)  clouds.  The  same  relationship 
is  true  for  clouds  of  maritime  ''ersus  continental  origin.  This  latter  effect  alone 
may  result  in  a  difference  in  droplet  concentration  by  a  factor  of  five  in  clouds  of 
the  same  genera  (Page  112  of  Reference  31).  A  further  generalization  is  that 
clouds  with  large  droplet  concentrations  usually  consist  of  small  droplets,  whereas 
those  with  small  concentrations  contain  large  drops.  Thus,  layer  clouds  contain 
droplets  larger  than  those  of  the  more  active  convective  clouds,  and  droplets  of 
the  maritime  clouds  are  larger  than  those  of  the  continental  clouds.  This  general¬ 
ization  is  borne  out  by  measurements  (Reference  31). 


High-altitude  clouds,  the  cirrus  clouds,  and  the  tops  of  some 
cumulus  clouds  are  composed  of  ice  crystals.  Natural  ice  crysta.  j  assume  a 
variety  of  shapes  (Reference  32).  The  predominant  crystal  form  in  high -altitude 
clouds  is  the  hexagonal  column,  which  is  typically  500^  long  with  a  length-to- 
breadth  ratio  varying  from  one  to  five.  Although  such  ice  crystals  can  hardly  be 
considered  spheres,  scattering  from  long,  circular  cylinders  (which  should  be 
similar  to  cloud  ice  crystal  scattering)  has  been  shown  to  be  remarkably  similar 
to  that  from  spheres  (Reference  30).  The  radius  of  the  cylindrical  cross  section 
is  the  important  dimension  of  the  cylinders  as  far  as  scattering  is  concerned. 
Thus,  within  the  limits  of  knowledge  of  the  physics  of  ice  clouds,  Equation  95  may 
be  applied  to  estimate  the  scattering  coefficient  of  ice  clouds  as  well  as  water 
clouds,  if  "a"  is  taken  as  the  radius  of  the  ice  columns.  The  average  value  of  the 
parameters  required  for  the  specification  of  cloud-scattering  attenuation  coeffi¬ 
cients  of  various  cloud  types  and  for  fog  is  summarized  in  Table  1. 


Table  1:  Cloud  Particle  Sizes  and  Number  Densities 


Cloud  Type 

Median  Particle 
Concentration 
(cm -3) 

Median  Particle 
Radius 
(microns) 

Fair-Weather  Cumulus 

200 

6 

Cumulus  Congestus 

70 

15 

Stratus 

80 

9 

Cirrus 

5 

30 

Fog 

1C 

25 

According  to  Equation  95,  the  scattering  attenuation  coefficients  applicable  to  the 
various  cloud  types,  as  specified  by  the  pax’ameter  values  given  in  Table  1,  are 
as  summarized  in  Table  2. 
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Table  2:  Scattering  Coefficients  of  Clouds 


(km'^  of  cloud  depth) 


Fair-Weather  Cumulus 
Cumulus  Congestus 
Stratus 
Cirrus 


The  -values  summarized  in  Table  2  agree  well  with  measured  values  obtained 
by  Zabrodskii  and  Morachevskii  (Reference  33).  It  must  be  remarked,  however, 
that  the/jv  -values  shown  represent  more  or  less  ''average"  conditions;  i.e.  , 
variations  of  plus  or  minus  a  factor  of  10  from  cloud  to  cloud  on  a  day-to-day 
basis  are  the  rule  rather  than  the  exception. 

Defining  an  opaque  cloud  as  one  in  which  the  transmittance  is 
0.01  or  less,  as  suggested  py  Gates  and  Shaw  (Reference  34),  thicknesses  of  the 
above  cloud  types  required  to  reduce  transmittances  to  this  point  are  shown  in 
Table  3. 

Table  3:  Cloud  Thicknesses  Required  For 
Transmittances  Equal  To  0.01 

Cloud  Thickness 

Cloud  Type  (meters) 


Fair-Weather  Cumulus 
Cumulus  Congestus 


Stratus 


Cirrus 


The  cloud  thicknesses  tabulated  above  represent  nominal  thicknesses.  Cloud 
droplet  size  distributions  and  concentrations  vary  widely  from  day  to  day  and  from 
location  to  location,  and  so  do  the  cloud  thicknesses  required  for  opacity.  But 
despite  this  variation,  it  is  still  safe  to  say  that  any  cloud  will  be  likely  to  negate 
the  performance  of  the  infrared  transmittance  measurements  of  the  hyfjrid,  active 
system. 


(1)  Use  of  rnicrowave  frequencies  in  clouds 

To  extend  the  active-system  capability  to  operation  under  cloudy 
conditions,  microwave  techniques  must  also  be  applied  to  the  determination  of  the 
first  two  terms  in  the  N  equation.  The  previously  specified  22-GHz  water-vapor 
band  will  suffice  to  determine  the  second  term,  aial  the  first  term  may  be  recovered 
from  transmittance  measurements  in  the  60-GHz  oxygen  band  (Figure  8).  The  36- 
GHz  window  may  be  used  as  the  reference  for  this  latter  measurement.  However, 
in  both  these  measurements  for  the  ^  term,  ^  the  term) 

an  uncertainty  is  introduced  into  the  final  estimate  because  the  absorption  coeffi¬ 
cients  at  22  and  60  GHz  no  longer  exhibit  the  proper  forms.  That  is,  what  is  de¬ 
sired  is  the  integrals  of  the  oxygen  and  water-vapor  partial  densities  alone,  and 
not  the  integrals  of  the  densities-to~temperature  ratios. 

But  consider  the  integral 


5 


Pa. 


The  equation  of  state  relates  partial  density  to  pressure  and  temperature  accord¬ 
ing  to 


i'A 


A* 


\ 


(96) 


where  is  the  molecular  weight  of  the  absorbing  specie,  R*  is  the  universal 
gas  constant.  The  hydrostatic  equation  further  states  that 
(ip  -  - 

=.  -  : -  . 

where  g  is  the  acceleration  of  gravity,  is  the  mixing  ratio  of  the  absorbing 
gas  of  partial  density  ^  ,  Finally,  the  differential  path  lengi;h  element  is  related 
to  the  vertical  element  according  to 

as  -PtOJi  (9H) 

where,  for  elevat'on  angles,  C-  ,  greater  than  about  10  degrees,  f  Ct)  is  very 
nearly  equal  to  cst  4  ,  At  lower  elevation  angles,  (  (*  )  takes  on  a  greater 


47 


ATTENUATION  (db  /  km) 


A=10  4  2  1.5  1.0  0.8  0.6  0.4  0.2  0.15  0.1cm 


Figure  8;  THEORETICAL  VALUES  OF  ATMOSPHERIC  ATTENUATION  BY 
OXYGEN  AND  UNCONDENSED  WATER  VAPOR  AT  SEA  LEVEL  FC 
A  TEMPERATURE  OF  20°C.  THE  SOLID.  CURVE  GIVES  THE  ATTENUATION 
BY  WATER  IN  AN  /TM05PHERE  CONTAINING  1 -PERCENT  WATER 
MOLECULES  (p=7.5  g/m^.  THE  DASHED  CURVE  IS  THE  ATTENUATION 
BY  OXYGEN  (REFERENCE  24). 
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complexity  due  to  the  curvature  of  the  Earth.  E'  uation  98  into  Equation  97  ami 
division  of  both  sides  of  the  equation  by  T  yields 


T  3 


if: 

T 


(99) 


Integrating  appropriately, 

'O 


\  <_l  )p 


T 


(100) 


Po 


The  integrations  are  to  be  performed  from  the  point  S  -  0  (where  the  pressure 
is  )  to  the  point  S  (where  the  pressure  is  )■ 

If  the  reciprocal  of  temperature  as  defined  by  Equation  96  is 
now  substituted  into  the  right-hand  side  of  Equation  100,  then 


Pi 


—  - 


u-». 


f^e) 


r  o 

The  law  of  the  mean  states  that 

Pi 


(101) 


v>  = 

la 


C  Pi 

J  P) 

-JLs _ 

^  d(>^  f  ) 
fo 


\ 


ci  l. 


(loJ) 


S 


d  e. 


where  the  brackets  denote  the  average  value  Therefore, 

'  A  d(z.,  iJd_  \  ^ 


(io:q 


The  integral  on  the  right-hand  side  of  Equation  103  is  that  of  the  right-hand  side 
Equfi.'on  101: 
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'a 


-y  \^p 

Y 


(101) 


Pc  '  f's 


The  term  (Z^  -  Z^.)  in  Equation  103  can  be  found  by  the  method 
used  in  meteorology'  to  find  pressure  altitude.  This  involves:  (1)  substitution  for 
density  from  Equation  9G  into  Equation  97,  and  (2)  assuming  linearity  in  temixjra- 
ture  so  that  the  mean  value  can  be  moved  through  tlie  integral.  Thus, 


_£p 


A 


A  z- 


and,  therefoj'e, 

LU^  t\ 

<  I ; 

Equfition  106  into  Equation  103  then  yields 


il^-A 

^n\  I 

\Ps  / 


F’c  ^ 


iik,  ..a.  \ 

<17 

From  Equation  101,  it  therefore  follows  that 


<  !  >  \^ 


^  ..  \  iLl 

’rj 


(103) 


(106) 


(107) 


(lOH) 


The  right-hand  side  of  Equ;ition  lOs  is  precisely  the  right-hand  side  of  Equation 

100.  Therefore, 

S 

/^ds  *■  \  ^  Az  (109) 

T 


A  C<= 
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or,  from  Equation  93 


ds 


<T>f(0 


(110) 


At  first  glance,  the  rather  elai. orate  derivation  Just  presented 
seems  almost  superfluous  because  the  final  result  (Equation  110)  seems  almost 
obvious.  However,  the  "mean"  temperature  must  be  carefully  defined  to  permit 
its  removal  from  the  integral  of  the  microwave  absorption  coefficient.  Specifi¬ 
cally,  the  mean  must  be  taken  w'ith  respect  to  the  natural  logarithm  of  pressure. 

Nevertheless,  on  the  basis  of  Equation  110,  Equation  SI  permits 
the  statements  that 


<T>-6(' 


L 

60  GHz. 


) 


(lil) 


Z2C.  Hz 


) 


(113) 


where  the  angled  brackets  now  denote  the  mean  process,  it  is  on  Uie  ba.sis  of  the 
above  equations  that  N  can  be  specified  completely  in  term.i  of  micrusvave  trans¬ 
mittance  measui'ements.  The  additional  uncertainty  (above  that  awUiinable  from 
iPifrared  measurements)  enters  because  of  the  need  to  estimate  the  additioivil 
jxirameter  <T>  .  Auxiliary  meteorological  date  (radiosondes)  are  thus  rt<juired 
to  complete  the  N  specification. 

B.  l  asslve  N -measurement  system 

teaching  cooperative  targets,  the  requirement  is  for  the  jussive  .\'- 
measurjmont  system.  The  integral  of  the  ab.sorption  coefficient  is  reeovereti 
from  the  inversion  of  the  second  term  "'n  the  right-hand  side  of  the  radiative  trans¬ 
fer  exjuation  (Equation  59).  In  this  case,  represents  :he  intensity  of  the  ther¬ 
mal  radiation  incident  on  a  radiometer  (a  sensitive  radiation  thermometer)  that 
originates  in  the  atmosphere  itself.  Bcsrausc  the  spectral  radiance  of  .such  emitted 
energ-,’  is  as  dictated  by  Planck's  black-or-gray-ljody  radiation  law,  the  lempc'ra- 
tures  tvpical  of  the  atmosphere  dictate  that  the  greater  contribution  cbmes  from 
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the  lower  frequency  (longer  wavelength)  radiation.  In  fact,  emitted  energies  are 
negligible  at  the  active-system  infrared  wavelengths  previously  specified.  There¬ 
fore,  to  ensure  measurable  signals  and  to  provide  the  greatest  degree  of  all- 
weather  capability,  consideration  should  be  given  to  developing  a  passive  N- 
measurement  system  operating  at  the  microwave  frequencies  in  the  absorption 
bands  previously  specified. 

(1)  Total  absorption  versus  frequency 

Because  of  the  finite  widths  and  pressure-broadened  line  shapes  of 
atmospheric  constituent  absorption  bands,  the  total  (total  throughout  the  depth  of 
the  atmosphere)  integrated  absorption  changes  rapidly  with  frequency  as  the  fre¬ 
quency  is  scanned  from  regions  outside  the  band  to  the  resonant -band  frequency. 
For  example,  for  frequencies  near  the  center  of  the  60-GHz  o.xygen  band,  the  ab¬ 
sorption  is  nearly  complete.  In  the  frequency  region  from  55  to  50  GHz,  the  ab¬ 
sorption  decreases  as  the  frequency  decreases  and  approaches  the  wings  of  the 
band.  Because  of  the  basic  connection  at  any  frequency  between  absorption  and 
emissivity,  this  effect  can  be  seen  in  the  plot  of  brightness  temperature  versus 
frequency  shown  in  Figure  9.  Brightness  temperature  is  related  to  accord¬ 
ing  to 


Near  the  center  of  the  band  the  atmosphere  is  essentially  a  blackbody  radiating 
at  a  temperature  nearly  equal  to  that  at  the  surface  of  the  Earth.  CH'er  the  rela¬ 
tively  small  frequency  interval  from  55  to  50  GHz,  a  brightness  change  of  about 
150°K  takes  place,  and,  in  the  steepest  portion  of  the  curve,  the  slope  is  about 
60°K  per  GHz.  The  data  contained  in  Figure  9  were  calculated,  according  to  Van 
Vleck's  oxygen-absorption  coefficient  from  actual  radiosonde  data  obtained  at 
Dakar,  Senegal  (Reference  27).  Similar  calculations  based  on  Gutnick's  mean 
midlatitude  water-vapor  profile  (Reference  35)  as  applied  to  a  standard  atmos¬ 
phere  are  shown  in  Figure  10.  Note  again,  the  rapid  change  in  brif^tness  tem¬ 
perature  in  the  interval  of  16  to  22  GHz.  The  slope  at  the  steepest  part  of  the 
curve  is  about  15 ®K  per  GHz.  The  difference  in  shapes  of  the  curves  shown  in 
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the  lower  frequency  (longer  wavelength)  rad’ation.  In  fact,  emitted  energies  are 
negligible  at  the  active-system  infrared  wavelengths  previously  specified.  There¬ 
fore,  to  ensure  measurable  signals  and  to  provide  the  greatest  degree  of  all- 
weather  capability,  consideration  should  be  t,  en  to  developing  a  passive  N- 
measurement  system  operating  at  the  microwave  frequencies  in  the  absorption 
bands  previously  specified, 

(1)  Total  absorption  versus  frequency 

Because  of  the  finite  widths  and  pressure-broadened  line  shapes  of 
atmospheric  constituent  absorption  bands,  the  total  (total  throughout  the  depth  of 
the  atmosphere)  integrated  absorption  changes  rapidly  witli  frequency  as  the  fre¬ 
quency  is  scanned  from  regions  outside  the  band  to  the  resonant -band  frequency. 
For  example,  for  frequencies  near  the  center  of  the  60-GHz  o.xygen  band,  the  ab¬ 
sorption  is  nearly  complete.  In  the  frequency  region  from  55  to  50  GHz,  the  ab¬ 
sorption  decreases  as  the  frequency  decreases  and  approaches  the  wings  of  the 
band.  Because  of  the  basic  connection  at  any  frequency  between  absorption  and 
emissivity,  this  effect  can  be  seen  in  the  plot  of  brightness  temperature  versus 
frequency  shown  in  Figure  9.  Brightness  temperature  is  related  to  I^  accord¬ 
ing  to 


Near  the  center  of  the  band  the  atmosphere  is  essentially  a  blackbody  radiating 
at  a  temperature  nearly  equal  to  that  at  the  surface  of  the  Ea^'th.  Ov’er  the  rela¬ 
tively  small  frequency  interval  from  55  to  50  GHz,  a  brightness  change  of  about 
150“K  takes  place,  and,  in  the  steepest  portion  of  the  curve,  the  slope  is  about 
60TC  per  GHz.  The  data  contained  in  Figure  9  were  calculated,  according  to  Van 
Vleck's  oxygen-absorption  coefficient  from  actual  radiosonde  data  obtained  at 
Dakar,  Senegal  (Reference  27).  Similai’  calculations  based  on  Gutnick's  mean 
midlatitude  water-vapor  profile  (Reference  35)  as  appHed  to  a  standard  atmos¬ 
phere  are  shown  in  Figure  10.  Note  again,  the  rapid  change  in  brightness  tem¬ 
perature  ill  the  interval  of  16  to  22  GHz.  The  slope  at  the  steepest  part  of  the 
curve  is  about  15®K  per  GHz.  The  difference  in  shapes  of  the  curves  shown  in 
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Figure  10:  BRIGHTNESS  TEMPERATURE  AS  A  FUNCTION  OF  FREQUENCY 
IN  THE  22  GHz  WATER- VAPOR  BAND 
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Figures  9  and  ICis  attributable  to  the  different  pressure  broadening  effects  active 
at  the  different  frequencies. 

The  suggestion  is  therefore  immediate  that  these  portions  of  the  oxy¬ 
gen  and  water-vapor  spectra  may  be  useful  for  a  ground-based  passive  probe  of 
the  lower  atmosphere  for  these  two  contituents.  This  is  so  because  the  large 
changes  in  brightness  temperature  over  the  indicated  frequency  intervals  satisfy 
the  dominant  requirement  for  ground-based  remote  probing  systems;  i.e. ,  the  \ 

changes  in  the  Tg  with  frequency  are  large  enough  to  permit  a  nondegenerate  | 

set  of  equations  to  be  obtained  for  inverting  the  bri^tness-temperature-versus-  | 

frequency  curve  to  determine  the  temperature -height  and  water -vapor -height  pro-  j 

files  of  the  lower  atmosphere. 

(2)  Inversion  of  the  radiation-transfer  equation 

In  terms  of  the  brightness  temperature  observed  when  looking  up 
into  the  atmosphere,  the  radiative-transfer  equation  (Equation  59)  for  the  passive 
mode  (I  .  =  0)  is  (Reference  36) 

Tg  ^  ds  (114) 

where  Tg  is  the  bri^dness  temperature,  T  the  absolute  temperature,  and 
can  be  either  the  oxygen  or  water-vapor  absorption  coefficient.  Equation  114  is  | 

nonlinear  and  difficult  to  solve  because  the  absorption  coefficient  that  occurs  in  | 

the  Kernel  is  a  function  of  path  through  its  dependence  on  temperature.  At  a  | 

single  frequency,  an  infinite  number  of  temperature  profiles  exist,  T(8),  so  that  the 
same  intensity  of  radiation  is  given  at  the  top  or  bottom  of  the  atmosphere.  Thus, 
the  recovery  of  a  unique  distribution  will  require  a  set  of  integral  equations;  r.s  j 

usual,  the  requirement  is  for  n  equations  to  define  n  variables  in  the  atmos¬ 
pheric  representation. 

A  description  of  the  atmosphere  under  this  requirement  has  been 
attempted  in  various  ways.  Becau  j  a  set  of  equations  like  Equation  114  is  mathe¬ 
matically  unstable  (Reference  37),  other  approaches  must  be  tried.  One  such 
approach  has  been  to  represent  T(8)  by  some  particular  analytic  function  and  then 
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to  use  techniques  based  on  the  technique  of  Lagrangian  multipliers,  such  as  a 
least-squares  fit  (Reference  38),  in  an  iterative  solution.  Some  examples  of 
functions  that  have  been  used  are:  (1)  truncated  algebraic  and  trigonometric 
series  (References  39  and  40),  (2)  exponential  (Reference  41),  and  (3)  Tschebyscheff 
polynomials  (Reference  27).  Thus,  in  an  algebraic  function,  the  series 

T  -  *■  AjZ  i-  ,  (115a) 

would  be  substituted  into  Equation  114,  and  the  coefficients  Aq,  Aj,  Ag,  ... 
would  be  determined.  Similarly,  for  the  trigonometric  series: 

T  ^  Aq  A I  s' 2.  +•  2.Z  -t-  ,  .  ^  (115b) 

fc-r  Tschebyscheff  polynomials: 

T  ■=A^+A/^^C2)  t-  .  .  (115c) 

The  accuracy  of  replication  of  the  atmospheric  temperature  struc¬ 
ture  is  limited  by  the  number  of  terms  in  the  series  expansion  and  the  complexity 
of  the  atmospheric  structure.  For  example,  a  simple  atmosphere  with  a  constant 
lapse  rate  can  easily  be  represented  by  the  first  two  terms  in  Equation  115a: 

TU)  -  t-A,z 

Here,  Aj  can  be  identified  with  the  lapse  rate  of  temperature,  whereas  T^  is 
the  ground  temperature. 

When  the  atmospheric  structure  is  more  complex,  a  larger  number 
of  terms  must  be  used  in  the  series.  This  leads  to  the  requirements  for  a  large 
number  of  measurements. 

(A)  Oxygen  profile 

From  Equation  66b,  the  oxygen  microwave  absorption  coef¬ 
ficient  can  be  expressed  as  a  function  of  pressure,  temperature,  and  frequency 
as  follows: 

2 

^0  ^  ^  (116) 

2  '  N  ,1 

o«k4 

hi  it 
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The  terms  can  be  defined  as 

^  L8S73  tib/kni 

—  .  /k»« 

for  :icentratlons  of  oxygen  in  the  atmosphere.  T  is  the  absolute  temperature 
in  degrees  Kelvin,  p  is  pressure  in  millibars,  and  »^is  the  frequency  in  GHz. 
The  summation  is  to  be  performed  for  the  odd  integers  from  1  to  45.  For  each 
integer  there  are  two  associated  resonant  frequencies  (Table  4). 

Table  4:  Resonant  Oxygen-Transition  Frequencies 


N, 

N+.  GHz 

N-,  GHz 

1 

56.2648 

118.7505 

3 

58.4466 

62.4863 

5 

59.5910 

6). 3061 

7 

60.4348 

59.1642 

9 

61.1506 

58.3239 

11 

61.8002 

57.6125 

13 

62.4112 

56.9682 

15 

62.9980 

56.3634 

17 

63.5685 

55.7839 

19 

64.1272 

55.2214 

21 

64.6779 

54.6728 

23 

65.2240 

54.1294 

25 

65.7626 

53.5960 

27 

66.2978 

53.0695 

29 

66.8313 

52.5458 

31 

67,3627 

52.0259 

33 

67.8923 

51.5091 

35 

68.4205 

50. 9949 

37 

68. 9478 

50.4830 

39 

69.4741 

49.9730 

41 

70.0000 

49.4648 

43 

70.5249 

48.9582 

45 

71.0497 

48.4530 

These  resonant  frequencies,  first  predicted  by  Van  Vleck,  were 
subsequently  measured  up  to  N  =  25.  At  each  Integer  quantum  number,  the  reson¬ 
ant  frequencies  have  been  calculated  up  to  N  =  45  by  Miller  and  Townes  (Reference 
42)  and  by  Tsao  and  Cumutte,  (Reference  43).  Excellent  agreement  has  been 
observed  between  the  calculated  and  observed  resonant  frequencies.  For  example, 
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the  standard  deviation  of  the  difference  between  the  observed  and  calculated 
frequencies  for  Zummerer  and  Mizushima's  work  is  1.35  x  10"^  GHz  out  of 
approximately  60  GHz  (Reference  29). 

The  term  5^  (consisting  of  terms  involving  ^  the  line  width 
parameter,  the  magnetic  dipole  moment  of  the  Nth  transition,  and  the  structure 
function  shown  in  Equation  66h^is  given  by 


r:  _  _ 

f  -  Az^ _ _ 

P  _  A  2/ 

V  M  +- 1 

Aiy.  ~  jq  ^ 


A2^ _ 


-  zO^ 4-r>(  +  i)^2N  t-i) 


The  line-width  factor,  Az/,  is  generally  assumed  to  be  propor¬ 
tional  to  pressure  and  temperature.  Meeks  and  Lllley  (Reference  44)  give  an 
expression:  ^ 

Ax^(  ^  a  p  21  +-0.74?k»J^^j  (118) 

where  a  and  b  are  constants.  The  factors  0.21  wd  0.78  will  be  recognized  as 
the  percents  by  volume  of  oxygen  and  nitrogen  in  the  atmosphere.  The  constant 
a  is  approximately  1.4626  MH*.  prr  ^UThe  "constant"  b  apparently  has  two 
values:  In  the  hlg^-pressure  range  at  pressures  greater  than  350  millibars, 
b  =  0.25;  at  pressures  less  than  25  millibars,  b  =  0.75.  Between  these  two 
pressures  there  is  a  region  of  uncertainty,  which  Meeks  and  Lllley  have  covered 
by  letting 
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h  -0.25 

where  h  is  the  altitude  of  computation,  Hj^  is  the  pressure  altitude  of  350  milli¬ 
bars,  and  H2  is  the  pressure  altitude  of  25  millibars.  Other  effects,  such  as 
Zeeman-splittir.g  and  Doppler-broadenlng,  occur  above  30  km  (lOmb). 

Since  the  recovery  of  information  about  the  atmospheric  struc¬ 
ture  remotely  depends  on  knowledge  of  its  precise  definition  is  a  matter  of 
some  importance.  The  simplest  method  of  accurately  determining  the  a  and  b 
parameters  would  be  with  concurrent  radiosonde -radiometer  measurements. 

Polynomial  approximations  to  the  atmosphere's  vertical  struc¬ 
ture  are  limited  in  accuracy.  Westwater  and  Abbott  (Reference  41)  have  given 
some  errors  that  can  be  associated  with  application  of  such  function  to  the  inver¬ 
sion  of  Equation  114.  Five  iterations  were  used  in  fitting  varying-degree  poly¬ 
nomials  to  Tg  data  obtained  at  12  measuring  frequencies.  Profile  errors  asso¬ 
ciated  with  polynomial  degree  are  given  in  Table  5. 

Table  5:  Errors  in  Profile  Parameters 
Versus  Polynomial  Degree 


Polynomial 

RMS  Error 

RMS  Error 

Degree 

(p,  mb  (fixed  hei 

3 

6 

3.6 

4 

4.2 

1.7 

5 

2.2 

-0.5 

A  further  difficulty  results  when  finite  random-brightness - 
temperature-measurement  errors  are  Introduced.  The  presence  of  such  errors: 

(1)  limits  the  number  of  terms  in  the  series  of  expansion,  and  (2)  introduces 
inaccuracies  In  the  derived  temperature.  Figure  11  is  from  Westwater  and 
Abbott's  error  analysis,  in  which  deviations  of  constant  magnitude  were  added  to 
the  postulated  measured  Tg's. 

The  actual  deviation  shewn  in  Figure  12  for  brightness-temperature 
errors  of  0*K,  ±1.25*K,  and±2,0*K  are  more  Instructive  as  to  the  actual  mathe¬ 
matical  processes,  in  terms  of  pressure  coordinates,  the  differential  equation 
that  generates  Equation  114  can  be  expressed  as 
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RMS  BRIGHTNESS  TEMPERATURE  MEASUREMENT  ERRORS  (®K) 


0  5  10  J  5  20  25 

RMS  KINETIC  TEMPERATURE  ERRORS  (®K) 


Figure  11:  RMS  ERROR  IN  KINETIC  TEMPERATURE  IN  THE  LAYER  0-10  km 
AS  RELATED  TO  RMS  MEASUREMENT  ERRORS  IN 
BRIGHTNESS  TEMPERATURE 
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Fioort  12:  DEVIATION  OF  CALCULATED  TEMPERATURE  PROFILE 
FROM  TRUE  PROFiU  FOR  VARIOUS  ALTERNATING 
ERRORS  IN  BRIGHTNESS  TEMPERATURE 
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If  it  is  assunied  that  T  is  some  function  of  p,  then  Equation  120  is  a  Riccati 
differential  equation  (Reference  45),  which  is  difficult  to  solve.  Each  Riccati- 
type  differential  equation  generates  a  different  class  of  functions  as  its  solution. 
The  most-studied  Riccati  equations  are  those  that  generate  such  functions  as 
exponential,  hyperbolic,  circular,  Legendrian,  and  Hermitian  types  as  their 
solution. 


Most  of  these  functions  are  periodic  or  nearly  periodic.  It  is 
likely  that  the  quasiperiodic  residual  errors  of  Westwater  and  Abbott's  results 
shown  in  Figure  12  arc  the  result  of  remainder  ter.ns.  These  residuals  are  the 
difference  between  the  4th-degree  polyTiomlal  approximation  and  the  true  function 
that  was  implicitly  generated  by  assuming  a  functional  relation  between  pressure 
and  temperature  of  the  form 


where  k  is  a  constant  and  r  is  a  variable  exponent. 

To  avoid  the  difficulties  associated  with  a  functional  repre¬ 
sentation  of  the  atmosphere,  the  model  shown  la  Figure  13  was  developed.  This 
model  postulates  mat  me  atinos|)here  can  r>e  divided  into  (n-1)  pressure  layers  that 
may  vary  In  thickness,  where  n  is  an  as  yet  untletermined  number.  Each  layer 
is  to  be  sufficiently  small  so  that  the  "mean”  thermodynamic  temperature  as  defined 


in  Equation  106  can  be  used  in  liquation  120  to  give  the  brightness  temperature 
at  the  bottom  of  the  layer  with  arbitrary  precision  as  compared  to  that  computed 
from  any  meteorologically  realizable  temporaturo-pressure  distribution. 

Under  this  restriction,  the  temperature  variable  in  Equation 
120  becomes  a  constant.  Terms  can  be  rearranged  and  integrated  for  the  layer 
bounded  by  pressures,  ip  ,f  i)  for  which  there  are  corresponding  brightness 
temperatures  a  ^ 


.?£i£,,\r  5  .  O'"' 


The  left-hand  side  of  Equation  121  can  be  integrated  iiL.-nedi- 


ately  to  give 


The  integral  equation  (121)  may  also  be  expressed  In 
closed  form  for  pressui'es  less  than  350  millibars,  where  b  is  a  constant  (see 
Equation  119),  since  given  by  Equation  118  can  be  expressed  as  an  integrable 
fuaciion  of  pressure  and  the  mean  temperature.  For  example,  consider  the  term 
(Equation  117b)  that  is  one  of  the  terms  that  make  up  Sj^.  For  simplicity, 


Ltj' ^  clLom 

The  integral  can  be  evaluated  directly  to  give 


■-f  r,  , 

/V  (n^  n~i) 


(125a) 


Similar  Integrations  can  be  performed  on  the  terms  Fj^  (Equations  117c  and 
117d).  Let  them  be 


dp  ■-  f . 


(125b) 


( r>  n-t) 


U.r^.O 


(125c) 


The  bracketed  term  on  the  right-hand  side  of  Equation  122  can  then  be  expressed 
as  the  logarithm  of  a  product  that,  since  it  is  an  exponential  tu.^  *:on,  becomes 
just  the  produLt  of  the  argument:  i.e. , 
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r  ^  f 

(r,(v  .|) 


-E. 


•  o,  ^ 


a26) 


The  expression  for  the  briefness  temperature  on  one  side  of 
a  layer  at  a  given  frequency,  in  terms  of  the  brightness  temperature  on  the  other 
side,  and  the  absorption  function  a^  at  the  mean  temperature  of  the  layer 
can  be  found  by  combining  Equations  126  and  122  to  give 


(127a) 


Note  that  in  the  above  equation  is  "basically"  the  parameter  desired  for 

estimation  of  the  contribution  of  oxygen  to  line-integrated  refractivity.  That  is, 
a^^  ^  represents  the  exponential  of  the  integral  of  the  oxygen  absorption  coefficient 
between  levels  n  and  n-1,  where  the  pressures  are  pn  and  . 

The  nomenclature  of  Equation  127  may  be  Illustrated  from 
Figure  13,  the  radiative  transfer  model.  For  instance,  !t  n  =  2,  then  Equation 
127  is 


T  ^  -  r 

_ ij.!. _ ~  ^  (127b) 

The  single  subscript  refers  to  a  specified  pressure;  for  Instance, 
Pj  is  the  surface  pressure  and  Pg  is  the  second  pressure.  Similarly  j  and 
Tjj  2  brig^itness  temperatures  at  pressures  Pj^  and  pg.  The  double  sub¬ 

script  refers  to  the  characteristics  of  a  whole  layer.  For  example,  T2  ^  and 
32  I  are  respectively  the  mean  temperature  and  the  absorption  function  bounded 
by  pressures  p^  and  P2. 

The  pressure  interval  (p^^  -  P^.j)  over  which  the  integration  is 
valid  has  been  ascertained  by  using  five  meteorolog^ically  realistic  atmospheres. 
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The  expression  for  the  bri^tness  temperature  on  one  side  of 
a  layer  at  a  gi^'sn  frequency,  m  terms  of  the  brightness  temperature  on  the  other 
side,  and  the  absorption  function  a^  at  the  mean  temperature  of  the  layer 
can  be  found  by  combining  Equations  126  and  122  to  gl\ 


- '5'e,o-i  ^  yt  (127a' 

Note  that  a^^  in  the  above  equation  is  ''basically"  the  parameter  desired  for 
estimation  of  the  contribution  of  oxygen  to  line-integrated  refractivity.  That  is, 

represents  the  exponential  of  the  integral  of  the  oxygen  absorption  coeffic’ent 


a 

n,  n-1 

between  levels  n  and  n-i,  where  the  pressures  are  pn  and 

The  nomenclature  of  Equation  127  may  be  illustrated  from 
Figure  13,  the  radiative  transfer  model.  For  instance,  if  n  =  2,  then  Equation 
127  is 


<\.> 


3  1 


—  a. 


-2,1 


(127b) 


■  '  3,2 

The  single  subscript  refers  to  a  specified  pressure;  for  instance, 
Pj  is  the  surface  pressure  and  Pg  is  the  second  pressure.  Similarly  Tg  ^ 

Tg  2  brightness  temperatures  at  pressures  pj^  and  pg.  The  double  sub¬ 

script  refers  to  the  characteristics  of  a  whole  layer.  For  example,  Tg  i  and 
02  2  are  respectively  the  mean  temperature  and  the  absorption  function  bounded 
by  pressures  pj^  and  Pg. 

The  pressure  interval  (p^^  -  Pjj.j)  over  which  the  integration  is 
valid  has  been  ascertained  by  using  five  meteorolog(ically  realistic  atmospheres. 
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The  comparison  was  made  by  calculating  the  brightness  temperatures  from 

Equation  121  in  a  standard  computer  routine  with  those  calculated  from  Equation 

127,  The  calculations  from  Equation  121  were  baaed  on  a  Runga-Kutta  variable 

step  solution  to  the  differential  equation.  The  steps  varied  from  10“^  millibar 

to  1  millibar.  It  was  found  that  the  T„  could  be  calculated  with  an  error  less 

J3 

than  0.25“K  for  pressure  increments  of  300  millibars  or  less.  These  calcula¬ 
tions  were  restricted  to  frequencies  for  w^lch  the  absorption  and  emission  above 
300  millibars  was  negligible  to  eliminate  the  problem  of  the  variable  b  (Equation 
119).  This  problem  could  be  obviated  by  stepwise  integration  of  the  integral 
equation  (122)  to  compute  the  absorption  functions  directly.  For  purposes  of 
later  computations,  however,  the  partial  derivative  of  a^^  with  respect  to 
Tn  is  required,  and  it  was  desirable  to  perform  the  integration.  Further, 
profiling  to  300  millibars  was  apparently  as  hi^  as  was  feasible,  because  of 
the  masking  by  the  absorption  in  the  atmosphere  below  300  millibars  of  the  radi¬ 
ation  from  levels  above  this  pressure.  This  results  in  a  great  loss  of  sensitivity. 

The  absorption  function  for  various  mean  temperatures  is 
shown  in  Figures  ll*,  15  and  l6.  These  data,  in  conjunction  with  measured 
brightness  temperatures,  permit  recovery  of  integrated  oxygen-absorption 
coefficients.  The  method  is  as  follows  (illustrated  for  a  three-layer 
atmosphere;  the  number  of  layers  is  determined  by  the  number  of  operating 
radiometer  frequencies). 

At  a  given  irequency,  the  briefness  temperature  at  the  bottom 
of  any  layer  can  be  expressed,  from  Equation  127  in  terms  of  brightness  temper¬ 
ature  at  the  top  of  the  layer,  the  mean  temperature  through  the  layer,  and  the 
absorption  function  through  the  layer.  For  a  three-layer  atmosphere,  the  three 


equations  required  to  do  so  are 

Tg,, 

a) 

+ 

''5,5 

4  (128) 

'ai  “  L' 

1  f 

c) 
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Figure  14:  ABSORPTION  FUNCTION  VERSUS  MEAN  TEMPERATURE 
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Tq  2  3  ^  eliminated  from  the  above  simultaneous  equations  to  give 


+  <S  ,■)[,- (a,,,)] <r,y  <T^,> 


This  equation  can  be  explained  in  physical  terms.  The  factor 
I  )(tCj  3^  ^f54  ^  brightness  temperature  at  the  bottom  of  the 

model  Td  i  contributed  by  the  flux  at  the  top  of  the  model,  \vhich  is  attenuated 

Ji,  1 

by  the  three  intervening  layers.  If  there  is  no  source  (such  as  the  Sun  or  Moon) 
in  the  antenna  pattern,  this  brightness  temp)erature  would  be  on  the  order  of  0” 
Kelvin.  The  term  [  1  ^  Is  the  self-emission  (since  1  minus 

absorptivity  is  equal  to  the  emlssivlty)  of  the  topmost  layer.  When  this  is  mul¬ 
tiplied  by  (<X3^X<a.^  I )  ,  the  attenuation  due  to  the  two  layers  below  it  gives  the 
contribution  to  Tg  fiom  the  topmost  layer.  Similar  attributes  can  be  assigned 
to  the  other  two  terms  in  the  equation. 

With  a  three-frequency  system,  three  equations  like  Equation 
129  differ  in  the<Tg>  coefficients  to  the  extent  the  a’s  are  frequency-dependent. 
These  three  equations  form  a  unique  set  of  transcendental  equations  that  must 
be  solved  by  numerical  methods.  An  iterative  approach  must  be  used  in  which 
the<lr3  ?  'b  for  each  layer  are  "gpuessed"  in  the  first  approximation.  Based  on 
these  gfuesses  's  the  appropriate  values  of  the  a's  are 
substituted  into  the  set  of  equations.  The  process  is  repeated,  narrowing  the 
choice  of  <T3> 's  until  the  matrix  of  the  set  of  equations  is  solved.  The  's 

so  defined  can  then  be  used,  in  conjunction  with  the  measured  brightness  temper¬ 
atures,  in  the  manner  presciibed  by  Ekiuation  127  to  determine  the  integrals  of  the 
absonXion  coefficients  for  each  layer. 

Strictly,  Equations  125a,  b,  and  c  should  also  be  integrated  with 
respect  to  frequency.  This  integration  can  be  accomplished.  However,  the  a's 
have  been  found  to  be  smooth,  almost  linear  functions  with  respect  to  frequency; 
thus  the  value  at  the  center  of  the  operating  band  can  be  used  because  it  closely 
represents  the  mean  value  that  would  result  from  the  integration.  Also,  this 
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integration  has  not  been  done  because  the  manner  in  which  the  integration  will 
have  to  be  performed  depends  on  the  filter  characteristics  of  the  electronics. 

Further,  the  antenna  beam  pattern  will  have  to  be  considered  because  this  factor 
will  somewhat  wei^t  the  results. 

The  frequencies  chosen  for  this  study  —  51.2,  52.  52.8,  and 

53.3  GHz  —  were  chosen  for  their  proximity  because  of  bandwidth  considerations 
in  the  radiometer  design.  They  are  centered  between  absorbing  lines  for  maxi¬ 
mum  smoothness  of  the  absorption  functions.  If  more  detailed  profiling  near  the 
ground  is  desired,  the  mean  temperatures  obtained  from  the  above  frequencies 
could  be  used  together  vith  frequencies  more  strongly  absorbed. 

(B)  Water-vapor  profile 

The  approach  to  surveying  the  water-vapor  profile  from  brightness- 
temperature  measurements  would  be  similar  to  that  of  the  cp®"  of  oxygen.  That 
is  inversion,  in  the  manner  Just  presented,  of  Equation  114  relative  to  the  22-GHz 
water-vapor  absorption  coefficient  would  permit  generation  of  a  ne’v  set  of  inte¬ 
grated  absorption  coefficients  as  functions  of  mean  temperature,  vleasured 
brif^tness  temperatures  performed  in  the  wing  of  the  22-GHz  bar  d  v  juid  then 
establish  the  set  of  equations,  represented  by  Equation  129,  required  to  define 
the  mean  temperature  in  anj*  atmosf^erlc  layer.  The  result  would  then  define 
the  integrated  absorption  coefficient  for  that  layer.  ; 

In  the  case  of  absorption  in  the  vicinity  of  22  GHz,  however,!  a 
complication  arises  because  the  wing  of  the  60-GHz  oxygen  band  extends  well 
into  the  22-GHz  frequency  region.  The  result  is  that  oxygen  absorption  can 
Influence  brightness  temperatures  in  the  22-GHz  vicinity.  For  large,  total, 
water-vapor  contents  (l.e.,  large  water-vapor  absorption),  such  as  would  be 
realized  in  looking  at  the  total  depth  of  the  atmosphere,  such  effects  ^re  negli¬ 
gible.  This  is  the  point  emphasized  by  Menius,  et  al,  in  their  discussion  of  the 
MARCOR  water-vapor  radiometer.  However,  if  consideration  is  to  be  given  to 
profiling  the  atmospheric  water-vapor  content  is  to  be  considered  (and  it  must  if 
passive  radiometric  techniques  are  to  be  applied  to  provide  slant-range  correc¬ 
tions  for  targets  Imbedded  in  the  atmosphere),  then  in  atmospheric  la^rs  of  low 
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absolute  humidity  the  oxygen  absorption  in  tae  vicinity  of  22  GHz  can  attain  signifi¬ 
cant  proportions.  This  fact  is  demonstrated  by  Figures  17  and  18,  which  show 
both  the  oxygen  and  water-vapor  absorpticm  (calculated  according  to  Equation  69), 
assuming  saturation  at  standard  surface  temperature  and  pressure  (Figure  17) 
and  saturation  at  stai^dard  300-mb  temperature  and  pressure  (Figure  18).  For 
example,  notice  that  under  surface  conditions  the  20-GHz  water-vapor  absorption 
is  about  an  order  of  magnit’'H(>  greater  than  the  ox}^en  absorption;  however,  at 
300  mb  the  water-vapor  absorption  ic  greater  than  that  of  oxygen  by  about  a  factor 
of  only  two. 

Nevertheless,  the  problem  seemingly  presented  here  is  readily 
solved  when  it  is  recognized  that  extinction  effects  are  additive  (Equation  88).  The 
integrated  layer  absorption  coefficients  derived  from  22-GHz  brightness -tempera¬ 
ture  measurements  can  thus  be  corrected  for  the  oxygen  effects  simply  by  sub¬ 
tracting  the  adjusted  (for  frequency)  absorption  coefficient  data  obtained  with  the 
55-GHz  profiling  radiometer.  Because  of  this  need,  it  is  safe  to  say  that  accurate 
water-vapor  profiling  cannot  be  accomplished  unless  corresponding  oxygen  data  is 
supplied  with  a  55-GHz  profiling  radiometer. 

(B)  Additional  comments 

It  is  often  stated  that  yie  dry -air  contribution  to  line -integrated 
refnictivity  can  be  attained  'Svlth  sufficient  accuracy  from  ground-level  measure¬ 
ments  and  available  knowledge  of  the  variation  of  atmospheric  density  with  alti  ¬ 
tude,  "  (Reference  25).  If  this  Is  true,  then  the  need  for  an  oxygen  radiometric 
system  is  questionable.  The  aee<'l  for  the  oxygen  measurement  capability  is  ade¬ 
quately  established  in  the  profiling  mode  of  operation;  but  beyond  that,  consider 
that  Equation  127,  when  the  entire  depth  of  the  atmosf^ere  is  considered,  is  essen¬ 
tially  that  presented  by  Menius,  et  al  (Reference  25),  for  analysis  of  the  MARCOR 
data.  Note  especially  that  equal  weights  are  given  to  both  the  "mean"  temperature 
and  the  brightness  temperature  for  their  effect  on  the  specification  of  integr;ited 
absorption  coefficient.  Therefore,  concern  over  accuracies  in  brightness - 
temperature  measurement,  as  related  to  accuracies  in  integrated  absorption 
specification,  is  just  half  the  problem.  Equal  concern  must  be  given  to  the  accu¬ 
rate  specification  of  t.T>.  Accuracy  in  <,T>  to  the  same  degree  as  the  accuracy 
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ABSORPTION  COEFUCIENT  (NEPERS/cm) 


FREQUENCY  (GHz) 


Figur*  17:  WATER-VAPOR  AND  OXYGEN  ABSORPTION  IN  A  SATURATED 

ATMOSPHERE  AT  STANDARD  SURFACE  TEMPERATURE  AND  PRESSURE 
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ABSORPTION  COEFFICIENT  (NEPERS/cm) 


Figur*  18:  WATER  VAPOR  AND  OXYGBM  ABSORPTION  IN  A  SATURATED 
ATMOSPH^tE  AT  300  mb  (Standard  T«fi>p«rotart) 


in  Tb  simply  cannot  be  attained,  ”from  ground-level  measurements,”  but  can 
only  be  attained  radiometrically  ^rom  oxygen  absorption  measurements  in  the 
manner  outlined  in  the  paragraph  entitled  ”The  Oxygen  Profile."  That  is,  the 
<T>  input  to  any  water-vapor  brightness-temperature  analysis  must  be  supplied 
from  the  data  attained  throu^  oxygen -absorption  measurements. 

C.  Summary 

Several  radiometric  refractometer  configurations  are  required  if  the 
need  for  N-measurement  is  to  be  satisfied  in  all  ranging  activities.  A  combina  - 
ticm  infrared/microwave  active  system  will  provide  the  most  accurate  N  data. 

This  is  the  system  best  suited  for  use  in  missile  test  range  operatlcms,  insofar 
as  such  operatioos  are  characterized  by:  (1)  the  availability  of  cooperative  tar¬ 
gets,  and  (2)  conformance  of  test  schedules  to  ideal  weather  conditions.  In  using 
the  1.26-yiw  oxygen  and  the  1.37-^  water -vapor  infrared  band,  and  the  22 -GHz 
microwave  water  band,  all  three  terms  of  the  integrated  Smith -Weintraub  refrac- 
tivity  equation  are  measured  in  their  precise  functic«ial  form . 

To  extend  the  operational  capability  of  the  active  system  to  cloudy  con¬ 
ditions,  microwave  absorption  in  the  55-GHz  oxygen  band  can  be  used  to  approxi¬ 
mate  the  dry -air  term  of  N,  and  the  22-GHz  band  can  be  used  to  provide  data 
relating  to  both  water-vapor  terms.  However,  in  this  approach  some  loss  of 
accuracy  is  incurred  because  an  independent  measure  of  a  suitably  defined  "mean" 
temperature  must  be  acquired  for  aa  adjustment  of  the  oxygen  and  water-vapor 
data.  That  is,  the  microwave  data  must  be  multiplied  by  this  mean  temperature 
before  tliat  data  can  he  accepted  as  defining  the  first  two  terms  of  the  N  equation. 

For  the  case  of  noncooperative  targets,  the  passive  mode  of  operation 
is  mandatory;  i.e. ,  N  must  be  determined  by  measuring  atmospheric  thermal 
emissions  In  the  microwave  oxygen  and  water-vapor  b.ands  previously  specified. 
Moreover,  multifrequency  capabilities  must  be  provided  in  each  band  to  permit 
a  degree  of  N  profiling  — a  mandatory  capability  if  the  passive  system  is  to  bo 
applied  to  tracking  noncooperaOve  targets  imbcdiled  in  the  atmosphere. 
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G.  ERROR  ANALYSIS 


To  estimate  residual  errors  that  remain  in  N  or  in  slant-range  data  after 
the  line-integral  refractometers  Just  discussed  are  applied  to  typical  ranging  op¬ 
erations,  certain  assumptions  must  be  made  regarding  the  atmosf^ere  and  the 
measurement  process  and  Instrumentation.  This  will  consist  largely  of  postula¬ 
tions  of  simplifi^u  models  for  the  refractometers  and  atmosphere,  which  will 
permit  at  least  a  semblance  of  analytical  procedure  in  an  otherwise  intractable 
problem.  The  final  results  can  be  given  no  more  of  an  interpretation  other  than 
that  they  represent  reasonable  estimates  of  what  actuality  will  eventually  show. 

A.  Basic  N-error  equation 

Denoting  Integrated  values  by  an  overbar.  Equation  79  is  written  as 

There  are  uncertainties  in  (lach  element  of  the  right-hand  side  of  Equation  130. 
The  k-constants  and  u  are  uncertain  to  the  extents  Indicated  in  Equations  74 
and  77,  respectively;  there  will  also  be,  as  In  any  physical  measurement,  un¬ 
certainties  in  the  specifications  and  ^/t).  If  all  these  uncertainties 

are  assumed  to  be  unoorrelated,  the  corresponding  uncertaicty  in  N  can  be  deter¬ 
mined  from  the  total  differential  of  N: 


4.^au;- -4-  — 

expressed  in  terms  of  finite  increments  as  \r/ 

=txAk,  *■  b/^k^rcAk3  -r 


(131a) 


(131b) 


In  the  above  Incremental  equation,  the  coefficients,  vblch  are  to  be  evaluated  with 
reference  to  some  set  of  "standard"  conditions,  are  defined  by 
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6.  ERROR  ANALYSIS 

To  estimate  residual  errors  that  remain  in  N  or  In  slant-range  data  after 
the  line-integral  refractometers  just  discussed  are  applied  to  typical  ranging  op¬ 
erations,  certain  assumptions  must  be  made  regarding  the  atmosphere  and  the 
measurement  process  and  insiromentation.  This  will  consist  largely  of  postula¬ 
tions  of  simplified  models  for  the  refractometers  and  atmosphere,  w4iich  will 
permit  at  least  a  semblance  of  anal3l:ical  procedure  in  an  otherwise  intractable 
problem.  The  final  results  can  be  given  no  more  of  an  interpretation  otlier  than 
that  they  represent  reasonable  estimates  of  what  actuality  will  eventually  show. 

A.  Basic  N-error  equation 

Denoting  integrated  values  by  an  overbar.  Equation  79  is  written  as 


^*0 


(130) 


There  are  uncertainties  in  each  elemeut  of  the  ri^t-hand  side  of  Equation  130. 
The  k-constants  ando  uncertain  to  the  extents  indicated  in  Equations  74 

and  77,  respectively;  there  will  also  be,  as  in  any  physical  measurement,  un¬ 
certainties  in  the  specifications  and  ^vTt).  If  all  these  uncertainties 

are  assumed  to  be  uncorrelated,  the  corresponding  uncertainty  in  N  can  be  deter¬ 
mined  from  the  total  differential  of  N: 


i  'c 


(131a) 


expressed  in  terms  of  finite  increments  as 


\r/ 

Abi  —  ^Ak,  ^  T  e 


(131b) 


In  the  above  incremental  equation,  the  coefficients,  vhleb  are  to  be  evaluated  with 
reference  to  some  set  of  "standard"  conditions,  are  defined  by 
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dwc 

JN  — 

dN 

^W3nr  j 

Q  d 

^2. 


(132) 


With  imcorrelated  uncertainties,  the  variance  in  N  is  simply 

In  the  above  equation,  notice  that  the  sum  of  the  terms  in  the  first  set  of 
brackets  on  the  ri^t-hand  side  represents  the  contribution  to  the  variance  in  N 
arising  from  basic  uncertainties  (the  Ak's  and  A ‘^(>2)  in  the  integrated  refractivity 
equation  itself.  As  such,  it  defines  the  fundamental  limit  of  the  accuracy  to  which 
N  may  be  specified  by  measuring  subsidiary  meteorological  parameters.  In 
other  words,  reduction  of  the  incremental  terms  in  the  second  bracket  by  using 
improved  sampling  techniques,  can  be  pushed  just  so  far.  Consider,  first,  the 
evaluation  of  the  terms  included  in  the  first  brackets. 

B.  Variance  in  N  established  by  equation  uncertainties 

The  standard  conditions  required  for  assignment  of  values  to  the  constants 
a  through  d  in  Equation  133  may  be  taken  as  those  defined  by  the  U.  S.  Standard 
Atmosi^ere.  This  atmosphere  has  an  equivalent  thickness*,  L,  in  the  zenith  di¬ 
rection  of  8.4  km,  an  oxygen  partial  density  of  2.50  x  102  gm/m“^,  and  a  tem¬ 
perature  of  288®K.  Also,  consider  a  standard  atmospheric  water-vapor  content 
of  2. 0  pi"tipitable  centimeters.  The  relationship  between  precipitable  cenitmeters, 
w  and^^  is 

W  4  5i 

p  -  '  lO  '  (134) 

^Equivalent  thickness  is  defined  by  the  hei^t  to  which  a  homogeneous  atmosphere 
would  extend  at  temperature  T  and  pressure  p. 
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when  W  is  given  in  precipttable  centimeters,  and  the  path  length  L  is  given  in 
meters.  In  this  homogeneous  atmosphere,  the  integrated  partial  densities  ma." 
be  considered  simply  as  the  product  of  density  by  path  length.  Taking  the  k-  and 
W02 -uncertainties  as  given  in  Equations  74  and  77,  the  standard  conditions  listed 
above,  and  the  atmospheric  path  as  the  equivalent  atmospheric  path,  the  numerical 
values  of  the  factors  appearing  in  the  first  set  of  brackets  on  the  right-hand  side 
of  Equation  133  are  as  given  in  Table  6.  Acc''nting  the  definition  of  air  mass  as 
the  ratio  of  the  path  length  tlirou^  the  atmosphere  at  any  given  angle  to  the  path 
length  tov.nrd  the  zenith  (Reference  46),  the  values  are  listed  as  applying  to  one 
standard  air  masa. 

Table  6:  Values  of  the  Constants  a,  b,  c,  and  d,  and  of  the  Error 
ComponentaJc]^UU<2'^^>  ^  N  x'ms  Error 

Equation,  Through  One  Standard  Air  Mass 


Constant 

Numerical  Values 

Constant  Error  Component 

Error  Component 

a 

1.02  x  lo"^ 

±  3.7  X 10"® 

Ak]^ 

b 

2.0  X  10^ 

±  3. 92  X  10"2 

Akg 

c 

6.95  X 10^ 

:fcl4.3 

Akg 

d 

-  9.75  X 10® 

±  5.0  X 10"® 

The  data  in  Table  6  permits  the  estimates,  shown  in  Table  7,  of  the  contri¬ 
bution  of  each  of  the  error  components  to  the  total  variance  in  N. 

Table  7:  Contributions  ot  the  Error  Components  to  the  Total  Variance  in  N 

^  Total  N 

^Variance  Total  N  rms 
Per  Air  Error  Per 
Mass  Air  Mass 

1.43  xlO^  6.15  xlO^  9.88  xlO^  4.88  xlO^  7.59x10^  8.71  xlO^ 

ihe  data  contained  in  Table  7  define  the  inherent  limitation  to  the  accuracy 
attainable  in  line-integrated  refractivity  measurement  through  application  of 


Contributions  to  Total  Variance  due  to: 
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oxygen  and  water-vapor  optical  depth  meaaurement  data  to  Equation  79.  That  is, 
assuming  errorless  determinations  and  <I^/t)  could  be  made,  the 

Inferred  N  data  would  be  subject  to  an  rms  error  in  N  of  about  8, 71  x  10^  meters 
per  air  mass,  or  to  about  0,04  percent  in  N.  To  such  precision  then,  the  equa¬ 
tion  for  line- Integrated  refractlvity  becomes,  from  Equation  79  with  the  k-  and 
WQg-values  to  three  significant  figures, 

'lS=0.‘?44^  -h  (135) 

'U>  \  r  /  * 

C.  Variance  in  N  established  by  radiometric  measurement  errors 

The  terms  in  the  second  set  of  brackets  on  the  rig^it-hand  side  of  Equa¬ 
tion  133  determine  the  contribution  of  the  uncertainty  in  N -specification  due  to 
errors  (the-A/ys)  in  measurement  of  ^>2*  and  Two  factors  may  be 

considered  as  contributing  to  such  measurement  errors;  viz. ,  instrumental 
"noise, "  and  atmosi^ere-induced  errors.  Consider  first  the  latter  factor,  as¬ 
suming,  in  effect,  "perfect"  instrumentation.  This  will  establish  the  atmospheric 
"noise  level"  to  which  the  instrumentation  may  be  designed. 

(1)  Modified  transmission  equations 

The  basic  premise  of  the  radiometric  approach  to  the  measurement 
of  the  integrated  partial  deitotty  data  demanded  by  Equation  135  is  that  the  attenu¬ 
ation  (or  emission)  of  energy  in  an  atmospheric-constituent  absorption  band  pro¬ 
ceeds  in  portion  to  the  total  amount  of  the  constituent  in  the  psdh  over  which  the 
energy  is  propagated.  To  Isolate  energy  in  the  absorption  bands  of  interest, 
narrow  band-pass  interference  filtera  may  be  used  in  Infrared  regions  of  the 
spectrum,  and,  at  radio  wavelengths,  suitable  band-pass  amplifiers  may  be  used 
to  accomplish  the  same  result.  But  where  such  mechanisms  may  be  extremely 
efficient,  in  general  the  passed  bands  will  be  of  a  width  tO  require  that  they  be 
considered  as  delineating  still  an  unresolved  spectral  abso^ion  line  structure. 

A  description  of  transmittance  in  such  a  band  requires  a  more  complicated  trans¬ 
mission  function  (i.e.,  explicit  definition  of  Equation  80)  than  the  simple  Bouguer- 
Beer  Law  relationship  of  Equation  84  (Reference  47).  That  is,  in  the  development 
of  a  function  to  describe  atmosF^eric  transmittance  over  a  group  of  unresolved 
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lines,  it  must  be  recognized  that  the  net  effect  represents  the  summing  up,  over 
the  frequency  interval  defined  by  the  band-width  of  the  filter,  of  the  individual 
effects  of  a  multitude  of  lines,  each  contributing  to  absorption  according  to  Equations 
66  &  67. The  transmitted  spectral  radiation  is  therefore  not  uniform,  since  in 
the  region  of  absorption  lines  it  is  decreased,  and  at  the  centers  of  intense  lines 
it  is  possible  for  the  absorption  to  be  complete.  Thus,  as  the  thickness  of  the 
transmitting  gas  is  changed,  there  is  a  change  not  only  in  the  intensity  of  the  radi¬ 
ation,  but  also  in  its  spectral  composition,  which  predetermines  the  peculiarity 
of  the  absorption  laws  in  spectral  measurements  of  unresolved  structures. 

In  general,  an  isolated  group  of  lines  in  an  unresolved  band  is  characterized 
by  differing  mutual  distributions,  degrees  of  overlapping,  integral  intensities, 
half-widths,  and  values  of  effective  collision  diameters.  For  such  a  system  a 
general  solution  to  the  problem  of  calculating  the  band  transmittance  is  not  pos¬ 
sible,  since  it  is  peculiarly  complex,  and  awkward,  and  requires  a  knowledge  of 
the  individual  ch'«i-acter  of  each  line.  Accordingly,  vaiious  model  systems  have 
been  proposed,  approaching  in  greater  or  lesser  degree  to  actual  molecular 
spectra.  Elsasser,  for  example,  proposed  a  system  of  lines  of  equal  intensity 
and  width  occurring  with  equal  distances  between  lines,  an  adequate  representa¬ 
tion  of  the  oxygen  spectrum  (Reference  48).  In  his  statistical  model.  Goody  intro¬ 
duced  a  system  of  randomly  distributed  lines,  the  intensities  of  which  are  similarly 
randomly  distributed,  as  a  representation  of  the  water-vapor  spectnun  (Reference 
23).  From  these  models,  the  following  transmission  functions,  defined  for 
Lorentizian  line  shapes,  may  be  stated.  First,  from  the  Goody  model,  the  trans¬ 


mission  in  a  water-vapor  band  is 


_ y 

/I  ^  V  • 


and  from  the  Elasser  model  the  oxygen  transmission  is 
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where  i  denotes  the  average  line  intensity  in  the  spectral  interval,  d  the  mean 
spacing  between  lines. 

(2)  Interpolation  equations 

Ekiuations  136  and  137  can  be  approximated  (References  49  and  47) 
in  the  following  fashion.  For  large  water-vapor  optical  depths  ("strong"  absorption) , 
the  term  t/^/TrAV  in  the  demoninator  of  the  argument  of  the  exponential  in  Equa¬ 
tion  136  is  much  greater  than  unity,  and  the  equation  reduces  the 

^  Zu  ^  (138) 

where  Cj^  denotes  a  generalized  absorption  coefficient  incorporating  the  line  in¬ 
tensity  and  line  spacing  parameters.  Conversely,  at  small  optical  depths 
(.^/ttAz/^1,  and 


4  ^ 


3  ruLf 


In  the  intermediate  case  when  1, 

^  (140) 

In  the  case  of  oxygen,  "strong"  band  absorption  is  associa.'.d  with  pronounced 
"overlapping"  of  the  regularly  spaced  lines.  Thus  the  strong  absorption  is  de¬ 
fined  when  the  term  2ir(ltVd  in  the  arguments  of  the  transcender' al  functions  on 
the  right-hand  side  of  Equation  137  are  much  laiger  than  unity.  When  such  is 


the  case, 


^  TV] 


fo„ 


In  the  weak  absorption  case,  2  <tAi^S  «1,  and 

For  the  intermediate  case 


7^  ^  I 
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Keep  in  mind  that  Equations  138,  139,140,  141,  142,  and  143  are 
basically  interpolation  equations  that  express  their  respective  general  transmis¬ 
sion  laws  in  terms  of  three  successive  approximations.  At  just  what  value  of 
transmittance  a  given  form  of  one  of  the  interpolations  takes  over  from  another 
cannot  be  stated  absolutely.  There  obviously  must  be  a  smooth  transition  from 
one  to  the  other. 

(3)  Line -broadening  errors 

A  peculiarity  to  be  noted  in  Equations  138  through  143  is  that  in  the 
"square-root  law"  regimes  of  both  water-vayor  and  oxygen  transmission,  T  is 
a  function  not  only  of  the  integrated  densities,  but  also  of  the  line  width, ^  y. 

It  is  likely,  moreover,  that  in  performing  the  transmission  measurements  to  ob¬ 
tain  the  necessarj'  data  for  computing  line -Integrated  refract!  vity  in  ranging  oper¬ 
ations,  a  good  part  (if  not  the  bulk)  of  such  measurements  will  be  performed  at 
transmittance  values  where  the  square-root  laws  are  applicable.  Thus,  in  gener¬ 
al,  the  inversion  of  the  T-data  to  yield  ^-data  will  involve  relations  of  the  type 


In  the  manner  of  Equation  133  then,  the  variance  of  the  uncertainty  In^  -specification 
can  be  estimated  from  an  expression  of  the  form 


Two  effects  are  thus  seen  to  contribute  to  errors  in  />data;  those  arising  from 
T-mensurement  errors,  and  those  stemming  from  uncertainties  in  the  specifica¬ 
tion  of  the  line -width  parameter,  ^  V  • 

Consider  first  that  the  latter  source  of  error  in  ^-specification 

is  detersiined  by  the  pressures  of  both  the  absorbing  and  foreign  gases  of 

the  atBosphere,  by  the  teaperature,  and  by  the  effective  collision  diasH 

eters  of  the  absorbing  and  foreign  gas  aolecules.  These  collisions  di- 

aaeters  are  deterained,  too,  relative  to  the  pressures  of  the  gases  and  not 

by  the  nuaber  of  collisions  (Reference  ii7).  The  net  result  is  that  in  teras 

of  pressure  and  teaperature,  4yean  be  expressed  as 
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p.T) 


where  the  subscripts  denote  the  values  of  the  various  parameters  under  some  set 
of  standard  conditions. 

Equation  146  indicates  that,  due  to  the  dependence  of  ^on  p  and  T, 
the  line-width  parameter  also  is  a  function  of  position  in  the  atmosphere.  There¬ 
fore,  the  integration  of  partial  density  called  for  to  effect  an  evaluation  of  N  must 
really  be  an  integration  of  the  product  of  partial  density  by  line  width;  i.  e.  , 


o 

^  ids 


The  integration  called  for  in  Equation  147b  can  be  accomplished 
explicitly  only  if  one  is  willing  to  assume:  first,  fjj  explicit  form  for  f  ;  second, 

it 

horizontal  homogeneity  and  an  exponential  variation  uf  both  pressure  and  temper¬ 
ature  with  altitude  so  that  an  explicit  expression  for  f  as  a  function  of  position  in 

2 

slant -path  operations  may  be  specified;  and  third,  that  the  mixing  ratio  of  remains 
constant  throughout  the  atmosphere.  When  this  is  done,  the  integration  called  for 
yields  an  expression  for  an  equivalent  optical  depth  (defined  as  a  function  of  slant- 
path  end  point),  surface  conditions  of  temperature  and  pressure,  and  slant -path 
elevation  angle  (see,  for  example,  Reference  50);  i.e. . 


On  the  basis  of  an  equation  of  the  type  depicted  by  Equation  148  . 
an  explicit  functional  relationship  could  be  derived,  in  the  manner  of  Equation  131. 
for  the  uncertainties  in  ^-specification  attributable  to  uncertainties  in  p  and  in  T, 
and  thus  Ay.  However,  all  assumptions  required  to  arrive  at  such  an  equation 
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are  questionable  at  best.  Regarding  the  case  of  water -vapor  absorption  espec¬ 
ially,  the  assumption  of  the  constant  mixing  ratio  is  particularly  bad.  On  a  sta¬ 
tistical  basis,  the  constant  mixing  ratio  might  at  least  approach  something  near 
actuality,  but  on  a  day-to-dt^r  operational  basis,  atmospheric  water  vtupor  is 
usually  found  to  occur  in  well-defined  layers  at  various  heights  throughout  the 
depth  of  the  atmosphere. 

Rather  than  go  the  route  of  the  explicit  error  analysis  then,  an 
ai^al  is  made  to  direct  observation  as  to  what  the  effects  of  pressure  and  tem¬ 
perature  line  broadening  effects  might  have  on  the  accuracy  of  p  measurement. 
Foster,  Volz,  and  Foskett  (Reference  51)  have  developed  a  photoelectric  spectral 
hygrometer  for  the  ESSA  for  determining  integrated  water-vapor  densities.  The 
Iqrgrcmeter  monitors  the  relative  radiant  intensity  of  direct  solar  energy  trans¬ 
mitted  in  a  water-vapor  absorption  band,  and,  in  the  manner  specified  previous¬ 
ly,  compares  the  absorbed  energy  to  the  energy  received  in  a  nearby  atmospheric 
window.  The  instrument  is  mounted  on  a  clock-driven  equatorial  telescope  mount, 
so  that  the  device  tracks  the  Sun  from  sunrise  to  sunset.  The  device  was  calibra¬ 
ted  through  comparisons  of  the  hygrometer  outputs  to  measures  of  precipitable 
water  vapor  calculated  from  standard  radiosonde  data.  In  this  way.  pressure  and 
temperature  line -broadening  effects  were  "calibrated  in"  to  the  system.  Having 
calibrated  their  system,  subsequent  operation  has  shown  that  changes  in  total  pre¬ 
cipitable  water,  as  detected  by  the  spectral  hygrometer,  were  r^nsistently  compar¬ 
able  to  those  changes  Indicated  by  successive  radiosonde  uieasurements,  regardless 
of  the  altitude,  (hence,  the  temperature  and  pressure)  at  which  the  changes  occurred. 
Subsequent  attempts  were  made  in  an  effort  to  obtain  a  quantitative  appraisal  of  the 
line -broadening  effect  on  /^-determinations,  with  the  results  that  the  magnitude  of 
such  effects  were  found  to  result  in  errors  not  normally  in  excess  of  2  or  3  percent 
of  the  total  amount  of  water  vapor,  with  the  ext4.eme  error,  under  conditions  of 
extreme  water-vapor  layering,  being  on  the  order  of  10  percent.  As  indicated 
earlier,  these  estimates  may  also  be  considered  representative  of  the  accuracy 
to  which  the  measurement  of  the  integrated  ratio  of  water  vapor  density -to -temper¬ 
ature  may  be  performed  at  microwave  frequencies.  These  estimates,  based  on 
the  findings  of  Foster,  et  al,  are  also  in  agreement  with  those  of  Kondrat'ev,  et  al 
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(Reference  52).  Because  the  oxygen -mixing  ratio  is  constant  throughout  the 
depth  of  the  atmosphere,  the  layering  problem  is  largely  eliminated  as  a  source  of 


-specification  uncertainty.  Horizontal  heterogeneity  still  will  give  an  eqiava- 
layering  effect;  'nierefore,  although  the  oxygen  data  will  be  much  more  precise 


than  that  of  water  vapor,  an  error  of  from  ±  0. 5  to  1  percent  must  still  be  accept¬ 


ed.  The  latter  estimate  is  substantiated  by  the  work  of  Westwater  (Reference  27). 


(4)  N-errors  related  to  transmittance  measurement  error 


Proceeding  now  to  the  contribution  of  T -measurement  errors  to  uncer¬ 


tainties  in  ^-specification  (Equation  145),  consider  the  implications  of  the  various 
forms  of  the  water-vapor  and  oxygen  transmission  laws  (as  stated  in  Equations  138 
through  140,  and  141  through  143,  respectively)  as  to  the  relation  between  Ar  and 
A^.  Defining  the  relative  errors  in  ^-specification  and  T-measurement  as5  (/>)  * 
£^/p  and  |  (r)  =  At/t,  respectively,  the  following  relations  are  obtained  by 
differentiation  of  Elquatlons  138  and  139,  and  141  and  142; 
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Equations  lii9  through  152  are  presented  in  Figure  19.  The  solid 
line  depicts  the  ratio  of  the  relative  error  in  the  integrated  water -vapor  density 
specification  to  the^relatlve. error  In  water-vapor  transmittance  measurement. 
The  dashed  line  depicts  the  same  ratio  for  the  case  of  oxygen.  The  seg¬ 
ments  in  the  center  of  the  curves  denote  the  "transition"  regions,  as  implied  by 
Equation)  140  and  143.  between  the  "strong"  and  "weak"  absorptions  of  Equa¬ 
tions  138  and  141,  and  140  and  143,  respectively.  The  curves  indicate  that  over 
a  wide  range  of  transmittances,  the  relative  error  in  ^-specification,  as  related 
to  the  relative  error  in  r -measurement,  is  not  very  sensitix'e  to  the  actual  mag¬ 
nitude  of  the  measured  transmittance.  A  rule  of  thumb  may.  therefore,  be  es- 
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tablished  for  the  purpose  of  the  error  analysis  to  follow;  for  transmittances  rang¬ 
ing  from  about  10  to  70  percent,  the  relative  error  in  J)  will  be  about  twice  that 
in  T.  Also,  because  of  the  similar  nature  of  the  line -broadening  effects  at  micro- 
waves,  this  rule  of  thumb  may  be  expected  to  apply  to  the  case  of  microwave  as 
well  as  infrared  absorption. 

The  T-error  increment  now,  in  term.s  of  the  actual  measurements 
performed  by  a  differential  absorption  spectrometer  system,  is  obtained  by 
differentiation  of  Equation  94: 
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But  also  from  Equation  (94) ,  notice  that 


and 
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Using  Equations  154  and  155  appropriately  in  Equation (153) ,  the  latter  trans¬ 
forms  to 
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The  terms  in  U.o  brackets  of  Equation  (156)  define  relative  changes  in  r<  eived 
sensing  and  reference -band  energies  attributable  to  effects  other  than  energy'  ab¬ 
sorption  by  the  atmospheric  constituent  being  sensed.  If  the  above  percentage 
chaiigcs  are  equal,  the  right-hand  side  of  the  equation  is  zero,  and  no  e.  or  in 
absorption  transmittance  measurement  will  result.  This,  of  course,  is  the  reason 
for  employing  a  reference  band  in  the  first  place.  However,  several  processes 
that  are  wave  length -dependent  could  act  to  attenuate  both  the  reference  and  sensing- 
band  energy  intensities,  th’is  making  the  right -band  side  of  Equation  (156)  nonzero, 
with  the  resultant  degradation  of  ^-measurement.  Such  processes  might  include 
(1)  differential  attenuation  caused  by  the  atmospheric  continuum  spectnun,  and  (2) 
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scintillation  effects.  These  two  effects  can  again  be  considered  to  act  in  an  un- 
corrclatcd  manner  in  the  measurement  degradation. 

(A)  Effects  of  the  continuum  spectrum  — Infrared  continuum 

The  infrcred  continuvim  spectrum  is  established  by  three  pro¬ 
cesses:  (1)  scattering  by  the  atmospheric  aerosol,  (2)  absorption  by  the  particles 
of  which  the  aerosol  is  composed,  and  (3)  absorption  due  to  wings  far  removed 
from  the  resonant  frequency  line  absorptions  of  all  the  various  absorbing  atmos¬ 
pheric  constituents.  Because,  however,  these  processes  combine  to  yield  the 
continuum  effect,  the  continuum  transmission  can  be  described  by  the  Bouguer 
law, 


At  7^  =  e'w\  ^  (157) 

where  represents  the  continuum  transmittance,  c'  the  continuum  extinction  co¬ 
efficient.  and  m  the  air  mass.  With  regard  to  c',  it  is  significant  to  note  that 
v  ithin  the  limits  of  the  Bouguer  law,  the  ratios  in  the  brackets  on  the  right-hand 
side  of  Equation  156  constitute  the  exact  definitions  of  the  extinction  coefficient 
(Reference  53).  Accordingly,  Equation  156  may  be  expressed  as 


That  is,  the  relative  error  ’.a  the  absorption  transmittance  ascribable  to  the  con¬ 
tinuum  extinction  effect  is  given  simply  by  the  difference  between  the  sensing  and 
reference  band  continuum  extinction  coefficients. 

Gates  (Reference  .54)  presents  the  relationship  between  the  con¬ 
tinuum  extinction  coefficients  and  wavel*;ngth  shown  in  Figure  20.  Also  Indicated 
are  the  ordinates  locating  the  1.2'^  oxygen  band,  the  1.32  water -vapor  band, 
and  the  1. 22^14.  window.  With  the  window  serving  as  the  reference  band  for  the  other 
two  absorption  bandsy^the  differences  between  the  absorption  and  reference -band 
continuum  e.xtinction  coefficients  are  estimated  as  being  on  the  order  of  10  for 
both  the  water-vapor  and  oxygen  bands.  Estimating  the  /^-specification  errors 
as  twice  these  T-measurement  errors,  it  can  still  be  said  that,  with  the  reference- 
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scintiilatiop  effects.  These  two  effects  can  again  be  considered  to  act  in  an  un- 
ccrrelated  manner  in  the  measurement  degradation. 

{A)  Effects  of  the  coRtiriuma  spectrum  — Infrared  continuum 

The  infrared  continuum  spectrum  is  established  by  three  pro¬ 
cesses:  (1)  scattering  by  the  atmospheric  aerosol,  (^)  absorption  by  the  particles 
of  which  the  aerosol  is  composed,  and  (3)  absorption  due  to  wings  far  removed 
from  the  resonant  frequency  line  absorptions  of  all  the  various  absorbing  atmos¬ 
pheric  constituents.  Because,  how'ever,  these  processes  combine  to  yield  the 
continuum  effect,  the  continuum  transmission  can  be  described  by  the  Bouguer 
law. 

7^  =  c'^y\  ^  (157) 

where  represents  the  continuum  transmittance,  c'  the  continuum  extinction  co¬ 
efficient.  and  m  the  air  mass.  With  regard  to  c',  it  is  significant  to  note  that 
w  ithin  the  limits  of  the  Bouguer  law,  the  ratios  in  the  brackets  on  the  right-hand 
side  of  Equation  156  constitute  the  exact  definitions  of  the  extinction  coefficient 
(Reference  53).  Accordingly,  Equation  156  may  be  expressed  as 
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That  is,  the  relative  error  In  the  absorptions  transmittance  ascribable  to  the  con¬ 
tinuum  extinction  effect  is  given  simply  by  the  difference  between  the  sensing  and 
reference  band  continuum  extinction  coefficients. 


Gates  (Reference  54)  presents  the  relationship  between  the  con¬ 
tinuum  extinction  coefficients  arid  wavei»'ngth  shown  in  Figure  20.  Also  indicated 
are  ihe  ordinates  locating  the  1. 2“^  oxygen  band,  the  1..32  /«  water -vapor  band, 
and  the  1.22^  window.  With  the  window  serving  as  the  reference  band  for  the  other 
two  absorption  bands^f  the  differences  between  *he  absorption  and  reference -band 
continuum  e.xtinction  coefficients  are  estimated  as  being  on  the  order  of  10  for 
both  the  water-vapor  and  oxj'gen  bands.  Estimating  the  /^-apecification  errors 
as  twice  these  T-measurement  errors,  it  can  still  be  said  that,  with  the  reference- 
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THE  INFRARED  CONTINUUM  EXTINCTION  COEFFICIENT  AS  A  FUNCTION 
OF  WAVELENGTH 


to -absorption -ratio  mode  of  operation,  the  continuum  extinction  effects  (as  com¬ 
pared  to  line -broadening  effects)  on  the  degradation  of  the  integrated  vater-vapor 
and  oxygen  density  data  are  negligible. 

(B)  Effects  of  the  continuum  spectrum — Microwave  continuum. 

At  microwave  frequencies,  the  continuum  spectrum  of  greatest 
concern  is  that  presented  by  the  occurrence  of  clouds.  Deirmendjian  has  calcu¬ 
lated  cloud  extinction  coefficients  as  a  function  of  microwave  wavelength,  assum¬ 
ing  representative  cloud  models  as  regards  droplet -size  distributions  and  concen¬ 
trations  (Reference  55).  The  results  are  shown  in  Figure  21.  Also  indicated  are 
the  ordinates  denoting  the  water -v8^)or  absorption  band  and  its  reference,  and 
the  oxygen  band  and  its  reference. 


Contrary  to  Figure  20,  the  data  contained  in  Figure  21  Indicate 
a  pronounced  dependence  of  extinction  coefficient  on  wavelength.  However,  from 
the  figure,  the  difference  in  extinction  coefficients  at  the  wavelengths  of  the  water- 
vapor-absoiiition  and  reference  bands  is  calculated  as 


0.003  -0.00/  =  0  Ooz 


r. 


and  that  between  the  oxygen -absorption  and  reference  bands  as 


(159) 


At 
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0.0 za  -  o  ooZj  -  0.0/ S  (160) 

Even  multiplying  the  right-hand  side  cf  Equation  159  by  two,  in  accordance  with 
the  dictates  of  Figure  19,  the  relative  error  in  attributable  to  the  cloud 

continuum  extinction  is  still  an  order  of  magnitude  less  than  the  previously  speci¬ 
fied  error  associated  with  line-broadening  effects.  Therefore,  the  error  may  be 
m  glee  ted. 
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Figure  21:  COMPUTED  EXTINCTION  CURVES  FOR  CLOUD  AND  PRECIPITATION  AS  A 
FUNCTION  OF  THE  WAVELENGTH  (After  DeirtTiend|ian,  1963) 


In  the  case  of  oxygen,  however,  twice  the  extinction  coefficient  difference 
implies  a  relative  error  in  /?^  measurement  of  better  than  1  percent.  The 
oxygen  continuum  extinction  error  therefore  should  not  be  dismissed,  but  rather 
attempts  should  be  made  to  actually  correct  the  data.  Suppose,  as  suggested  by 
Deirmendjian,  the  cloud  continuum  extinction  coefficient  is  expressed  as  (Refer¬ 
ence  55) 

(i( h)  -  A  X  ^  (161) 

where  A  is  a  parameter  depending  on  the  physical  composition  (i.e. ,  cloud-drop¬ 
let-size  distribution,  concentration,  etc.)  of  the  actual  clouds  extant,  and  q  is 
a  constant.  Goldstein  (Reference  56)  suggests  the  value  q  =  2.  Suppose  further 
that  the  ratio  between  the  two  microwave  reference  signals  is  acquired  along 
with  the  sensing- to- reference  signals.  Then,  in  line  with  Equation  155 
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where  Kj.g£  denotes  the  ratio  of  the  intrinsic  source  energy  at  the  first  reference 
wavelength  to  that  at  the  second  reference  wavelength,  and  1^.  and  denote  the 
corresponding  continuum  transmittances.  From  the  Bouguer  extinction  law  it  also 
follows  that 
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The  right-hand  side  of  Equation  163,  may,  by  virtue  of  Equation  161,  be  redefined 
to  yield 
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Thus,  measurement  of  the  ratio  of  the  reference- oignal  intensities  suffices  to 
evaluate  the  parameter  A  required  for  estimating  the  continuum  extinction  co¬ 
efficient,  in  accordance  with  Equation  161,  for  any  wavelength.  The  right-hand 
side  of  Equation  158  is  thus  quantified,  and  the  numerical  value  of  the  percentage 
error  in  T'  measurement  is  established.  The  data,  as  defined, by  Equatioa 
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155,  may  therefore  be  corrected  appropriately,  probably  to  within  10  percent  of 
the  value  given  in  Equation  160.  The  residual  error  in  would  thus  be  some- 
thing  more  like 

ZS  f' 

^  <3*00/8 

r. 

One  further  advantage  may  be  gained  from  measuring  the  ratio  of  the  microwave 
reference  signals.  The  percentage  error  in  oxygen- absorption  measurement 
therein  established  defines  the  magnitude  of  the  cloud  continuum  transmittance. 
In  turn,  this  transmittance  is  a  function  of  (and  therefore  provides  a  measure 
of)  the  li<iuld  water  content  of  the  cloud  (Reference  56).  Knowledge  of  the  cloud's 
li(]uid  water  content  will  permit  estimation  of  the  contribution  to  line-integrated 
refractlvlty  established  by  the  cloud  droplets  themselves. 

(C)  Effects  of  mean-temperature  specification  errors 

When  microwave  oxygen  measurements  are  performed,  the 
acquired  data,  to  yield  the  correct  form  of  the  oxygen  term  in  the  N-equation, 
must  be  operated  on  in  accordance  with  Equation  112.  That  is,  the  integrated 
oxygen- absorption  coefficient  must  be  multiplied  by  the  mean  temperature  to 
obtain  the  desired  integrated  oxygen  density: 


4  ^  m 


An  additional  source  of  error  in  specification  thus  becomes  evident;  consider 


An  error  in  ^T)  measurement,  therefore,  also  contributes  to  an  error  in 

Specification.  Specifically  the  <^T^-  caused  relative  error  in  is  pre¬ 
cisely  that  of  the  T^  error  Itself;  i.  e. , 

In  effect.  Equation  169  constitutes  an  additional  term  that  must  be  added  to  the 
right-hand  side  of  Ecjuatlon  146  to  completely  specify  the  error  in  ^  measure¬ 
ment  under  cloudy  conditions.  The  magnitude  of  Equation  169  would  more  than 
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likely  contribute  an  additional  percentage  of  error  in  .  That  is,  the  mean- 
temperature  specification  error  would  result  in 

-0.01  (170) 

The  above  estimate  implies  a  /^(<T>)  error  of  order  ±  S^K. 

(D)  Scintillation  effects 

Intensify'  (and  phase  and  amplitude)  scintillations  of  signals 
received  from  remote  energy  sources  are  caused  by  the  scattering  of  direct 
waves  by  random  inhomogeneities  in  the  atmospheric  refractive  index  field. 

To  good  approximation,  the  inhomogeneittes  can  be  treated  as  equivalent  opti¬ 
cal  elements,  and  the  relative  signal  fluctuations  caused  them  explained  by 
the  Bom  single-scattering  solution  to  the  Schroedinger  wave  equation.  Using 
this  approach,  Wheelon,  for  the  case  of  a  plane  wave  falling  on  a  turbulent  medium 
in  which  is  immersed  an  upward-looking  receiver  (Figure  22),  presents 
the  following  expression  for  the  variance  of  the  relative  intensity  fluctuation 
(Reference  57) 


c 


(171) 


where  L  is  the  atmospheric  path  through  turbulence  (Figure  22),  )\  is  the  wave¬ 
length  of  the  propagated  energy;  where  k  now  denotes  the  wave  number  of  the 
turbulent  refractive  index  eddies,  S(k)  denotes  the  spectrum  of  that  turbulence. 

Note  once  more  that  the  left-hand  side  of  Equation  171  defines  precisely  the 
» 

terms  on  the  rig^t-hand  side  of  Equation  156,  which  are  now  required  to  esti¬ 
mate  scintillation  effects  on  7^  measurement. 

The  function 

attains  its  physical  significance  through  recognition  that  kM  is  the  scattering 
parameter-^- A (where  $  with  1  the  linear  dimension  of  the  turbulent 

eddies)  appropriate  to  each  eddy  size  in  the  hierarchy  described  by  S(k).  The 
eddies  that  are  large  compared  with  ^  give  Fresnel  scattering,  whereas 
large  values  of  7^  correspond  to  Fraunhofer  scattering.  Equation  172  in 
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22:  SCATTERING  BY  REFRACTIVE  INDEX  INHOMOGENEITIES  AS  THE 
CAUSE  OF  SCINTILLATION  IN  THE  IDEALIZED  LINE-OF-SIGHT 
PROBLEM  OF  A  SOURCE  VIEWED  UPWARD  THROUGH  A 
TURBULENT  MEDIUM 
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combination  with  171  establishes  the  following  limiting  expressions  for  the  variance 
in  intensity  fluctuations: 

z  I 


The  Fresnel  approximation  of  Equation  173  will  always  be  valid  for  A  of  infrared 
wavelengths,  whereas  the  Fraunhofer  approximation  of  Equation  174  must  be 
considered  if  /\  is  of  microwave  wavelengths.  Thus,  the  variance  of  the  relative 
intensity  fluctuations  at  infrared  wavelengths  is  <to  good  approximation)  given  by 
the  cube  of  the  path  through  turbulence,  times  the  integral,  over  all  wavenumber 
magnitudes  of  the  third  moment  of  the  refractive  index  spectrum.  No  explicit 
dependence  on  wavelength  is  Indicated.  At  microwave  wavelengths,  however, 
such  a  wavelength  dependence  is  indicated,  with  the  variance  in  relative  intensity 
varying  as  the  reciprocal  of  wavelength  squared  times  turbulence  path  length  to 
the  first  power  times  the  integral  of  the  first  moment  of  the  refractive  index 
spectrum.  Each  case  will  be  considered  separately. 

Scintillation  at  infrared  wavelengths 

For  analytical  purposes,  a  Gaussian  refractive  Index  correlation  function 

_  ^ 

may  be  assumed.  The  associated  spectrum  is 

5  Ck)  (176)  > 

where  is  the  variance  of  the  refractive  index  lluctuations4<n  •  The  result¬ 

ing  elementary  integration  called  for  by  Equation  173  may  now  be  performed  to 


yield 


'  r" 


where  1^  denotes  the  turbulence  length  scale.  The  relative  intensity  variance  is 
thus  related  to  the  variance  in  the  refractive  index.  The  standard  deviation  of 


relative  intensity  is  thus 
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and  if  the  refractive  index  fluctuations  are  assu'red  to  be  normally  distributed, 
an  average  intensity  fluctuation  can  be  calculated  as 

where  D(n),  the  mean  deviation,  is  taken  as  4/5  the  standard  deviation.  Accord¬ 
ing  to  Equation  156  then,  scintillation-caused  transmittance  errors  are  dependent 
on  the  difference  between  sensing  and  reference- band  mean  deviations  in  refrac¬ 
tive  index.  Specifically, 

^  [  0(^1  ~  OCnl  ]  (ISO) 

As  compared  to  temperature  and  pressure  influences,  water  vapor  only  slightly 
affects  refractive  index  at  infrared  wavelengths.  Edlen  (Reference  58)  expresses 
refractive  index  for  any  atmospheric  state  in  terms  of  n  in  standard  air  as 

Incremental  changes  in  n  due  to  turbulence  in  the  pressure  and  temporature 
fields  would,  therefore,  be  estimated  from 


The  lunction  f(p,  T)  is  not  wavelength-dependent,  so  the  ratio  of  n-fluctuations 
at  different  wavelengths  would  simply  ue  as  the  ratio  of  standard-air  refractive 
indices.  Mean  deviations  would  show  the  same  relationship  so  that 


and,  therefore. 


-  D(nl  _  0(n),.,  / -  1  ) 


Considering  the  sensing  and  reference  infrared  wavelengths  "as  specified" 
Edlen's  normal  dispersion  formula  indicates  that  the  refractivity  ratio  in  the 
parenthesis  on  the  right-hand  aide  of  Equation  184  will  differ  from  unity  in  about 
the  fourth  decimal  place.  Therefore, 
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At  visible  wavelengths,  Reiger  reports  a  turbulence-generated  standard  deviation 
in  n  of  about  3  x  10“®  (Reference  59).  This  estimate  is  based  partly  on  measure¬ 
ment  and  partly  by  inference  from  scintillation  data.  The  mean  deviation  of  re¬ 
fractive  index  will  also  be  on  the  order  of  10"®,  and  the  final  order- of- magnitude 
estimate  for  the  diR'erence  in  infrared-sensing  and  reference- band  mean  deviations 


is  thus 


i-  k 

- 


0  /o 


In  view  of  the  order  of  magnitude  defined  above,  scintillation  effects,  through 
suitable  instrument  design  (i.  e. ,  simultaneous  sampling  of  all  bands  through  a 
single  aperture)  can  be  reduced  to  negligible  proportions  in  spite  of  the  dependence 
of  intensity  scintillations  on  the  3/2-  power  of  the  atmospheric  turbulence  path. 

For  if  turbulence  is  considered  to  be  confined  mainly  to  altitudes  below  10km, 
then  even  at  a  slant-path  elevation  of  1  degree,  L  will  only  be  on  the  order  of 

9 

10  km.  The  turbulence  length  scale,  1^,  is  on  the  order  of  10m  (Reference  60). 
The  coefficient  of  the  mean  deviation  difference  in  Equation  180  is  thus  on  the 
order  of  10^.  Thus 

0  ,c->' 

The  implication  contained  in  the  order-of- magnitude  Equation  187  is  borne  out 
by  the  work  of  Yates  (Reference  6),  wherein  be  could  detect  no  scintillation 
dependence  on  wavelength  in  the  infrared  portion  of  the  spectrum.  His  work, 
moreover,  involved  using  one  wavelength  near  the  N2O  absorption  band,  where 
anonialous  dispersion  effects  might  be  expjected  to  augment  the  mean  deviation 
differences  stated  in  F  ’>ation  186.  However,  no  such  effeci.  was  noted. 

2.  Scintillation  at  microwave  wavelengths 

With  the  same  Gaussian  correlation  function  assumed  for 
turbulence  in  the  microwave  refractive  index,  the  integration  called  for  in 


Equation  174  is  pjerformed  to  y<''ld 
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'  'yK-u.t-/e  />  ' 

Again  with  the  mean  deviation  in  n  taken  as  4/5  the  standard  deviation, 


(Tj 


Therefore,  from  Equation  156 

Mt  A  _ 


l)  (190) 

^  /«,  V  ^  y  l_  ^'fcrti  ^  r*/'  J 

Arguing  as  before  in  developing  Equation  183,  the  mean  deviation  in  the  sensing- 


band  refractive  index  can  be  expressed  in  terms  of  D(n)^^^  times  the  ratio  of 
refractivity  at  the  two  wavelengths.  Therefore, 

fqty\  _  (x^n’/’-ALj Qc„)  '  ]  ,19, 


yu.-wm.tfc. 


There  is  no  normal  dispersion  at  microwave  frequencies,  so  the  only  factor 
that  makes  the  refractivity  ratio  in  Equation  191  different  from  unity  is  the 
effect  of  anonudous  dispersion.  Sullivan  (Reference  62)  pomts  out  that  such 
anomalous  dispersion  effects  result  in  differences  in  refractivity  across  an 
absorption  band  on  the  order  of  10"^.  The  ratio  of  refractivities  at  frequencies 
within  and  outside  such  bands  will  be  on  the  order  of  10“'*.  Equation  191  may 


thus  be  approximated  as 

m) 


rr> 


A  1 

I 


A  .  A 

rr,  5«-.»5 J 


With  the  sensing  and  reference  wavelengths  Indicated  in  Figure  21,  the  term  in 
brackets  on  the  rirbt-hand  side  of  Equation  192  is  on  the  order  cf  10^. 
Therefore,  again  accepting  the  value  10~^  is  the  mean  deviation  in  refractive 


index  at  the  reference  wavelength,  it  may  be  stated  that 

(“)  0 

'  4.  ^  / 

Postulating  the  values  for  1^  and  L  mentioned  in  conjunction  with  the  analysis 
of  the  Fresnel  scattering,  the  final  order- of- magnitude  estimate  is  obtained  as 


•A<’ 


shown  below: 


■)  0  io 
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Again,  the  if/T)--  specification  error  associated  with  f'  measurement  errors 
of  the  orders  shov/n  above  are  negligible  in  comparison  with  the  errors  arising 
from  line-broadening  effects.  The  fact  that  scintillation  has  even  less  impact 
on  T'  measurement  at  microwave  wavelengths  than  at  infrared  wavelengths  (in 
spite  oi  the  explicit  microwave  scintillation  wavelength  dependence)  is  ascribable 
to  tiie  iiiucn  weaker  effect  oi  turbulence  path  length  at  microwave  frequencies. 

(5)  Considerations  of  passive  operational  node 

The  above  error  analysis  is  derived  when  considering  an  active 
system.  Because  of  the  basic  connection  at  any  frequency  between  absorpt’on 
and  emission,  however,  the  results  can  be  applie  i  directly  to  the  case  of  passive 
systems,  except  for  the  scintillation  effect.  Scintillation  phenomena  do  not  affect 
sky  brightness  temperature  measurements,  and  consideration  of  them  may  be 
dismissed  entirely.  Line-broadening  effects  on  emission  (i.  e. ,  sky  brightness 
temperature)  measurements,  however,  are  the  same  as  on  transmittance 
measurements.  And  insofar  as  only  microwaves  are  utilized  in  the  passive 
mode,  the  continuum  effects  are  similarly  the  same.  This  follows  because  the 
microwave  continuum  is  that  preseated  by  clouds,  and  the  process  that  establishes 
it  is  absorption. 

(6)  Cumulative  N  error 

,  Assuming  the  line-broadening,  the  continuum  extinction,  the  scin- 
tillation;  and  the  ('T^  -specification  errors  in  specification  to  be  uncorrelated, 
then  the  total  rms  v^lue*  of  or  of  a(^)  can  be  considered  to  consist  of  the 
square  root  of  the  sum  of  the  variances  in^  orA(£)  generated  by  each  effect 
individually;  i.  e. , 


^  /j.  line  fTuum 


The  above  e'qaatiop  is  an  expansiojfof  Equal  .-on  145,  and  di-fines  the  individual 
contributing  effec.ts^t^errort  ir^./’.or^f’/''.).  Th..-  first  t//o  tgrms  on  the  right- 
hand  side  of  Equation  195  are  active  contributors  to  emrs,  regardless  of  the 
conflg>i ration  of  the  radiometric  system;  1.  e,  <  .activ^^l^Wid  iR/mfcrbvyave 


(195) 


system,  all- microwave  active  system,  or  the  passive  system.  The  third  term 
contributes  error  only  in  the  active  mode  of  operation,  and  the  fourth  term  only 
when  microwave  absorptions  are  measured  to  determine  values  oiJ>  from 
(f/T)  data.  The  last  term  thus  applies  to  the  cloudy- weather  active  system  or 
to  the  passive  system.  The  magnitudes  of  the  various  contributions  to  the  total 
variance  inf  or  (f/T)- specification  are  summarized  in  Table  8, 

The  variance  shown  in  Table  8,  estimated  according  lo  the  dictates 
of  Figure  19  from  the  previously  specified  relative  eri-ors  -measurement, 
are  listed  by  contributing  effects  applicable  to  the  different  radiometric  measure¬ 
ment  system  configurations.  The  total  variance  listed  at  the  bottom  of  each 
column  represents  the  column  sum  in  accordance  with  Equation  195.  Ihe  rms 
value  is  simply  the  sq-.iare  root  of  total  variance.  The  associated  rms  line-inte¬ 
grated  refractivity  errors  for  one  standard  air  mass  are  listed  in  Table  9.  The 
values  listed  were  computed  in  accordance  with  Equation  133,  after  defining  N 
for  one  standard  air  mass  in  accordance  with  the  standard  conditions  listed  in 
Paragraph  6.  B.  and  Equation  135.  ^ 

Note,  from  fable  8,  that  the  line -broadening  effect  is  the  biggest  con- 
tributor  to  N  errors  in  the  hybrid  system.  In  both  the  all -microwave  active  and- 
passive  systems,  the  biggest  contributor  is  the  effect  of  error  in  the  <(T^ 
-specification. 

D.  Residual  errors  Iq  slant  range  and  elevation  angle 

# 

(1)  Dasic  residual  error  equations 

In  Paragraph  2.  A.  ,dS  was  considered  the  atmospheric  contribution  to  slant  range. 
Considering  that  the  uncertainty  in  the  specification, of  constitutes  'the  residual 
slant  range  error,'  then  from  Equation  9,  by  dlffereptiation,  .  ^  . 

A  (AS)  (196) 

where tbs  error  of  estimate  of  the  atmr  spheric  range  contribution,  denotes 
‘  x-estdual  err<Mf. 

•  ». 

,  H' 
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Tablp  9:  N-speci£icaticn  Errors  (per  Unit 
v'  Standard  Air  Mass)  Ascribable  to  At¬ 
mospheric  Effects  on  Different  Radio- 
metric  3y8tem  Configurations 


<IMS  N- error 
pel  Air  Mass 

Absolute  Error 
(Meters) 


Radiometric  ^stem 

Active  Hybird  Active  Microwave 

.  (Clear- Weather)  (Cloudy  Weather) 


Passive 


1.32  X  10 


4 


3.  54  X  10 


4 


3.54x  10 


4 


Relative  Error 
(Percent) 


0.54 


1.47 


1.47 


In  a  similar  manner,  considered  the  atmosphoric  contribution 

to  elevation  angle.  The  residual  elevation- angle  error  in  interferometric 
tracking  operations  (Figure  3)  can  be  obtained  from  differentiation  of  Equation  18 
The  result  is,  in  finite  increment  form; 

A(A(-)=  (197) 

where  dieA{A^'s)  denote  residual  errors  in  the  interferometer  angles. 
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(2)  Error  estimates 
(A)  Slant  range 

The  magnitude  of  the  residual  slant -range  error  is,  by  Equa¬ 
tion  (196),  linearly  related  to  the  magnitude  of  the  error  in  line -integration  re- 
fractivity.  The  N  error,  furthermore,  is  linearly  related  to  the  air  mass  value. 
Air  mass  is  simply  the  length  of  the  atmospheric  path  traversed  by  a  ray  (eman¬ 
ating  from  a  source  entirely  outside  the  atmosphere)  in  reaching  the  Earth,  meas¬ 
ured  in-  terms  of  this  path  when  the  source  is  at  the  zenith.  The  air  mass  value  to 
be  assigned  the  atmospheric  path  between  a  radar  and  a  target  thus  varies  with 
(1)  target  elevation  angle,  and  (2)  target  altitude  so  long  as  the  target  is  within  the 
atmosphere.  The  air-mass  variation  with  elevation  angle  io  shown  in  Figure  23. 
From  the  data  contained  in  Figure  23,  the  air-mass  variation  with  target  altitude 
can  be  estimated  f^o^l 

b- 


(198) 


where  m(  6^.h)  is  the  air-mass  value  of  the  atmospheric  path  in  the  direction  6^ 

to  target  height  h, ,  m(6  )  is  the  air  mass  for  the  total  atmospheric  path  air  mass 
1  0 

in  the  direction^Q,  f»(hi)  is  the  atmospheric  pressure  at  height  hj  andp^  is  the  sur¬ 
face  pressure.  The  value  ni(€j^)  is  the  total  atmospheric  path  air  mass  in  the  direc- 
tione  j,  where  the  local  line-of -sight  elevation  at  the  target,  is  given  by 


€ ,  =  90  ~  COS 


[ 


■1 


(199) 


The  geometry  of  the  problem  is  as  shown  in  Figure  24.  Equation  198  simply  states 
that  air  mass  between  radar  and  target  is  given  by  the  difference  between  the  total 
air  mass  and  the  air  mass  above  the  target. 

The  geometrical  relationships  between  target  altitude,  slant  range, 
and  elevation  angle  are  shown  by  the  set  of  dashed  curves  in  Figure  25.  The  curves 
are  labeled  by  target  altitude  in  kilon  oters,  the  ordinate  labels  are  slant  ranges  in 
miles,  and  the  abscissa  labels  are  elevation  angles.  The  top  of  t.he  effective  atmos¬ 
phere  is  considered  to  be  at  an  altitude  of  40  kilometers.  With  the  curves,  knov/- 
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AIR  MASS  (DIMENSIONIESS) 
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Eigure  23: 
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AIR  MASS  AS  A  FUNCTION  OF  ELEVATION  ANGLE 
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Top  of  the  AfiTiosphere 


THE  GEOMETRY  RELATING  TO  SPECIFICATION  OF  AiR  MASS 
AS  A  FUNCTION  OF  TARGET  ALTITUDE 


SQUARE  RESIDUAL  SLANT-RANGE  ERROR  NOMOGRAM 


ledge  of  any  two  range  parameters  permits  determination  of  the  third.  For  ex¬ 
ample,  the  slant  "ange  at  a  degree  elevation  angle  to  a  target  at  10-kilometers 
altitude  would  be  about  34  miles,  ^his  is  determined  by  fo  ’owing  the  10  degree 
abscissa  to  its  intersection  with  i.’  ^  10-’  iometer  target  altitude  and  noting  the 
ordinate  that  also  intersects  this  point,  similarly,  a  tar;  H  at  an  altitude  of  40 
kilometers  (the  top  of  the  atmosphere)  100  slant-range  mileb  away  would  be 
viewed  at  an  elevation  angle  about  13  depre  s.  "^his  ib  determined  by  follow¬ 
ing  the  100-mile  ordinate  to  its  inters  'oi  wim  the  40  kilomei  r  target  altitude 
and  noting  the  value  of  the  ansc'ssa  thai  'so  Intersec  s  the  point. 

The  geometi  oal  rt’ation  hips  defin  '  by  t.  'i  dashed  curves  in 
Figure  25  were  used  first  to  defii.e,  in  ccordan  ;  with  Figure  23  and  Equation 
198,  the  air-mast  value  to  be  assi,  ned  th,  ?lant  path  i*'  ny  direc'^ion  lo  any  alti¬ 
tude  .  Then,  invoh’ng  the  linear  .  el  tionsuip  bet  een  air  mass  and  N-measurement 
error,  the  defined  a^r  mass  values  'e  L.ultipliei  by  the  at  olute  ri  s  1  alues 
per  unit  air  mass  contained  in  Table  9  t  detern.ine  die  cui  resp  mding  residu  1 
slant -range  errors.  These  dat  i  are  plotted  as  ti^e  solid  curves  in  Figure  25. 

Figure  25  can  thus  be  used  as  a  nomogram  for  determining  the  i.  sidual  slant- 
range  error  associated  wit>  "'ny  receiver -to -target  geometery.  That  is,  the  inter¬ 
section  of  any  abscissa  with  .re  solid  and  the  dashed  target  altitude  dt  ermines 
the  slant  range  to  th>.  target  and  the  residual  error  that  woi'ld  result  if  hat  slant 
range  were  measured  with  a  radar  ana  corrected  with  radiometric  ally  acquired  N- 
data.  For  example,  suppose  a  target  at  an  altitude  of  40  kilometers  /ns  vi  wed 
at  a  10  degree  elevation  angle.  Th''  i.  tersection  of  th.  lu  degree  abscissa  "Ith 
the  dashed  40-kilometer  target  altitude  fixes  the  true  rang,,  at  ubout  130  mileb, 
and  Its  ^section  with  the  40 -kilometer  solid  curvet  es^U'blishrs  nhout  a  G.  2- 
fe.'t  I'esidual  slant -range  error.  Therefore  the  radar-rneabured  range  woulu  be 
S  =  130  miles  a^O.  2  feet. 

The  data  contained  In  Figure  25  were  replotted  'n  a  different 
fashion  in  Figure  Ordinate  labels  are  slant  range  In  miles,  al  scissa  labels 
are  rojs  residual  slant -range  errorb  In  feet,  and  the  family  of  cu  v  ‘s  is  labeled 
by  slant -range  elevation  angle.  For  example,  Uie  slant-range  lesirual  error 
associated  with  a  slant  range  of  30  miles  to  a  target  at  20  degrees  eu-vation  would 
be  about  0. 1  foot.  This  .'alue  is  obtained  by  following  the  30 -mile  slant -range 
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ordinate  to  its  intersection  with  Uie  20  degree  curve  and  noting  the  abscissa  that 
also  Intersects  this  point.  For  the  s^me  geometry,  the  same  residual  error  could 
have  been  obtained  from  Figure  25.  Figure  26  illustrates,  however,  that  once  the 
target  leaves  the  atmosphere,  the  residual  error  is  independent  of  slant  range. 

For  example,  a  target  viewed  at  a  20  degree  elevation  angle  100  slant -range  miles 
away  is  already  beyond  the  atmosphere  (i.e. ,  at  an  altitude  in  excess  of  40  kilo¬ 
meters).  The  residual  error  has  already  attained  its  maximum  value  (about 
0. 107  feei)  and  will  stay  at  that  value  regardless  of  the  increase  in  slant  range. 

It  must  be  emphasized  that  the  residual  slant-range  errors 
specified  in  Figures  25,  26,  and  27  can  only  be  considered  as  representative  of 
actuality  because  they  are  based  on  consideration  of  a  standard  air  mass.  To  re¬ 
capitulate;  a  standard  air  mass  has  been  defined  as  one  in  which  pressure,  temp¬ 
erature,  and  equivalent  thickness  are  as  defined  by  the  U.  S.  Standard  Atmosphere, 
and  as  cne  in  which  there  is  a  total  water  content  of  2.  0  precipitable  centimeters. 

A  real  atmosphere  that  departs  from  this  standard  wilt,  of  course,  generate  resi¬ 
dual  errors  different  from  the  ones  specified. 

It  is  further  emphasized  that  the  values  pretented  are  those  that 
would  be  realized  if  slant-range  data  were  corrected  by  using  the  hybrid-active- 
system  line  integral  refractometer.  Use  of  either  the  all -microwave  active  sys¬ 
tem  or  the  passive  system  would  result  in  residual  errors  about  2.  5  times  greater 
than  those  indicated  above.  For  the  latter  system,  the  realization  of  a  residual 
slant-range  error  only  2.  5  times  that  attainable  with  the  hybrid  active  system  would 
be  strictly  true  only  for  targets  entirely  outside  the  atmosphere.  Targets  imbed¬ 
ded  in  the  atmosphere  raise  the  problem  of  how  well  the  contribution  to  N  from 
the  atmosphere  in  front  of  the  target,  the  only  contribution  that  acts  to  establish 
A  (AS)  can  be  distinguished  from  that  behind  the  target  with  systems  of  limited  pro¬ 
filing  capability.  Of  course,  as  emphasized  In  Pa*  ograph  5.  B.  (2)(A),  the  more 
frequencies  used  in  the  profiling  radiometer,  the  more  detailed  the  measured  pro¬ 
file.  However,  such  a  philosophy  can  be  carried  just  so  far.  purely  on  the  basis 
of  economics.  With  just  a  four -channel  radiometer,  which  peniits  the  determin¬ 
ation  of  K  to  four  different  levels,  the  .issumptlon  simply  of  an  exponeatial'de- 
crease  of  line -integrated  refractivlty  with  height  between  levels  should  result  in 
no  more  than  a  lO-percert  increase  in  the  relative  error  in  N  estimated  to  imy 
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point  within  that  layer.  Such  an  additional  uncertainty  would  require  multiplying 
the  residual  errors  applicable  to  atmosphere -imbedded  targets  shown  in  Figures 
25,  26,  and  27  by  a  factor  of  three  as  opposed  to  2.  5. 

(B)  Elevation  angle 

To  determine  the  residual  elevation-angle  error  resulting  from 
the  residual  errors  in  slant -range  measurement,  the^-terms  of  Equation  197 
must  first  be  evaluated  as  dictated  by  Equations  11  and  12.  !^ecifieally,  this 
would  entail  evaluations  of  such  expressions  as 


A  bewildering  complexity  in  the  dependence  of  ^  {A^)  on  the  geometry  of  the 
problem  is  demonstrated.  That  is,  different  combinations  of  the  three  interfer¬ 
ometer  slant  ranges  —  S^,  S^,  and  — with  their  corresponding  slant  resi¬ 
duals  — )•  ),  aridZ^(AS  )  — would  result  in  different  values  for  A(A6). 

Ox  ^ 

The  specification  of  a  particular  combination  for  quantitative  evaluation  of  Equa¬ 
tion  197  is  strictly  aroitrary,  and  the  results  derived  must  be  accepted  as  being 
only  representative  as  to  the  order  of  magnitude  in  the4(A^-)  that  might  be 
realized  in  actual  ranging  operations. 

In  making  such  a  specification,  (strictly  as  a  matter  of  conven¬ 
ience  to  reduce  the  computational  complexity  of  the  problem)  consider  elevation- 
angle  measurements  to  locus  of  points  everywhere  equidistant  from  the  interfero¬ 
metric  slave  stations  and  at  a  constant  elevation  angle  from  the  mother  station. 

In  this  case  S  =  S  =  S,  <p  ~  (p  =^,  and  It  may  be  assumed  thatzl(^  ) 

1.  2  1  2  1  ’  2 

If  it  is  further  considered  that  B  »  B  »  B,  then  Equation  197  reduces  to 

X  A 

-  2.  5L^2(p  (202) 
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Furthermore,  Equation  18  may  be  written 


I  ~  (p  (iJ03) 

or 

COS^G  -2to5^C^  (204) 

Equation  202  thus  defines  the  elevation  angle  in  terms  of  the  interferometer  angle 
p'  as  defined  by  either  equation  11  or  12.  Equations  201,  202,  and  204  thus  permit 
a  parametric  study  to  determine as  a  function  of  slant  range  and  slant -rarge 
residual  for  the  geometry  specified.  The  results  for  B  =  20  miles  are  shown  in 
Figure  28.  Ordinate  labels  are  slant  range  in  miles  as  measured  from  the  mother 
station;  the  abscissa  labels  are  residual  elevation  angles  in  microradians.  The 
family  of  curves  is  labelled  by  elevation  angle.  The  data  contained  in  Figure  28  is 
even  less  general  than  those  contained  in  Figures  25,  26,  and  27  because  not  only 
does  the  standard  air  mass  restriction  apply  in  this  case,  but  the  special  geometry 
invok  i  for  the  parametric  study  must  also  be  kept  in  mind.  In  spite  of  these  re¬ 
strictions,  it  is  still  safe  to  state  that  interferometrically  determined  elevation 
angles  will  be  correct  to  orders  of  tens  of  microradians  if  the  determining  slant- 
range  data  are  corrected  through  radiometric  N  -measurement  to  the  values  shown 
in  Figures  25,  26,  and  27. 

E,  Summary 

Regarding  the  determination  of  N  from  spectrographic  measurements  of 
integrated  oxygen  and  water -vapor  absorption  coefficients,  two  factors  cause  un¬ 
certainties  in  the  final  N  -specification.  The  first  is  related  to  uncertainties  in  the 
refractivity  equation  itself;  the  second  is  related  to  spectrographic  measurement 
errors.  The  former  effect  is  negligible  compared  to  the  latter. 

Spectrographic  measurement  errors  stem  from  (1)  uncertainties  in  the 
functional  relationships  describing,  for  any  portion  of  the  electromagnetic  spectrum, 
transmittance  as  a  function  of  bandwidth;  (2)  pressure  and  temperature  line -broaden¬ 
ing  effects;  (3)  continuum  extinction  effects;  (4)  scintillation  effects;  and  (5)  effects  of 
errors  in  mean-temperature  specification. 
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Item  1  above  establishes  the  fundamental  measurement  error,  and  its  effect  is 
felt  regardless  of  the  radiometric  systers  configuration.  The  remaining  items 
are  classed  as  atmospheric  effects.  Of  these  latter  items,  Items  2,3,  and  4 
affect  spectrcgraphic  measurements  performed  with  the  hybrid  active  system; 
items  2,  3,  4  and  5  apply  to  the  all-microwave  active  system;  items  2,3,  and  5 
apply  to  the  passive  system. 

The  reference -to -sensing  signal  ratio  mode  of  operation  greatly  reduces 
the  effect  of  Items  2  and  3  on  mec.surement  error.  In  descending  order: 
line -broadening  and  me  an -temperature  specification  effects  result  in  the  greatest 
percentage  of  error  in  the  ultimate  N -specification;  continuum  extinction  effects 
are  next  in  importance;  scintillation  effects  can  be  dismissed  from  consideration 
entirely.  Under  typical  operating  conditions,  the  total  effect  of  all  error  contri-’ 
butors  is  that,  with  the  hybrid  system,  N-specification  will  be  correct  to  from 
0, 5  to  1.  0  percent,  and  with  the  all-microwave  active  or  passive  system  the  N  - 
error  will  amount  to  from  1. 0  to  2. 0  percent  of  N,  These  N-errors  result  in 
hybrid-system  slant -range  errors  of  about  1. 0  feet  rms  for  targets  outside  the 
atmosphere  viewed  at  a  1  degree  elevation  angle.  This  error  will  be  Increased  by 
about  a  factor  2. 5  wlth’the  all-microwave  or  passive  system.  The  associated  re¬ 
sidual  elevation-angle  errors  will  be  on  the  order  of  tens  of  microradians. 

Residual -error  magnitudes  decrease  with  decreasing  slant  range  for 
targets  imbedded  in  the  atmosphere.  In  this  case,  however,  the  passive  system 
accuracies  will  be  degraded  by  an  additional  10  percent  over  that  realized  by  the 
all-microwave  system,  due  to  limited  passive -system  profiling  capabilities. 


Section  II 

INSTRUMENTATION 


1.  BACKGROUND 

The  considerations  relating  to  the  radiometric  measurement  of  line -integrated 
refraetivity  discussed  in  Section  I  served  to  establish  basic  design  features  for  the 
various  line-integral  refractometer  configurations.  Specifically,  measurements 
of  received  energy  intensities  in  absorptions  bands  must  be  referenced  to  similar 
measurements  performed  in  nearby  atmospheric  windows;  this  is  done  to  minimize 
degradation  of  the  acquired  data  by  continuum  extinction  effects.  Furthermore, 
the  active-system  measurements  must  be  performed  simultaneously  and  through 
a  single  aperture  to  eliminate  scintillation  effects.  With  these  general  features 
in  mind,  the  follo'ving  instrument  designs  are  offered  as  guides  to  ultimate  instru¬ 
ment  development. 

2.  INFRARED  ACTIVE  SYSTEM 

A.  Boeing  system 

(1)  Description 

An  infrared  system  incorporating  the  basic  design  feabares  listed 
above  has  been  developed  and  qualitatively  field  tested  by  The  Boeing  Company 
(Reference  63).  Briefly,  energy  from  a  continuous  spectrum  source  is  collected 
at  a  receiver  through  a  single  aperture.  The  field  is  collimated  with  a  reflecting 
optical  system.  Temperature -compensated  lead  sulfide  cells  are  used  to  sense 
the  intensity  of  the  collected  energy  In  three  bands  of  the  infrared  spectrum : 
one  centered  at  a  wavelength  of  1.27  microns,  the  second  at  1.73  microns,  the 
third  at  2. 12  microns.  The  first  wavelength  lies  in  the  oxygen  band  previously 
identified,  the  second  is  in  the  wing  of  the  1.86-mlcron  water-vapor  band,  and  the 
third  is  an  atmospheric  window.  Wavelengths  are  isolated  with  narrow-band 
germanium  interference  filters.  The  required  band-energy  ratios  are  effec¬ 
tively  determined  by  passing  the  appropriate  detector  outputs  through  logarithmic 
signal  converters  and  then  through  a  differential  c  iplifler  to  a  recorder.  A 
more  detailed  description  of  the  system  is  given  in  Appendix  1 . 
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(2)  Performance 

Examples  of  data  collected  with  the  Boeing  infx’ared  system  are 
as  shown  in  Figures  29,  30,  and  31,  Shown  in  Figure  29  are  simultaneous  record¬ 
ing  of  dry -bulb  temperatures  (the  left-hand  trace)  and  of  line -integrated  oxygen 
density  (the  right-hand  trace)  measured  over  a  2 -kilometer  propagation  path. 

The  scale  factors  are  as  Indicated  on  the  figure.  ITote  that  temperature  increases 
to  the  right,  and  ^  increases  to  the  left.  The  latter  recorder  response  was 
selected  because  density  is  inversely  proportional  to  temperature.  Therefore, 
a  given  change  in  the  temiierature  trace  will  be  reflected  in  a  corresponding 
change  in  the  ojiygen  trace.  The  general  trends  of  the  two  traces  are  well  cor¬ 
related,  but  the  smoothed  effect  in  the  integrated  value  is  also  readily  apparent. 
This  effect  demonstrates  conclusively  why  estimates  of  line -integrated  refractivity 
from  point-measurement  data  are  so  uncertain. 

___  Figure  30  again  shown  simultaneously  recorded  temperature  and 

^  traces,  along  with  a  wind  direction  and  speed  record.  The  notable  feature 
is  the  correlation  of  the  "humps"  in  the  bottom  portions  of  all  four  traces.  The 
change  in  wind  direction  and  speed  obviously  heralded  warm -air  advection 
(as  evidenced  in  the  temperature  trace)  — a  circumstance  tliat  is  also  reflected 
in  the  trace.  Again,  the  integrated  effect  is  not  as  great  as  the  single - 
point  effect. 

Figure  31  shows  simultaneous  dry -bulb  and  integrated  water - 
vapor  density  traces.  As  expected,  very  little  correlation  is  noted  b(jtwfcon  these 
two.  The  traces  were  obtained  in  early  morning  after  a  foggy  night.  Therefore, 
the  slight  warming  trend  indicated  at  the  bottom  of  the  temperature  trace  was  re¬ 
flected  in  the  integrated  vapor  density,  or  absolute  humidity,  trace.  The  effect 
is  attributable  to  the  Increased  capacity  of  the  air  to  hold  more  water  vapor  at  the 
increased  temperature. 

The  purpose  of  the  tests  conducted  during  which  the  data  shown 
in  Figures  29,  30,  and  31  were  acquired  was  not  one  of  instrument  calibration, 
but  rather  to  establish  the  basic  worth  of  the  instrument  philosophy  and  design. 

For  this  purpose  the  tests  were  successful.  That  is,  it  was  conclusively 
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demonstrated  that  (1)  integrated  oxygen  and  water-vapor  density  data  are  readily 
attainable  by  spectrographic  means,  and  (2)  that  the  serusing-to-refe  ente  band 
ratio  mode  of  operation  serves  effectively  to  negate  the  potentiallv  adverse  efieets 
of  continuum  extinction  and  scintillation  on  data  quality. 

B.  Range  operations  system 

The  Boeing  infrared  system  just  described  defines  the  basic  design  for 
the  development  of  the  system  applicable  to  range  ojxirations.  For  sach  oikum- 
tions,  however,  it  is  tentatively  suggested  that  the  sarr.oling  wavelengths  be  1.2(i  } 
microns  for  the  oxygen  band,  1.32  microns  for  the  water-vapor  band,  and  1.22 
microns  for  the  window.  Such  relocation  of  the  water-vapor  sensing  wavelength 
will  likely  ensure  transmittances  high  enough  for  measurable  signals  even  at  the 
extreme  ranges  likely  to  oe  encountered  in  ranging  operations.  The  relocation 
of  tlic  window  wavelength  is  for  closer  grouping  of  the  three  wavelengths  for  the 
maximum  compensation  for  continuum  extinction  effects  (see  Kquation  158). 

The  system  can  be  further  upgraded  by  (1)  increasing  the  size  of  the  collecting 
optics,  and  (2)  using  cooled  radiation  detectors.  The  former  procedure  re.sults 
in  the  direction  of  a  greater  amount  of  energy  onto  the  detector,  the  latter  in 
increased  detector  detectivities  (Reference  61). 

C,  Expected  performance 

The  power  received  by  a  detector  from  a  remote  source  emitting  energ)- 
in  the  wavelength  interval  can  be  estimated  from 


P 


s 


J. 


(205) 


where  P  is  the  received  pwwer,  the  radhant  intensity  of  the  source,  .A-p  thi 
area  of  the  target,  Appj  the  aperture  area  of  the  collecting  optics,  /''^t  the 
transmittance  of  the  atmosphere,  /^opt  transmittance  of  the  collecting  optical 
system,  zAA  U'.c  wavelength  interval  of  the  band  pass  filter,  and  S  the  slant  rangi- 
Consider  the  following  parameter  values; 

2  -2  -1 
c-  5  X  10  watts  cm  sterad.  micron 


4  X  10^  cm^ 
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.  ^  ‘  10^  cm^ 

,,.-1 

_2 

j  X  10  microns 
1.4  X  10^^  cm^ 


The  above  set  of  parameters  relate  to  a  3 -inch -diameter  3000“C  target 
viewed  with  a  24 -inch -diameter  telescope  at  an  elevation  angle  of  1  degree  and 
at  a  slant  range  of  about  380  kilometers  (the  distance  to  the  top  of  the  atmosphere). 
Under  these  conditions  the  power  received  is  on  the  order  of  10~12  watts.  The 
detection  of  such  low  signal  levels  is  within  the  capability  of  current  cooled - 
detector  technology.  Taus,  on  an  average  "clear"  day,  fhs  infrared  system 
could  continue  to  fxmction  until  tracked  targets  left  the  atmosphere.  Thereafter, 
operation  would  switch  to  the  ail -microwave  system.  On  excessively  hazy  days, 
the  switch  to  all-microwave  operation  might  have  to  be  effected  sooner. 

3.  MICROWAVE  ACTIVE  SYSTEMS 

The  philosophy  of  the  microwave  active  systems  is  the  same  as  the  active 
infrared  system;  that  is,  received  powers  in  oxygen  and  water -vapor  bands  must 
be  referenced  to  power  levels  received  in  atmospheric  windows.  The  requirement 
is  for  dual -frequency  beacons  and  dual -frequency  receivers.  For  the  water-vapor 
measurements,  it  is  therefore  suggested  th^  3.5  GHz  be  used  as  the  reference 
frequency  and  17  GHz  as  the  sensing  frequency.  This  sensing  frequency,  lying 
in  the  wing  of  the  22 -GHz  water -vapor  band,  is  harmonically  related  to  the 
reference  frequency;  thus,  frequency  multiplying  techniques  can  be  used  in  both 
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beacon  and  receiver  to  simplify  instrument  design.  Similarly,  for  the  oxygen 
measurements,  11-GHz  can  be  used  as  the  reference  frequency,  and  44-GHz 
(in  the  wing  of  the  60-GHz  oxygen  band)  can  be  used  as  the  sensing  frequency. 
Here,  too,  the  reference  and  sensing  frequencies  are  harmon  »cally  relate^ 
again  leading  to  simplification  in  instrument  design. 

Both  the  beacon  and  receiver  designs  in  both  the  44-GHz  and  the  l7-GHz 
systems  have  similar  design  concepts.  The  beacons  will  differ  only  in  the  rf 
transmitter  power  requirements,  and  the  receivers  will  differ  only  in  mixers 
and  local  oscillators.  ITierefore,  beacon  and  receiver  designs  can  be  presentrd 
as  applicable  to  both  systems. 

A.  Dual -frequency  beacons 

The  basic  beacon  design  is  e^own  in  Figure  32.  In  general,  both  dual¬ 
frequency  beacons  will  be  pulse -modulated  units  featuring  solid-state  signal 
generation.  In  the  8.5/17-GHz  system,  signals  can  be  amplified  with  ti’aveling- 
wave  tube  amplifiers  (TWT);  in  the  11/44-CHz  system,  signals  cui  be  amplified 
with  pulsed  hi^-power  magnetrons.  Carrier  frequencies  can  be  derived  by 
using  a  veractor  multipl3ring  chain  to  up-convert  a  low-frequency  crystal- 
controlled  oscillator  output.  Hig^y  stable  local  oscillators  will  serve  both 
frequencies  in  each  system.  In  the  water-vapor  system,  a  frequency  doubler 
will  serve  to  boost  the  8.5-GHz  reference  frequency  to  the  17 -GHz  sensing 
frequency. 

The  modulators  will  be  transistorized  and  will  contain  oscillators ,  squaring 
circuits,  switch  drivers,  and  switches.  The  switches  will  be  capable  of  switching 
1-watt  power  levels  with  an  Insertton  loss  less  than  1  db,  on/off  isolation  of  40  db, 
and  switching  rise  times  of  10  nanoseconds.  Presently  available  diode  switches 
reliably  satisfy  these  requirements.  The  remaining  modulator  functions  are  not 
significant  design  impositions,  and  standard-type  circuits  will  be  used. 

The  bandpass  filters  can  be  of  the  waveguide ,  using  Tchebysheff  techniques 
to  achieve  minimum  ripple  and  insertion  loss  vdiile  retaining;  the  desired  cutoff 
characteristics.  Loads  mus* '' ''  equalized  at  the  filter  output  and  input  to  ensure 
minimum  mismatch  losses. 
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B .  Dual-fre<r>ency  receivers 

The  receiver  shown  in  Figure  33  is  a  two-channel,  superheterodyne, 

PRF -gated  subsystem,  PRF  gating  can  be  achieved  by  switching  off  the  LO 
signal  for  the  PRF  period  after  receiving  a  pulse.  At  the  end  of  the  period, 
the  "LO"  signal  comes  back  cn  and  the  receiver  is  again  on  for  the  next  pulse. 

A  delay  of  100  nanoseconds  should  be  incorporated  to  allow  the  full  pulse  to  be 
processed.  Boxcar  detectors  can  be  used  to  retain  the  data  during  the  interpuise 
periods.  Provision  should  also  be  made  for  observing  with  a  sampling  scope  the 
data  at  the  output  of  the  IF  amplifier  to  determine  the  influence  of  multiple  trans¬ 
mission  paths  or  noise. 

Bandwidths  of  the  receivers  should  be  wider  than  those  of  the  trans¬ 
mitters  so  that  differences  in  the  pulse  slqpes  between  short  and  long  transmis¬ 
sion  may  be  observed.  PRF  gating  can  be  used  to  remove  multipaths  via  antenna 
side  lobes  and  also  to  reject  environmental  noise  between  pulses. 

C.  Expected  performance 

Ihe  frequencies  specified  for  both  active  microwave  systems  were  chosen 
to  ensure  large  enough  atmospheric  transmittance  values  for  workable  received 
power  levels  at  low  elevation  angles  and  extreme  slant  ranges.  It  is  expected  that, 
even  for  targets  beyond  the  atmosphere  at  1  deeree  elevation  angle,  the  oxygen 
absorption  at  will  result  in  signal  losses  of  order  lOdb,  Thus,  in 

spite  of  inverse -square -law  radiation -spreading  losses,  targets  may  be  tracked 
to  much  greater  slant  ranges  with  the  microwave  systems  than  with  the  infrared 
systems,  because  microwaves  do  not  suffer  from  attenuation  by  atmospheric  haze. 

Because  the  wavelength  region  in  which  the  active  systems  operate  is  not 
represented  by  suitable  cw  power  sources,  the  systems  are  designed  for  pulsed 
operation.  Consequently,  the  designs  specified  in  Figures  32  and  33  largely 
incorporate  components  well  within  the  capability  of  current  microwave  technology. 
With  gain  stabilization  thereby  achieved  (on  the  order  of  10 db),  and  with  hori¬ 
zontal  polarization  to  minimize  forward  scattering  effects,  transmission  measure¬ 
ments  correct  to  1  percent  rms  should  be  achieved.  The  contribution  to  N  error 
due  to  instrument  noise  is  therefore  not  significantly  greater  than  that  established 
by  atmospheric  effects  (Table  8). 
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4.  PASSIVE  SYSTEMS 


The  radiometric  receiver  designs  at  20  and  53  GHz  are  basically  Dicke 
systems  (Reference  65)  using  superheterodyne  techniques.  At  20  GHz,  a  three- 
channel  receiver  system  with  a  window  reference  frequency  is  specified,  while 
at  53  GHz  a  four-channel  receiver  with  reference  is  required.  For  both  the  20- 
and  53-GHz  systems,  each  channel  will  be  capable  of  measuring  brightness 
temperature  to  ITK. 

A.  Oxygen  system 

The  design  of  the  53-GHz  radiometer  is  shown  in  Figure  34.  This 
system  is  designed  with  tunnel -diode  amplifiers  as  the  IF  amplifiers  so  that  two 
or  more  bands  can  be  handled  simultaneously.  Duplication  of  some  electronic 
components  can  be  avoided  with  this  technique. 

For  this  53-GHz  region,  the  noise  comparison  source  is  injected  at  the 
IF  because  ferrite  switches  at  53  GHz  would  only  have  isolation  in  the  oft-state  of 
approximately  20  db.  The  diode  switch  will  give  isolation  of  60  db  or  more  and 
will  not  modulate  the  local  oscillator.  The  system  will  require  accurate,  known 
conversion  loss.  Single -sideband  reception  in  the  superheterodyne  receiver  is 
necessary  so  that  cnly  the  absorption  frequency  region  about  53  GHz  is  selected. 
The  IF  is  selected  at  X-band  to  support  the  two  individual  channels  simultaneously 
because  the  300-MHz  requirement  on  the  bandwidth  is  somewhat  restrictive  and 
four  frequency  bands  cannot  be  supported  in  one  humel-diode  amplifier.  Although 
a  traveling-wave  tube  (TWT)  could  theoretically  support  the  four  bands,  the  noise 
figure  would  be  higher  and  the  gain  stability  worse.  Previous  studies  have  shown 
that  tunnel-diode  amplifiers  show  excellent  short-term  (time  ^  1  second)  gain 
stability. 

The  traveling-wave  tube  is  specified  as  the  output  stage  because  tunnel  diode 
amplifiers  saturate  at  a  level  given  in  dbm  as  1  db  compression  point  =  -(2Cj+10) 
dbm,  where  G  =  gain  in  db.  If  three-bands  are  amplified  in  the  system  simultan¬ 
eously,  the  compression  point  will  be  lowered  10  db  minimum.  Consequently,  the 
system  will  be  designed  to  enter  the  TWT  at  or  near  a  -60dbm  signal  level.  By 
properly  selecting  waveguide  size,  single -sideband  operation  can  be  achieved  by 
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operation  of  the  lower  sideband  in  a  v/aveguide  beyond  cutoff.  This  form  of  filtering 
will  modulate  the  impedance  of  the  local  oscillator,  since  filtering  is  accomplished 
before  the  signal  enters  the  mixer  and  does  not  affect  the  desired  signal.  The  in¬ 
sertion  loss  resulting  from  an  inc  reased  waveguide  length  for  the  signal  will  be 
much  less  than  that  experienced  in  a  lumped  filter  network.  P'igure  35  shows  the 
general  outline  of  this  technique.  Tentatively,  operating-frequency  bands  have 
been  specified  as  52.3,  52.8,  53.3,  and  54.8  GHz,  with  a  bandwidth  of  300  MHz 
per  band.  The  three  first  channels  are  handled  in  a  three-channel  radiometer. 

The  54.  S-GHz  band  must  be  handled  in  a  separate  channel  because  of  the  band¬ 
width  lirritations  of  tunnel -diode  amplifiers.  However,  the  second  channel  has 
the  capability  to  study  two  additional  bands  and  make  a  six -band  radiometer  if 
subsequent  need  arises.  The  bandwidth  of  the  system  can  handle  three  bands 
and  provide  amplification  with  Intermodulation  products  of  ^  50  db. 

The  system  is  specified  to  have  an  accuracy  ATg  =  1®K,  and  from  Strom 
(Reference  66) 


where 

Tg  =  System  noise  temperature 
W  =  Prediction  bandwidth  —  300  MHz 
t  =  Look  time 

Ihe  system  noise  temperature  will  be  approximately  SOOOIC.  Using  the  above 
data,  t  =  0.375  second. 

Field  results  usually  show  at  least  a  2:1  field  degradation  in  jd  Tg  measure¬ 
ments.  Consequently,  the  look  time  will  be  1  second  or  longer.  The  postdetection 
amplifier  must  be  capable  of  handling  low  frequencies  with  good  fidelity  and  low 
distortion.  Because  this  Is  such  an  Important  technical  area,  an  experimental 
study  was  undertaken  to  investigate  design  features.  Figure  36  shows  the  schematic 
of  a  prototype  postdetection  amplifier;  test  results  are  shown  in  Figure  37.  Re¬ 
presentative  sizes  for  this  prototype  postdetection  amplifier  are  1.5  by  0.6  by 
3.5  inches. 
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B.  Water-vapor  system 

The  18-GHz  water-vapor  system  is  shown  In  Figure  38.  It  differs  from 
the  53-GHz  oxygen  system  in  that  each  band  hasa  separate  receiver.  The  20  CHz 
system  uses  transistor  amplifiers  in  the  IF  to  supply  the  10-MHz  bandwidth  with 
the  operating  frequency  of  about  500  MHz.  Sideband  or  image  rejection  is  obtained 
by  using  a  low  pass  filter,  with  the  cutoff  at  the  local  oscillator  frequency,  mount¬ 
ed  on  the  antenna  side  of  the  local  oscillator. 

The  local  oscillator  is  tentatively  specified  to  be  a  solid-state  unit  with 
three  frequencies  tentatively  specified  as  19,  20,  and  21  GHz,  developed.  The 
switching  elements  are  specified  as  semiconductor  switching  diodes  to  reduce 
leakage  because  ferrite  switches  do  not  provide  high  isolation  at  the  higher  fre¬ 
quencies. 

For  the  individual  bands  in  Figure  38,  the  bandwidth  will  be  lOO  MHz , 

Tg  =  2600*K,  T3  =  1*K,  and  the  look  time  is  t  =  0.5  seconds,  which  again  wiU 
be  boosted  to  1  second  or  more  to  achieve  the  A  Tb  -  1“K. 

C.  Addltiog*'.!  comments 

Both  the  20-  and  53-GHz  systems  are  cut  off  at  the  detector,  unless 
e.xperience  gained  th  operating  systems  indicates  the  advisability  of  continuous 
readings  in  all  bands.  It  may  be  sufficient  only  to  scan  the  bands,  but  this  carnot 
be  ascertained  until  the  sj stems  are  field  tested.  Insofar  as  possible,  both  sys¬ 
tems  will  be  designed  around  solid-state  components  to  Increase  rell.'ibllity  and 
reduce  power  requirements. 

The  window  radiometer  is  identical  for  both  systems.  It  must  be  a 
separate  instrument.  ‘P'C  instrumentation  for  this  receiver  can  be  an  AlL  Type 
2392 -B  universal  radiometric  system,  which  provides  a  cost-effective  method 
of  obtaining  the  window  band  with  a  minimum  of  design  and  test  efforts.  Although 
some  of  the  system  a8p.M:ta  outlined  above  may  seem  overly  complex,  auch  com¬ 
plexity  is  necessary  to  achieve  desired  results. 

D.  Expected  perfomnance 

The  four-channel  oxygen  system  will  permit  four-layer  integrated  oxygen 
density  profiling,  whereas  the  three-channel  water-vapor  system  will  permit  only 
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a  three-layer  integrated  water-vapor  density  profile.  However,  the  vertical 
extent  of  die  effective  water -vapor  atnoMphere  is  much  less  than  that  of  oxygen; 
therefore,  three  layers  are  considered  sufficient. 

An  integration  time  of  1  s'xsond  is  required  to  effect  the  measurement  of 
brightness  temperature  for  each  layer;  thus,  some  of  the  real-time  capability 
is  lost.  However,  die  values  of  ^  and  P  to  be  inferred  from  such  measure - 

Je>^  Jy**r 

ments  are  applicable  to  layers  of  finite  vertical  extents ,  and  must  be  interpreted 
as  average  values  for  these  layers.  The  finite  times  taken  by  targets  in  moving 
throu^  such  layers  lessens  the  impact  of  the  lost  real-time  capability. 

Both  systems  will  provide  brightness  temperatures  correct  to  ±1*K,  and 
accuracies  in  integrated  absorption  to  better  than  1  percent.  The  further  contri- 
bution  to  N  error  from  instrument  noise  is  again  of  minor  importance. 

5.  SUMMARY 

The  three  different  types  of  radiometric  systems  required  to  satisfy  all  re¬ 
quirements  for  line -integrated  measurement  in  all  ranging  (^rations  are  within 
the  current  capability  of  both  infrared  and  microwave  technology.  Also,  systems 
can  be  designed  and  constructed  to  perform  the  N  measurements  in  real  time  to 
the  precision  allowed  by  the  Umifing  effects  of  atmospheric  noise.  Furthermore, 
Insofar  as  the  measurements  are  performed  in  real  time,  no  data  processing  is 
required,  and  the  information  obtained  can  be  introduced  directly  into  the  tracking 
radar  systems  readouts  so  that  the  corrected  range  parameter  values  are  the 
direct  outputs  of  the  radar  systems. 

The  designs  of  all  the  systems  just  presented  — hybrid  active,  all -microwave 
active,  and  passive  — must  be  accepted  as  tentative  engineering  model  designs 
rather  than  as  true  prototype -system  designs.  No  design  can  be  called  final  be¬ 
cause  there  is  too  much  empiricism  in  the  "lews"  invoked  to  describe  the  trans¬ 
mission  of  energy  through  the  atmosphere.  Such  final  design  can  come  only  from 
e^qperimentation  with  the  engineering-model  instrumentation. 
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Section  111 


CONCLUSIONS  AND  RECOMMENDATIONS 

1.  CONCLUSIONS 

A.  Range -data  correction 

The  interferometric  or  trilateraticm  mode  of  ranging  operations  isolates 
target  slant  range  as  the  basic  range  parameter.  Thus,  ray  retardation  is  the 
atmospheric  effect  that  causes  errors  in  ranging  data.  Accurate  evaluation,  from 
raal-time  measurement  of  line -integrated  refractivity,  K,  is  therefore  mandatory 
for  precision  ranging  operations. 

Estimates  of  N  from  point  measurements  of  refractivity  do  not  satisfy 
the  real-time  precision  requirements  of  ranging  operations.  The  problem  is  the 
inadequacy  of  the  representation  of  an  integrated  function  as  ^  finite  siun .  The 
direct  measures  of  the  integrated  function  from  xay  transit  time  differences  at 
two  radiation  wavelengths  — established  by  (1)  normal  dispersion  at  optical  wave¬ 
lengths,  and  (2)  anomalous  dispersion  across  atmospheric  constituent  absorption 
bands  at  radio  wavelengths  — also  do  not  provide  the  complete  answer.  The  nor- 
rnal-dispersion  technique  is  Inadequate  because  the  data  so  obtained  can  be  applied 
only  to  optical  ranging  devices.  The  anomalous-dispersion  technique  cannot  pro¬ 
vide  a  true  real-time  capability  in  N  measurement  because  of  excessive  data- 
processing  averaging  times.  Also,  both  dispersion  techniques  require  cooperative 
targets  in  their  implementation. 

B.  Solution 

Both  absorption  and  refraction  of  radiaticHi  proceed  in  proportion  to  atmos¬ 
pheric  compositicui,  thereby  providing  the  basis  for  an  alternate  approach  to  the 
measurement  of  N.  That  is,  at  optical  frequencies  ^corresponding  to  wavelengths 
less  than  20  microns)  the  form  of  the  atmospheric  input  parameter  to  absorption 
calculations  is  precisely  that  of  the  electronic  charge  distortion  terms  (as  Identi¬ 
fied  by  the  first  two  of  the  three-term  Smith -Weintraub  refractivity  equatiw)  of 
line -integrated  refractivity;  at  microwave  frequencies  the  form  of  the  absorption 
parameter  matches  that  of  tlie  N  water -vapor  dipole  orientation  term. 
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Radiation -absorption  measurements  performed  in  oxygen  and  water-vapor 
absorption  bands  at  infrared  frequencies,  and  in  a  water-vapor  band  at  microwave 
frequencies,  will  thus  provide  the  date,  required  for  calculating  line -integrated  re- 
fractivity  as  dictated  by  the  integrated  form  of  the  Smith -Weintraub  refractivity 
equation. 

TTiree,  distinct,  spectrographic  line-integral  refractometers  are  proposed 
to  satisfy  the  demands  of  all  types  of  ranging  operations.  For  missile  test-range 
iterations  — which  are  characterized  by  (1)  the  availability  of  cooperative  targets, 
and  (2)  conformance  of  test  schedules  to  ideal  weather  conditions  —  the  combined 
irfrared/microwave  active  system  is  offered  as  one  thai  provides  the  capability 
for  the  most  accurate  N  determinatiuus.  The  1.27-micron  oxygen  and  the  1.37- 
micron  water-vapor  band  maybe  used  for  the  infrared  measurements;  the  22-GHz 
water-vapor  band  may  be  used  for  the  microwave  measuremt  nt. 

To  extend  the  capability  of  the  active  stystem  to  all-weather  operations, 
an  all-microwave  system  is  envisioned.  The  22-GHz  water -vapor -band  measure¬ 
ments  can  be  used  to  evaluate  both  water-vapor  terms  of  the  N  equation.  Absorp¬ 
tion  measurements  in  the  55 -GHz  oxygen  band  may  be  used  in  place  of  the  1.27- 
micron  oxygen  measurement. 

A  passive  mode  of  operation  is  suggested  for  noncooperative  targets. 

Ill  this  mode  the  absorption  measurements  are  equivalently  performed  by  meas¬ 
uring  thermal  emissions  of  the  atmosphere  in  the  55 -GHz  oxygen  and  22-GHz 
water-vapor  bands.  Also,  multifrequency  operation  in  each  band  will  provide  a 
capability  for  a  degree  of  N  profiling 
C.  Expected  results 

By  referencing  absorption -band-radiation  intensity  measurements  to 
similar  measurements  performed  in  nearby  atmospheric -window  regions  of  the 
absorption  spectrum,  it  is  anticipated  ^hat  line -integrated  refractivity  may  be 
correctly  estimated  to  within  0.5  to  1.0  percent  when  using  the  hybrid,  combined  • 
infrared/microwave  line-integrai  refractometer.  With  the  all-microwave  sys¬ 
tem,  N  will  be  accurate  to  within  1.0  to  2.0  percent.  An  additional  0.5-percent 
error  will  arise  when  the  passive  microwave  system  is  applied  to  tracking  tar¬ 
gets  imbedded  in  the  atmosphere.  Such  accuracies  in  N  measurement  imply  that. 


137 


with  the  hybrid  system ,  residual  errors  in  slant -range  measurements  to  targets 
outside  the  atmosphere  viewed  at  a  1-degree  elevation  angle  will  be  about  1. 0- 
foot  rms.  With  the  other  systems  proposed,  the  slant-range  residual  error  will 
increase  by  the  factor  of  increase  in  the  basic  N  uncertainty.  The  associated 
residual  elevation-a  :-^le  errors  will  be  on  the  order  of  tens  of  microradians. 

D.  Instrumentation 

The  three  different  types  of  radiometric  systems  needed  to  satisfy  re¬ 
quirements  for  line -integrated  refractivity  measurements  in  all  ranging  operations 
are  within  the  current  capability  of  both  infrared  and  microwave  technology.  More¬ 
over,  systems  can  be  designed  to  perform  N  measurements  in  real  time  to  the 
precision  allowed  by  the  limiting  effects  of  atmospheric  noise. 

2.  RECOMMENDATIONS 

The  designs  for  the  three  basic  N  measurement  systems  — hybrid  active, 
all -microwave  active,  and  passive  — presented  in  this  study  are  not  true  proto¬ 
type  designs,  but  must  be  accepted  as  tentative  engineering -model  designs. 

There  is  too  much  empiricism  in  the  "laws"  invoked  to  describe  the  transmission 
of  energy  through  the  atmosphere  for  any  design  to  be  set  forth  as  final. 

It  is  recommended,  therefore,  that  the  instrument  designs  presented  be  im¬ 
plemented  as  engineering  models,  and  that  such  models  then  be  evaluated  to  de¬ 
termine  the  advisability  of  proceeding  to  the  true  prototype  system »  Such  an 
evaluation  will  necessarily  involve  two  aspects:  calibration  and  performance, 

A.  Engineering-model  calibration 

Because  of  the  uncertainty  that  exists  in  the  various  semiempeiical  ab¬ 
sorption  'laws"  invoked  to  describe  the  atmospheric  absorption  process  within 
different  atmospheric  constituent  absorption  bands  and  at  various  radiation  wave¬ 
lengths  (SectloD  I,  Paragraph  6.  C.  (1))  it  is  Impossible  to  perform  a  theoreUcml 
iistrument  calibration  to  any  degree  of  precision.  The  calibration  procedure 
must  consist  of  compering  the  instrument  outputs  against  corresponding  measure¬ 
ments  of  the  quantities  of  oxygen  and  water  vapor  that  exist  in  prescribed  sensing 
oaths.  A  calibration  procedure  using  laboratory  procedures  is  not  feasible  be¬ 
cause  of  the  excessive  path  lengths  required  to  duplicate  the  typical  quantities  ct 


oxygen  and  water  vapor  that  will  exist  in  the  atmospheric  paths  realized  when 
tracking  targets  beyond  the  atmosphere  at  low  elevation  angles.  At  a  1 -degree 
elevation  angle,  for  example,  such  an  atmospheric  path  would  be  equivalent  to 
about  a  200-kilometer  sea-level  path.  Also,  the  actual  form  of  the  absorption 
law  changes  as  the  path  length  (optical  depth)  increases  beyond  a  certain  undefin- 
able  limit.  The  long-path  calibratiai  cannot,  therefore,  be  expressed  as  a  sim¬ 
ple  multiple  of  the  short -path  calibration. 

(1)  Active -system  calibration 

Active -system  calibration  will  consist  of  comparing  simultaneous 
measurements  of  absorptlcm-band  transmlttances  and  partial  densities  of  oxygen 
and  water  vapor.  The  problem,  however,  is  that  conventional  meu^orological 
instruments  performed  the  measures  of  and  at  a  single  point,  whereas 
the  transmittance  measurement  represents  a  space -integrated  (i.e. ,  over  the 
prc^gatlon  path)  effect.  In  a  perfectly  homogeneous  atmosphere  there  would  be 
no  problem  in  estimating  integrated  partial-density  values  from  point  measure¬ 
ments  because  the  former  values  would  simply  constitute  the  latter  measurements 
multiplied  by  the  propagation -path  length. 

Complete  homogeneity  in  actual  test  situations  cannot  be  expected. 
Careful  test -site  selection  (e.g. ,  an  over-water  propagation  path  with  an  over¬ 
water  upwind  fetch)  can  result  in  the  closest  possible  approach  to  the  Ideal ,  but 
a  one-to-one  correlattcm  between  the  integrated  partial-density  value  and  the 
single -point  measurement  simply  cannot  be  expected.  Therefore,  the  calibration 
curve  must  be  determined  statistically.  That  is,  time-averaged  values  of  partial- 
density  data  must  be  obtained  at  several  points  along  a  callbratioa  propagation 
path.  These  time -averaged  values  then  must  be  averaged  again  against  each 
other  to  obtain  what  amounts  to  a  time -and -space -averaged  partial -density  value. 
Multlplylag  this  latter  value  by  the  propagation  path  length  should  yield  something 
comparable  to  what  the  active-system  transmissometer  responds  to  on  a  time- 
avciage  basis. 

Figures  39a  and  b  are  examples  of  the  calibration  procedure  l  ist  * 
described.  The  solid  curve  in  Figure  39a  represents  the  plot  of  half-hourly 
values  of  line -integrated  oxygen  density  (obtained  with  the  Boeing  Infrared  system) 
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a.)  Half-Hourly  Values 


Figure  39:  COMPARISON  OF  INTEGRATED  OXYGEN  DENSITY  AND  DRY-Bt'LB 

temperature  data 


measured  ever  a  2-kllometer  path  the  nlg^t  of  2  November  1965.  The  dashed 
curve  deplete  the  corresponding  half>hourly  values  of  dry-bulb  temperature  re¬ 
corded  at  one  end  point  of  the  propagation  path.  The  abscissas  are  labeled  in 
arbitrary  units  because  the  data  are  presented  for  comparative  purposes  only. 

A  correlation  between  the  ^  -data  and  the  temperature  data  is  readily  apparent. 

Figure  39b  depicts  90-minute  running  means  of  the  data  contained 
in  Figure  39a.  Note  that  the  simple  time -averaging  process  has  considerably 
"smoothed*'  both  curves,  and  that  from  about  10:00  j«to  4:00  a.m.  there  appeared 
to  be  even  a  (me-to-one  correlatioo  between  the  two  sets  of  data.  The  conformance 
of  the  two  sets  of  data  would  have  been  even  more  noticeable  if  additional  tempera¬ 
ture  measurements  had  been  performed  at  other  points  alcmg  the  path.  The  spread 
in  data  points  for  a  regression  analysis  would  thus  have  been  minimized,  and  the 
reliability  of  the  ultimate  calibration  curve  correspondingly  increased. 

The  process  just  described  should  be  completed  for  vaiious  path 
lengths  under  a  variety  of  weather  conditions  to  establish  the  calibration  for  the 
entire  range  of  path  lengths  (and,  hence,  in  transmittance)  to  be  encountered  in 
missile  tracking  operations.  Note  that  for  calibration  purposes  fast -response 
meteorological  sensors  are  not  required  because  averages  in  both  time  and  space 
are  to  be  performed  anyway.  Measurement  accuracy  is  the  desirable  instrumenta¬ 
tion  feature.  Partial  deasities  are  not  to  be  measured  directly,  but  are  to  be  in¬ 
ferred  (via  the  dictates  of  the  equations  of  state)  from  precision  measurements  of 
dry-bulb  and  dewpoint  temperatures  plus  atmospheric  pressure.  Boeing  recom¬ 
mends  that  at  least  one  set  of  calibration  measurement  be  performed  over  a  slant 
path  from  a  valley  to  a  mountain  top  in  an  attempt  to  determine  the  line -broadening 
pressure  correctian. 

(2)  Passive -system  calibration 

Calibratloo  procedures  for  the  passive  systems  would  be  similar  to 
those  for  the  acb.'e  systems.  That  is,  calibration  would  be  based  on  a  compari¬ 
son  of  radiometric  data  widi  radiosonde  tempeiature,  pressure,  and  humidity 
data.  Radiosonde  data,  however,  represent  app  roximately  the  amounts  of  water 
vapor  and  oxygen  in  a  vertical  column  of  the  atmosphere.  Hence,  to  establish 


the  passive  system  slant -path  calibrajloa ,  the  radiosonde  data  would  be  multiplied 
by  the  appropriate  air -mass  factor.  Again,  the  final  calibration  cunre  would  be 
established  on  a  statistical  basis,  with  a  regression  analysis  performed  oi^  a  large 
mass  of  acquired  radiosonde  and  radiometric  data. 

B.  Engineering-model  evaluation 

(1)  Active -system  evaluation 

Ihe  active -system  calibration  procedure  outlined  above  will  do  nothing 
to  establish  the  capability  of  the  active-system  line-integral  refractometer  for 
effecting  the  real-time,  short-term,  range-error  corrections.  Such  an  evaluation 
can  best  be  performed  through  actual  ranging  measurements.  That  is ,  assume 
that  a  radar  ranging  system  was  set  up  over  a  fixt'd,  measured  course ,  such  as 
between  geodetic  bench  marks.  The  electilcal  range  measured  by  the  radar 
would  shon’  short-term  fluctuation  errors.  Then,  if  the  calibrated  active  systems 
were  set  up  parallel  to  the  radar  path,  the  actlve-s^tem  data  could  be  used  to 
compute  N  for  the  path  and  to  correct  the  radar  data  continuously,  If  the  meas¬ 
ured  electrical  ranges  were  corrected  to  the  known  range  continuously  and  to 
within  the  precision  limits  of  the  radar,  then  the  refractometer  evaluation  would 
be  established.  If  a  residual  bias  error  remained,  an  error  in  the  refractometer 
calibration  could  be  expected.  A  random  residual  error  would  indicate  a  '*nolse" 
error  in  the  refractometer  system.  Simultaneous  meteorological  measurements 
would  not  be  required  in  th*"  evaluation  scheme.  Again,  such  tests  should  be 
performed  over  a  range  of  prupagatlan-path  lengths. 

(2)  Passive -system  evaluation 

Because  the  passive  system  monitors  only  the  self-emissioo  of  the 
sky,  the  electrical  range  measurements  to  bo  corrected  would  Imve  to  be  meas¬ 
ured  to  an  airborne  target.  High-performance  aircraft  would  suffice;  that  is, 
electrical  range  m^isurements  would  be  made  to  the  target,  and  these  nneasure- 
ments  would  be  corrected  with  the  passive -system  refractometer  data.  Again, 
if  the  correction  yielded  true  slant-range  data  that  was  correct  to  within  the  pre¬ 
cision  limits  id  the  radar,  the  capability  of  the  refractometer  would  be  established. 
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A  LONO-PATH  ABSORPTION  REFRACTOMETER 


W.  G.  Tank* 


INTRODUCTION 

Because  the  earth's  atmosphere  Is  generally  turbulent,  the  density  and 
water  vapor  content  of  air  vary,  at  every  point,  with  time.  This  heterogeneity 
of  atmospheric  composition  implies  that  the  speed  of  electromagnetic  signals 
throqgh  the  atmosphere  la  a  function  of  location.  Mathematically,  this  fact  can 
be  eiqiTessed  as  n  =  n(x),  where  n  deiKJtcs  the  normalized  propagation  speed, 
or  atmospheric  refractive  index.  The  conversion  factor  from  time  measitre- 
ment  in  systems  of  space  measurement  based,  as  are  electromagnetic  ranging 
devices,  essentially  on  the  simple  expression  that  distance  is  the  product  of 
speed  and  time,  must  therefore  be  the  space -averaged  propagation  speed. 


n  =  y  n(x)dx  , 
o 

where  L  denotes  the  propa^tion  path.  In  determining  n  with  standard 
meteorological  instrumentation,  one  is  always  limited,  physically,  to  determin 
ing  n<x)  at  some  finite  number  of  values  of  x  along  the  path.  The  integral 
expression  for  n  permits  the  interpretation  that,  even  though  errorless 
determinations  of  the  index  n  at  each  of,  say,  k  sampling  points  were  made, 
the  average  index  n  would  still  be  subjject  to  a  finite  sampling  error  to  the 
extent  that 
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L 
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n(x)dlx. 
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The  above  inequality  constitutea  the  current  limit  to  the  accuracy  attainable  in 
distance  measurement  using  electromagnetic  ranging  techniques.  This  paper 
proposes  a  refractive  index  measuring  system  under  development  at  Boeing 
capable  of  relaxing  that  limit. 


THE  INSTRUMENT  SYSTEM 
Theoretical  Background 


The  refractive  index  at  a  point  in  clear  air  can  be  written  as  the  sum  of 
two  terms, 


N 

n  =  A^Pl  +  Bp^. 

i=l 


(2) 


where  the  p's  refer  to  the  densities  or  concentrations  of  the  different  atmos¬ 
pheric  gases.  Specifically,  refers  to  the  water  vapor  density,  p.  to  the 
densities  of  the  other  constituents  (oxygen,  nitrogen,  argon,  etc. ).  The 
indicated  summation  is  to  proceed  over  all  these  latter  constituents,  thus 
providing,  according  to  Dalton's  Law,  a  measure  of  the  total  dry-air  density. 
Water  vapor  is  by  far  the  most  variable  constituent  of  the  atmosphere.  The 
dry  air  components  show  remarkable  constancy  in  their  relative  proportions. 
Hence,  knowledge  of  the  air  content  (i.  e. .  amount  per  unit  volume)  or  density 
of  Just  one  of  these  latter  constituents  constitutes  such  knowledge  of  all  the  rest. 

Accepting  Equation  (2)  and  ail  its  implications,  the  absorption  of  infrared 
radiation  by  different  atmospheric  gases  at  selected  wavelength  bands  now 
provides  the  basis  for  an  instrument  system  to  determine  n  as  a  truly  integral 
function  of  path  length,  for  the  total  amount  of  energy  absorbed  from  a  beam 
propagated  throu^  the  atmosphere  dc|)ends  only  on  the  amount  of  the  absorbing 
constituent  along  the  path,  and  not  on  its  distribution.  Hence,  monitoring  an 
Infrared  beam  in  a  water  vapor  absorption  band  and  In  an  absorption  band  of  one 
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of  the  dry  air  constitucnta  will  provide  all  the  information  required  to  compute 
refractive  indices  according  to  Equation  (2),  which  value,  because  of  the  actual 
physical  process  of  absorption,  will  then  precisely  represent  the  fl  defined  In 
Equation  (1). 

System  Description 

The  idea  of  determining  atmospheric  composition  by  measuring  atmospheric 
Infrared  absorptions  la  certainly  not  new,  the  idea  preceding  efficient  means 
of  implementation  by  nearly  forty  years  (Fowle^^^  ).  The  concept  is  simplicity 
itself,  requiring,  basically,  a  continuous  spectrum  radiation  source,  a  i.',(iiation 
detector,  and  means  for  effecting  narrow-band,  wavelength  isolation.  The  need 
for  a  precisely  calibrated  source  can  be  negated  by  measuring  received  energy 
in  two  different  but  adjacent  bands  of  the  infrared  spectrum,  one  of  which  is 
attenuated  1^  the  constituent  of  interest,  the  other  not.  The  desired  constituent 
density  measurement  is  then  provided  by  the  ratio  of  the  two  band  energies 
sensed,  bt  this  scheme  the  former  band  Is  considered  the  sensor,  the  latter 
the  reference.  To  determine  the  atmospheric  refractive  index  in  this  manner, 
then,  a  minimum  of  three  wavelength  bands  must  be  monitored  —  two  serving 
as  sensors  (one  for  water  vapor,  the  other  for  any  one  of  the  dry-air  gases), 
the  third  as  the  reference. 

The  Sensing  and  Reference  Bands 

The  choice  as  to  which  three  wavelength  bands  to  use  is  largely  arbitrary. 
The  three  should,  however,  be  grouped  as  closely  as  possible  together  in  order 
to  minimize  adverse  effects  on  system  accuracy  that  could  be  caused  by  differ¬ 
ential  scattering  effects  in  t  hazy  atmosphere  (more  on  this  later).  Furthermore, 
to  take  advantage  of  the  simplicity,  reliability,  and  obtainabillty  offered  by  lead 
sulfide  cell  Infrared  detectors  as  system  radiation  sensors,  all  wavelengths 
should  be  in  the  region  of  peak  sensitivity  of  such  cells.  Accordingly,  consider 
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Figure  1 ,  The  Near- Infrared  Solar  Absorption  Spectivm  (The  absorbing 
atmospheric  constituents  are  identified  at  the  top  of  the  chort 
above  the  different  absorption  bonds,  the  refroctometer  filter 
bands  by  the  three  dashed  abscissos.) 


Bolar  absorption  spectrum,  normalized  to  peak  transmission,  presented  by 

(2l 

Gates  '  and  shown  in  Figure  1.  The  abscissa  range  indicated  defines  the 
wavelength  region  of  peak  sensitivity  of  uncooled  lead  sulfide  cells  (Kruse, 
McGlauchlin  and  McQuistan 

The  spectrum  shown  in  Figure  1  is  characterized  by  four  strong  water 
vapor  absorption  bands,  four  "window"  regions,  a  weak  oxygen  and  a  weak 
carbon  dioxide  band.  Radiation  on  any  of  the  window  regions  could  be  moni¬ 
tored  to  supply  the  reference  signal ,  that  in  any  of  the  water  vapor  absorption 
bands  t>  effect  the  measurement  of  water  vapor  content,  and  that  in  the  oxygen  or 
the  carbon  dioxide  band  to  determine  the  dry -air  constituent  content. 
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Three  narrow  wavelength  banda  aatiafying  the  requirenents  of  an  absorp¬ 
tion  refractometer  are  thua  specified  as  indicated,  centered  at  wavelengths  of 
2. 12  p,  1.  73  p  and  1.  256  p  for  the  reference,  water  vapor  and  oxygen  or  dry- 
air  density  sampling  banda,  respectively.  The  water-vapor  sampling  band  is 
displaced  from  the  center  to  the  wing  of  the  strong  1.  86  p  absorption  band  for 
two  reasons.  First,  over  apy  sort  of  range  at  all.  virtually  all  the  1.  86  p 
radiation  would  be  absorbed  from  an  infrared  beam,  leaving  no  signal  to  effect 
the  desired  water  vapor  measurement.  Second,  energy  in  the  oxygen  band  suf¬ 
fers  some  water  vapor  absorption  by  the  wing  of  the  1. 37  p  water  vapor 
absorption  band.  Water  vapor  absorption  at  1. 73  p  is  equal  to  that  in  the 
o^QTgen  band.  Hence,  by  referencing  the  band  energy  of  the  latter  to  that  of  the 
former,  unwanted  water  vapor  absorptions  at  1.  256  p  can  be  negated.  The 
1. 73  p  water  vapor  sampling  band  itself  is  referenced  to  the  window  band. 

Receiver  Optica 

An  optical  ayatem  for  simultaneously  monitoring  a  beam  of  infrared 
radiation  at  the  three  wavelengths  just  specified  through  a  single  aperture  is 
shown  schenutically  in  Figure  2.  A  mechanically  chopped,  collimated  l)cam  of 
energy  from  a  remote,  broadband,  tungsten  filament  source  of  high  thermal 
inertia  is  directed  at  the  receiver.  The  Cassegrain  collecting  optics  focus  the 
Incoming  beam  onto  a  field-stop  aperture.  The  divergent  beam  then  piasscs 
through  a  series  of  mirror  beam  splitters,  the  itrst  of  which  has  a  transmlssiun- 
to-reflectlon  ratio  of  about  2;1,  the  second  a  ratio  of  1: 1.  The  beam-splitter 
series  thus  serves  to  split  the  collected  beam  into  three  components,  directing 
approximately  one-third  of  the  collected  energv  onto  one  of  the  three  indicated 
narrow  bandpass,  germanium  interference  filters.  Each  arm  of  the  split  and 
filtered  beam  then  impinges  on  a  mixer  cone  which  serves  to  distribute  the 
filtered  energy  evenly  over  the  surface  of  a  lead  sulfide  cell  Infrared  detector. 
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Figura  2.  Schematic  Diogram  of  the  Refroctometer  Receiver  Bead  Optici 

This  even  distribution  eiiminates  shot  noise  effects  that  could  be  caused  by  in¬ 
equalities  in  sensitivity  at  various  spots  on  the  detector  surfaces.  The  interior 
reflecting  surfaces  of  the  mixer  cones  are  gold-plsted  for  the  greatest  degree 
of  reflection  efficiency.  Source-to-receiver  alignment  is  accomplished  with  the 
aid  of  20X  spotting  telescopes  mounted  on  both  the  beam  and  field  collimating 
optics.  The  entire  receiver  head  is  attached  to  a  standard  3-1/2-  x  8-inch 
thread  transit  head  mount  for  utilization  on  a  standard  surveyor's  transit  tripod. 

The  Read-t  .t  System 

The  desired  comparison  between  reference  and  sensing  band  energies 
is  their  ratio,  because  such  a  comparison  reflects  only  disproportionate,  and 
not  proportionate,  energy  changes.  That  is.  system  response  would  not  re- 


reflect  a  change  in  the  general  level  of  ecnaed  energy  caused  by  the  presence  of 
haze  or  fog  in  the  light  beam.  The  desired  ratio  is  obtained  as  follows. 

Located  in  the  receiver  head,  directly  coupled  to  the  lead  sulfide  cell 
detectors,  are  high  input,  low  output  impedance  detector  preamplifiers  (Figure 
2).  The  amplified,  essentially  alternating-current  (as  a  result  of  the  source 
chopping)  detector  signals  are  then  fed  into  alternating-current,  capacitor- 
coupled  (to  eliminate  the  high  direct-current  cell  bias)  logarithmic  converters. 
These  converters  rectify  Ute  alternating-current  signals  and  give  direct-current 
outputs  proportional  to  the  logarithm  of  the  root  mean  square  value  of  the  input 
signals.  The  c'^nverter  outputs  are  then  combined  in  a  differentiating  network 
to  yield  a  recorded  signal  prt^rtional  to  the  difference  between  the  logarithms, 
and  hence  to  the  ratio,  of  the  detector  input  signals.  Schematically,  the  read¬ 
out  system  is  shown  in  Figure  3. 


FkQurs  3.  Theoreticol  Rsfroctometer 
Rsodout  Sytftn 
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ANTICIPATED  OPERATIONAL  CHARACTERISTICS 


The  inetrumenl  Juat  described  is  atlll  in  the  development  stage.  But 

experience  with  an  infrared  absorption  hygrometer  previously  developed  and 

[41 

applied  (Tank  and  Wergin  '  ).  and  operating  according  to  the  philosophy  and 
methods  of  the  proposed  refractometer.  permits  the  following  refractometer 
performance  characteristic  estimates. 

Range  Capabilities 

An  estimate  of  the  range  capability  of  the  proponed  refractometer  boils 
down,  essentially,  to  the  answer  to  the  following  question:  how  far  can  th£ 
source  be  back<.>'i  off  frem  the  receiver  such  that  the  radiant  energy  collected  is 
still  sufficient  to  generate  a  signal-to-noise  ratio  significantly  greater  than 
unity?  The  power  received  from  a  remote  source  emitting  a  directed  beam  of 
energy  in  the  wavelength  Interval  AA  is 
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(3) 


where 

I^  =  intrinsic  spectral  radiant  omittance  of  the  source, 
watts  cm"3  micron'* 

A  =  the  area  of  the  collecting  optics  aperture; 

=  bandwidth  of  the  wavelength  isolation  filters,  microns 
=  transmittance  of  the  optics; 

0 

r  -  transmittance  of  the  atmosphere,  and 

A 

L'  propagation  path  length. 

The  solid  angles  W  and  W  arc  defined  by  the  source  collimating  o|)tics. 

1  ^ 

Specifically,  is  the  solid  angle  subtended  by  the  collinuiting  optics  aperture 
at  the  source,  the  solid  angle  subtended  by  the  source  at  the  aperture  The 


four  leading  tarma  on  the  right  aide  of  Equation  (4)  ar^  thus  simply  system 
parameters.  Both  the  beam  and  field  collimators  used  in  the  refractometer 
consist  of  8-inch  aperture,  128-inch  (folded)  focal  length  Caesegrainion  re¬ 
flectors.  The  filter  bandwldths  are  0. 025  micron.  For  a  half-inch  source 
diameter,  then.  Equation  (3)  can  be  written,  utilizing  the  numerical  values  of 
the  system  parameters  and  taking  as  a  pessimistic  estimate  of  a  value  of 
0. OS,  as 


where  L  is  expressed  in  kilometers. 

When  the  source  emittance,  1^ ,  is  expressed  as  defined  in  Equation  (3), 
the  power  received,  W  ,  is  given  in  watts.  That  incident  radiation  on  an  un- 

A 

cooled  lead  sulfide  cell  which  produces  a  signal-to-noise  ratio  of  unity  is,  in 
the  wavelength  region  of  interest,  of  the  order  10  watts  (Kruse,  McGlauchlin 
and  MoQuistan^^^  ).  Therefore,  for  workable  signals  the  following  inequality 
must  hold: 


l.S  X  10 


■X  ? 


The  emlttanoe  of  a  1, 400^C  blackbody  (peak  emission  at  about  1.7  m)  will  be  of 

-2  -I 

the  order  of  10  watts  cm  micron  .  Utilizing  this  value  for  the  source  emlt¬ 
tanoe,  inequality  (6)  may  be  expressed  as 


The  range  at  which  the  refractometer  can  be  worked,  then,  depends  on  the 
atmo^>hsric  transmittance.  For  ranges  on  the  order  of  tens  of  kllometera. 
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the  transmitUince  may  (pesslmiBtically)  be  expected  to  be  on  the  order  of  10"^. 
’i'luH  estimate  considers  both  absorption  and  scattering  losses  in  tlic  pru|)agutcd 
1‘cam.  Hence,  the  workable  range  of  the  system,  too,  will  be  on  the  order  of 
tens  of  kilometers. 


jcled  Calibration 


It  is  expected  that  the  calibration  of  the  rcfractomctcr  will  prove  to  be 

[2l 

quite  predictable.  Gates  has  shown  that  throughout  the  spectral  region 
0.  67  p  to  2.  54  p  both  the  random  line  (water  vapor)  and  regular  line  (oxygen) 
absorption  band  transmissions  vary  as  the  square  root  of  the  total  amounts  of 
absorbing  material  in  the  propagation  path  according  to 


InT  =  Cm^^^  (7) 

In  the  above  Equation  (valid  for  transmissions  ranging  between  20  and  80 

percent),  T  represents  transmittance,  C  a  transmission  coefficient,  and  m  the 

total  amount  of  absorbing  constituent.  For  water  vapor,  this  latter  quantity 

is  usually  expressed  as  precipi table  centimeters,  w,  of  water  and  for  oxygen 

as  atmo-km,  i,  of  oxygen.  The  recorded  signals  from  the  refractomcior  may 

1/2 

thus  be  expected  to  vary  linearly  with  m  ,  as  past  work  with  the  infrared 

(4] 

hygrometer  indeed  verifies  (Tank  and  Wergin  ).  Over  the  range  of  w  and 
1  values  to  be  expected,  and  for  propagation  paths  ranging  in  length  from  5 
to  5C  kilometers,  the  !  and  w  calibration  erves  for  the  refractometer  are  as 
shown  in  Figure  4.  The  cur  es  are  calculated  as  dictated  by  Equation  (3)  and 
by  the  response  characteristics  of  the  previously  calibrated  infrared  hygro¬ 
meter.  The  dry-air  and  water  vapor  densities  can  be  recovered  from  the 
measured  parameters  according  to  the  following  relationship: 
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vnsy  (UffER  LANIlS)  AND  ttOMA  LABELS) 


X 

iil 

Figui*  4.  ThMfVttcal  lt«fracfQin«l»r  Colibfolton  CutvM 

where  and  p^  denote  the  dxy-air  and  water  vapor  deneitiee,  respectively, 
in  the  unite  shown.  When  both  L  and  w  are  given  in  kilometers  and  centi¬ 
meters,  respectively,  the  propapation  path  length,  L,  must  be  given  in 
kilometers. 

Sxpected  Accuracy 

In  estimating  the  accuracy  of  dry-air  and  water  vapor  measurement 
achievable  with  the  refractometer,  the  following  factors  must  bo  considered: 
electronic  noise  effects,  ambient  temperature  and  pressure  effects  on  ab¬ 
sorption,  differential  scattering  attenuation  effects,  aiui  the  variability  of 
atmospheric  carbon  dioxide  content.  The  manner  in  which  these  factors  act 


to  limit  lystcm  accuracy,  and  the  expected  magnitude  of  the  effect, 
can  be  presented  as  follows. 

Electronic  Noise  Uvfl  Umits 

Amplifier  and  recorder  noise  combine  to  give  a  syatem  noise  level  of 

something  less  than  10  microvolts.  Based  on  the  expected  system  respon- 

sivity  as  indicated  in  Figure  4,  such  a  signal  level  is  equivalent  to  optical 

•7  *6 

depths  of  o)(ygen  and  water  vapor  of  about  10  kilometers  and  10  centimeters, 
respectively.  For  ranges  of  orders  of  kilomoters,  such  magnitudes  imply 
electronic  noise  errors  in  oj^gen  and  water  vapor  densities  of  orders  one  part 
pep  hundred  million  and  one  part  per  billion,  respectively.  In  other  words,  the 
effect  is  negligible. 

Ambient  Temperature  and  Pressure  Effects 

The  absorption  of  energy  within  a  certain  wavelength  band  is  dependent  to 

some  extent  on  both  atmospheric  temperature  and  pressure.  Although  such 

effects  become  appreciable  for  drastic  changes  in  ambient  temperature  and 

pressure,  calibration  curves  established  at  given  pressure-altitudes  should 

remain  insensitive  to  those  changes  associated  with  the  normal  progression  of 

"weather.  "  This  statement  is  substantiated  by  past  work  with  the  Infrared 

f47l 

hygrometer  (Tank  and  Wergin*  ). 

DifferenUal  Scattering  Attenuation  Effects 

Scattering  of  near  infrared  energy  by  atmospheric  haze  particles  is  wave¬ 
length-dependent.  Specifically,  the  scattering  theory  predicts  that,  for  given 
particle  size  distributions,  the  shorter  the  wavelength  th"  more  efficient  the 
scatterer.  Consequently,  the  radiation  intensify  perceived  at  1. 73  ft  will  be 
disproportionately  less  than  that  perceived  at  2. 12  p,  and  that  at  1. 256  p  will 
be  less  than  that  at  1, 73  p  .  Baaed  bn  data  on  light  transmission  through  hazb 
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presented  by  Gibbons/  the  magnitude  of  the  effect  of  the  refractometer  out¬ 
put  signals  Is  as  shown  In  Figure  S.  Both  sets  of  ordinate  labels  are  normal Ir.ed 


Ordinate  labels  ore  normolized  to  propogotion  poth  length  in  kilometers.  The  upper 
set  of  righMiond  ordinate  labels  refers  to  atmo-kms  of  oxygen;  the  lower  set  to  pre~ 
cipitoble  ems  of  water. 

Figure  5.  Calculated  Error  Signols,  Couted  by  Differential  Scattering  Effects  in 
Atmospheric  Haze,  as  Functiom  of  Visual  Ronge 


to  the  propagation  path  length.  L.  in  kilometers.  The  corrections  implied 
are  negative  corrections,  insofar  as  the  equivalent  optical  depth  values  express 
apparent  excesses  of  oxygen  and  water  vapor  contents  attributable  to  differen¬ 
tial  scattering  effects.  But  again  the  effect  is  negligible,  because  for  |>ath 
lengths  of  tho  order  of  10  kilometers  the  scattering  noise  signals  arc  comparuble 
to  those  associated  with  electronic  noise. 
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The  VhrlebUity  of  Atmospheric  CO^  Content 

Dslton’s  law  relsting  total  pressure  to  partial  pressure  can  be  cx|)rcssod 
In  terms  of  density  as 

■■  ■'’n)- 

That  is,  total  dry-air  density  can  be  considered  the  sum  of  the  individual  at¬ 
mospheric  constituent  densities.  Individual  constituent  densities  can  be 
expressed  as 

Pj  =  MRjP,j,  .  (11) 

th 

where  M  denotes  the  1  constituents  mixing  ratio,  i.  e. ,  its  |)erccnt  by 
weight  of  the  total  dry-air  mass  per  unit  volume.  At  this  point,  the  philosophy 
of  the  refractometer  is  that  all  dry-air  constituents  exhibit  consbint  mixing 
ratios.  Though  generally  true  for  the  so-called  permanent  atmospheric  gases 
(nitrogen,  oxygen,  argon,  plus  helium  and  hydrogen  and  the  noble  gases),  the 
atmospheric  carbon  dioxide  content  does  show  considerable  variation.  This 
variation  is  due  to  a  variety  of  reasons,  chief  of  which  are  the  release  of  CO^ 
to  the  atmosphere  by  combustion  processes  and  the  exchange  via  biological 
processes  of  CO  with  soil  and  vegetation.  At  any  rate,  the  CO  mixing 

£  if 

ratio  is  more  property  expressed  as 

""co  ^  ^““co  ' 

2  2  2 

where  the  overbar  denotes  a  constant  "average"  mixing  ratio  and  the  A  a 
variable  component  brou^t  about  by  the  processes  must  mentioned. 

Utilizing  Equatlona  (11)  and  (12)  in  accordance  with  Equation  (10),  then 

m 

^T  “  ^T  Z  M  rj  +  A  M  R^q  .  (13) 

i-l  ^ 
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The  flummation  In  Equation  (12)  ia  equal  to  unity,  and  the  remaining  term  on 
the  right  side  is  thus  Interpreted  as  an  error  in  the  inferred  density,  as 
measured  by  the  refraetometer,  aeeribable  to  a  variable  atmospheric  CO 

£ 

contmt.  The  variation  in  CO2  mixing  ratio  con  attain  a  value  of  600  parts 
per  million.  Hence,  refraetometer  measured  densities  will  be  correct  to  this 
eame  ratio.  This,  then,  is  the  fundamental  limit  to  the  accuracy  of  density 
measurements  attainable  with  the  refraetometer. 

It  should  be  remarked  that  the  above  limit  does  not  apply  to  the  water 
vapor  density  measurement,  insofar  as  this  latter  quantity  is  measured  direct¬ 
ly.  Stystem  sensitivity  will  permit  water  vapor  densities  to  be  determined  with 
accuracies  of  the  order  of  one  part  per  million., 

CONCLUSIONS 

The  claims  made  regardiaf  the  performance  of  the  infrared  absorption 
reftaetomeler  Just  described  remain  to  be  proved,  insofar  as  the  first  proto¬ 
type  unit  is  still  to  be  tested  la  the  field.  But  in  light  of  past  experience  with 
an  absorption  hygroaaeter,  the  claim  as  to  the  ability  through  use  of  the  re- 
fraotometer  to  measure  dry-air  aad  water  vapor  densities  to  within  six  parts 
per  tea  thousand  aad  a  few  parts  per  million,  respectively,  is  not  considered 
unreasonable.  Such  accuracies  imply  that  the  atmospheric  refracth'e  index 
may  be  specified  to  seven  significant  figures.  And  what  is  Just  as  important 
relative  to  use  of  the  refraetometer  ia  conjunction  with  electromagnetic  ranging 
devloes  la  that  the  refractive  index  so  specified  will  represent  a  apace  averaged 
value  determined  as  a  truly  integral  function  of  path  length. 


163 


REFERENCES 


F.  E.  Fowle,  "The  Spectroscopic  Determination  of  Aqueous  Vapor," 
Astrophys.  J.  35,  149-162,  1912. 

D.  M.  Gates,  "Near  Infrared  Atmospheric  Transmission  to  Solar 
Radiation."  J.  Opt.  Soc.  Am.,  1299-1304,  1960. 

P.  W.  Kruse,  L.  D.  McGlauchlin  and  R.  B.  McQuistan,  Elements  of 
Infrared  Technology:  Generation,  Transmission,  and  Deteetlon, 

New  York.  John  Wiley  and  Sons,  Inc. , 

W.  G.  Tank,  and  E.  J.  Wergln,  "A  Long-Path  Infrared  Hygrometer, " 
in  Humidity  and  Moisture.  Vol.  I.  A.  Wexler,  ed.,  Nejy  York, 

Relnhold  Publishing  Co. ,  1964. 

M.  G.  Gibbons,  "Wavelength  Dependence  of  the  Scattering  Coefficient 
for  Infrared  Radiation  in  Natural  Haze,  "J.  Opt.  Soc.  Am. ,  172-176, 

1958. 


161* 


UNCLASSIFIED  - 


Security  ClaMification 


DOCUMENT  CONTROL  DATA  -  RAD 

ra«<wri>r  clmlllcMm  al  tula,  body  at  abalntt  and  tndaamt  mnalallan  mual  W  aniaiad  atian  lha  aaarall  rapart  la  claaaillad) 


t  oaic  >4*tlNC  ACTIVI'^V  (Caipaiaia  aulhar) 

Bc<iin);  Company 
S<?attle  Washington 


|2«  ACPOMT  tCCuniTV  C  LAfllFiCATlON 

Unclassified _ 

26  anou^ 

flA 


3  Rt  rOfIT  TITCI 


Improved  Sampling  Techniques 


I^AUTHOfKS;  rLMlfMM*.  Amiamm.  Mtil) 

Tank,  William  G.  Palmer,  Thomas  Y. 

Kreiss,  William  T. 

Rouley,  James  J. 


7«  .1.  MO  OR  RAOID 

188 

66 

A#.  OIIIOIM*TOf»>6  RCROIRT  MUM 

f  6-  ^J^KR^^ROMT  mo^5^  (Any  orfwr  aii6aw  Aa#  may  6«  aaal0tmd 

RADC-TP-66-I*07 

4  Oltcaierlvc  NOTIS  ITypa  at  tapan  and  laalualaa  daiaa) 

Final  Report 


•  aiaOMT  DATE 

August  1966 


l«  contract  or  orant  no. 

AF30(602)-3715 

A  RROJtCT  NO 

5579 

$'57902 

d 


10  A  VA  ILAalLITV/LIMITATION  NOTICU 

This  document  is  subject  to  special  export  controls  and  each  transmittal  to  foreign 
trovernments  or  foreign  nationals  may  be  made  only  vith  prior  approval  of  RADC(EyLI) 
GARB.  N.Y,  13)«1*0 


II  lusatsMeNTAav  motss  u.  (ponookino  uiutanv  activity 

Space  Surveillance  6  Instru.  Pranch 
Rome  Air  Development  Center 

_ Griffiss  APR  Mow  York  13440 _ 

IJ  ASSTRACT 

Line-integrated  rcfractivity  is  isolated  a--  the  basic  atmospheric  I 
parameter  to  be  determined  for  correctiont.  0  radar  ranging  data. 

A  basis  for  the  determination  of  integrated  rcfractivity  from 
measurements  of  the  absorption  of  infrared  and  microwave  radiation 
by  atmospheric  oxygen  and  water  vapor  is  then  established.  Three 
absorption-measuring  systems  are  next  described,  and  engineering-model 
design  specifications  are  presented.  Each  of  the  specifications 
offers  particular  advantages  to  particular  ranging  operations.  It 
is  estimated  that  incorporation  of  such  devices  into  ranging  systems 
will  permit  real-time  measurements  of  slant-range  and  elevation  angle 
to  accuracies  five  to  ten  times  those  currently  achieved. 


Vrtclo.r}sifie6 


IvccurUy 


Rsv  ecrrca 


f  Trcpcnphcric  nefriicticn 
M^cfr^ction  Correction  Tcchniqiics 


MOCK  (  nr 


ItOTCtJCVSJXJ 


1.  0SiG!WATl?®0  ACTlVTTYt  EaS«f  t!w  Gsa»e  f c:4A«s 
c!  (h«  C9SC7ftctct,  co^c«otr«cCcr.  er«e£e«.  D«poetaftra2  cf 
fsnM  activtly  or  oihmr  erscaiCMtioA  fc«rpor*<«  «uihof)  IkkcIsq 
(h«  report. 

:*.  REPOST  OTCUOTT  ClJiSES^ATOfl:  ots> 

fiKcttfliy  ciK»aifU«ilOA  oi  ttsm  rapctft*  la^UealK  «bet2i07 
''Ksauictoid  Dmk'*  Is  iadodvdl  Marklfl^  U  to  b«  lit  accord 
•ACS  «rith  oppfoprloeo  oocortty  roc«Ao)loos« 

CROUP:  AttfiMBiiEic  down^odinc  U  tpocUUd  ia  DcD  Di- 
rocUvo  9200. 10  ood  Aioaod  Pofcoo  Lsteotriol  ^lanornrii  fiater 
(ho  croi^  OfKrd>or.  Ateo,  whoo  otiptiCKblo*  t%9m  tliot  opftof  J 
surtOocs  hovo  b««o  moil  (m  Omp  S  oad  Oroop  4  m  •alliar>* 

iMd 

9.  RBPOStT  TmXi  Tjotm  C>»  cM^loto  rspott  to  «1S 
c^ttol  tottsro.  TlUm  to  i£l  coooo  itesSd  bo  aBaolooo&flo4 
ti  A  Mo&lacfxU  tl&lo  eoMot  bo  mloctod  wttlnol  cloooilleo* 
tion,  a£»w  ilUo  closoliUolioo  In  oil  ccpMolo  i«  poroaifcoolo 
lEcaeodioiAiy  foUowiap  Um  UUo* 

4.  DEXTCtPTlVE  HOm  tf  opproprioto*  ootSK  tSto  typo  oS 
r«pcf1,  o  &.«  Isoor&ai,  profyooo.  ooamory,  omool.  m  fln^ 

Ctvo  tha  (ocSukIto  dotos  orboo  o  v»oclfb:  rcporllsis  period  is 

cover  cd. 

5.  AtjriS3R:CS>  E#eor  tho  nmmeitt)  ol  csSboKo)  oo  clworii  eo 
or  lA  tho  report.  Boi«r  lost  ooMo,  ftraft  oobo,  bMA o  lomal. 

1/  sJlUory,  Alfesv  rordi  end  broncli  of  oorvtco*  Tto  mmo  of 
tho  prtncipAl  «)j|liof  U  oo  obooiolo  mloimaoi  rogilrononf 

6.  REPORT  DATD  Cotor  tbo  dofo  ef  tho  report  oo  doy, 
caoatlv  yoar,  or  noatK  7<oo  U  Boro  thoo  ooo  doto  oppoofo 
on  (ho  r«por(«  uoo  dots  of  pohlkotloo. 

7a.  total  tnnmn  op  PAOSIi  Tho  tstol  poco  cmas 

should  fellow  Qocuoi  po^Mtloa  prorodoros.  Lo..  oator  (ho 
nutobw  of  pspoo  cewtohiTt  Iwfcrwdeo, 

7b.  NUUUSR  07  REPERKNCFR  Enter  tho  tatsl  om^ror  of 
TcfcrcncAS  ettod  la  Iho  report. 

8a.  contract  or  ORART  NUK3CR:  If  spproprlcto.  oatev 
the  eppticAhl^BAnhor  of  (ho  cootroci  ce  roidof  which 
tho  ropoft  woo  wrtttOB 

ttJ,  ac,  til  &d.  PROJECT  NUlcaER:  Eolor  (ho  •^propdoto 
militory  dopodaool  IdmUflcoUoitf  oocb  so  proioct  oo^er, 
suhprojoct  QBtihor,  oyoloa  n—bero^  (aoh  miiobsi,  otix 
9«.  ORIOIRATOR'S  REPORT  IfUMBCRfS):  EcSer  tbo  ofO- 
clsl  report  nustbm  by  which  tho  doeuwont  will  bo  idoBtfiod 
And  controlled  by  (he  ortslnotuiQ  ortlvlty.  Tbit  atmbor  moot 
be  unique  lo  (hit  report* 

9t.  0Tt^E^!  REPORT  RUtSOEX^^:  U  Iho  report  boo  bo«a 
easicnod  sey  oCior  report  (HoObere  feirTior  by  (ho  er'^^nseor 
cr  by  (^9  tpensKft),  Also  cstor  ihio  nonborle). 
la  AVAtLAfSa-lTY/LOOTATlOW  KOTtCEa  Eater  ocy  Us> 
ttstiona  ce  fatbor  dUsocilnaUoa  of  tho  report,  other  thoo  thnccij 
B - 


IspoAod  by  cedT^^lly  ciaselTlcadao.  uAinp  otoodard  otatosiORto 
each  Ac: 


■y  chtAlA  coploo  o»  this 


(1)  **(^Aiif!od  foqmotora  i 
report  firoea  DDC.'* 

(2)  .**7oralsp  oiseo  one  urine  cat  oad  dloooBifkAtlcn  of  this 
report  by  DDC  U  M  ewtboriood*' 

(3)  **IX  IL  Oo^rOnMSOflt  OCOOCiOO  BAjr  ohtoia  coploo  of 
thio  report  dirocUy  hoB  DOC.  Othor  qaallflod  DDC 

BOOT*  ohoAl  lO^OOt  thlOBfh 


(4)  **19.  E  odlUsey  o«ooeloo  Boy  oUsio  coploo  of  thlo 

roBpeet  dkwctly  Rob  DDC  Othor  poollflod  oooro 
Ahoil  foqooot  throwgh 


(^)  **  All  diotrlheUoa  of  this  report  U  oootrolloA  Qub3* 
Uiod  DDC  Booro  oh«U  loqpstot  throeph 

tf  tho  report  heo  booa  fwraiohod  to  tbo  OfQco  of  Tochnlcsl 
Cerrlcoo,  DepBtBoB  of  CoBiaorco.  for  ooio  to  the  pObllc.  tndl^ 
cMo  Udo  foot  oad  owtor  tha  prtea.  If  kaowB 

IL  SVPPLDiEirrART  IdOTEBc  Uoo  for  odd  Uonol  osplAao- 
tory  ootOB 

IX  BPOlCaCMaieO  MnJTART  ACmTTY:  Eotortho  noBoof 
tho  dR^ortOBaUl  pcoloct  offIcB  or  laborotory  opOBOorlof  'par* 
fad  for)  tho  roooArch  oad  doeotopBoat.  Inclodo  oddroso. 

15.  ABSTRACT:  fotor  on  obolroct  flrlac  s  brlof  and  foctool 
OBOBory  of  tt»o  doewBoat  ladicolieo  of  Iho  report,  oroa  thouch 
It  aoy  olse  oppoor  oloowhoro  U  Iho  body  of  tbo  tochcicol  ro* 

U  AddlUoaol  Apoco  lo  ro<|alred.  o  contlnuAtioo  shoot  oKoll 


It  lo  hlchty  doolrsblo  thol  tho  obetroct  of  ctaoslflod  reports 
be  aoclsooiflod.  Bach  poroprspb  of  the  obetroct  oboll  ood  with 
ca  ladicotioo  of  (be  Bdlitary  oeewrity  clsoolflcatioo  cf  tbo  In- 
foTBcUoa  la  tha  parofyoi*.  repreoaated  oo  fTf).  ($}.  fC).  or  (U) 

TWre  lo  aa  Uatllatloa  oa  (ho  loacth  of  tho  obotroct  How- 
oeer.  tho  oappestod  loaptb  U  hoB  ISO  to  22S  srovdo. 

14.  EEY  WORDS:  Soy  wordo  ora  tochaleoUy  BooaiapfjS  toroui 
cr  ohort  phroooo  that  choroctortoo  o  report  oi^  Bay  bo  Bood  oo 
IndoK  oatvloo  for  cstoloplap  the  report.  Soy  words  must  be 
selected  oo  dtot  ao  oecarlty  eUooillcotioo  lo  reowlred.  IdooU' 
lUfB.  oach  BO  oqiiffwrsrnt  SMdsI  dooi^otloa.  trade  aoiso,  saiitory 
project  code  aos^a.  pBoyr^hlc  locotioa.  isoy  bo  mod  ao  key 
words  bat  wlU  be  followod  by  oa  iodicotioa  of  tochatcol  coo- 
loot.  The  Aoal^Bfootit  of  Uaks.  ruloo,  oad  weights  is  optlonel 


Unclaocif led 


S«^>ity  CluuiKicatioo 


